The 


Philosophical Magazine 


FIRST PUBLISHED IN 1798 
A Journal of Theoretical 


Experimental and Applied Physies 


Vol. 6 October 1961 No. 70- 
Eighth Series 


25s. Od., plus postage 
Annual Subscription £13 10s. 0d., payable in advance 


Printed and Published by 


AYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


THE PHILOSOPHICAL MAGAZINE 


Editor 
Professor N. F. Mort, M.A., D.Sc., F.R.S 


Editorial Board 


Sir Lawrence Brace, O.B.E., M.C., M.A., D.Sc., F.R.S. 
Sir Grorcr Tuomson, M.A., D.Sc., F.R.S. 
W. H. Taytor, M.A., D.So. 


AvuTHORS wishing to submit papers for publication in the Journal should 
send manuscripts directly to the Publishers. 


Manuscripts should be typed in double spacing on one side of quarto 
(8x10in.) paper, and authors are urged to aim at absolute clarity of 
meaning and an attractive presentation of their texts. 


References should be listed at the end in alphabetical order of authors and 
should be cited in the text in terms of author’s name and date. Diagrams 
should normally be in Indian ink on white card, with lettering in soft pencil, 
the captions being typed on a separate sheet. 


A leaflet giving detailed instructions to authors on the preparation of papers 
is available on request from the Publishers. 


Authors are entitled to receive 25 offprints of a paper in the Journal free 
of charge, and additional offprints can be obtained from the Publishers. 


The Philosophical Magazine and its companion journal, Advances in Physics, 
will accept papers for publication in experimental and theoretical physics. 
The Philosophical Magazine publishes contributions describing new results, 
letters to the editor and book reviews. Advances in Physics publishes articles 
surveying the present state of knowledge in any branch of the science in which 
recent progress has been made. The editors welcome contributions from 
overseas as well as from the United Kingdom, and papers may be published 
in English, French and German. 


| RO | 


A Note on Thermoelectric Power and Inelastic Scattering 


By A. M. Guinautty and D. K. C. MacDonatp 
Division of Pure Physics, National Research Council, Ottawa, Canada 


[Received March 18, 1961] 


ABSTRACT 


Recently Kasuya, Bailyn and de Vroomen have independently proposed 
that spin-dependent scattering of electrons by magnetic impurity ions would 
provide a mechanism to account for the ‘giant’ thermoelectric powers 
observed in certain dilute alloys at very low temperatures. A simplified 
model discussed here indicates directly that the essential requirement 
is that of two largely independent groups of conduction electrons whose 
mobilities are significantly different, but in which the groups are coupled by 
(spin-dependent) inelastic scattering. It also follows in general agreement 
with experiment that resistive anomalies are to be looked for when such 
“mutual electron-drag’ is present. 


A starting point for the solution of electron transport problems is the 
so-called ‘Boltzmann equation’ for the electron distribution function, f; 
and the Boltzmann equation may be written formally : 


(ir) ante, * (Fa, =° frat sae 


If the electron scattering is purely elastic, such as by lattice vibrations at 
‘high’ temperatures (k7' =hv) or by impurity ions whose scattering may be 
represented by a simple electrostatic potential, then it is argued on a quasi- 
free electron model that the collision term (which is generally an integral 
over all possible scattering processes (k, k’)) may be written : 


Dx, _ _ fuel #) ~ fol) 
(hy ce HORI) ere 
= 
where f,(Z) denotes the distribution function in equilibrium and # denotes" 
the energy of an electron having wave vectork. The ‘relaxation time’ 7 is 

a direct function of the scattering matrix elements appropriate to the 

collisions, but does not involve the Pauli exclusion principle (cf. Mott and 

Jones 1936, p. 260, Peierls 1955, p. 116). This means that any energy- 

dependence of + would arise directly from an energy-dependence of the 

matrix element for an electron making the transition (k, k’) in the absence 

of any other electrons; we shall assume in what follows that the matrix 

elements themselves are constant or at most relatively slowly varying 

functions of the total electron energy, E. 
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Quite generally, the absolute thermoelectric power, S, may be written - 
pe 1 of 0 of 0 3 
S= 7 | B-Do(h) Fear | | o(H) Sede eee (Say 


where, to quote Ziman (1959), o(Z) is “ the electrical conductivity which the 
metal would have, at absolute zero, if the Fermi level, ¢, came at the energy 
E”’. Ifeqn. (2) is valid, it may be shown (cf. Mott and Jones, loc. cit. p. 310, 
Wilson 1953, p. 205) that as a first approximation the thermoelectric power 
due to electron diffusion may be written : 


x ak?T dlno(E) 


3e oH (35) 
_ eT (@lnn(E) , alnv%(H) , alnz(Z) 
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where the derivatives are to be evaluated at the Fermi level, ¢. 

When the scattering is inelastic the situation becomes more complex, and. 
generally speaking the collision integral cannot be reduced rigorously to the 
form of eqn. (2). Consider however a simple inelastic collision problem in 
which electrons make purely inelastic collisions with ‘atoms’ each of which 
has a ground state and an excited state (both non-degenerate) separated 
by an energy, ¢; let us assume for the moment that the scattering matrix 
elements are strictly constant. In that case we write: 


df fH) exp (= e/h2l) a ae 
CE), sm i 1+exp (= €/kT) (1 (he(# + €) >) 


(fe (H + €)> 1 
Liiecerae "ae Tfexp(—</EP) (1—f,.(Z)), See 


- 


for transitions out of the electron state k to one of higher electron energy. 
and the reverse transition into the state k. A similar expression (replacing 
+ «by — e) must be added for the transitions between the state k and ones of 
lower electron energy. An elementary calculation now shows that if as a 
first approximation we replace ¢f,.(L + €)) by fy(# + «), then: 


d (BYE (FE 
| (F) = Sul oo ) (5a) 
where 
bostaed? exp ()+1 exp (y) +1 
pe eererne SIESMMRCICTTCETEE | 
E-¢. 


€ 
aoe eth 
The first term in 1/7’ derives from eqn. (4); the second term naturally 
derives from the second group of electron transitions. 

7, like 7 in eqn. (2) will be at most a slowly varying function of electron 
energy. We note that in this case the collision integral has been reduced to. 
the precise form of eqn. (2) (and we remark particularly that f, and J, are to 
be evaluated for the same electron energy, Z) but now the Pauli exelusion 
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principle has figured directly in producing an overall relaxation time, 7’, 
which now varies quite rapidly within the thermal ‘tail’ of the Fermi 
distribution. In fact if we rewrite eqn. (5b) as: 


1 
ey 


/ ~ 


¥ 
then 4 varies in the neighbourhood of £ when 8 ~ 1 such that: 


1) 0A 1 


Cea I were 
Siri ance ean ce 


(cf. sketch of A in the figure). Since the collision integral in this case 
has been reduced to precisely the form of eqn. (2), we can immediately 


Variation of A with y for various values of f (full lines). The two dashed curves 
show separately the behaviour for B=2 of the two terms in eqn. (5), 
which determine 1/7,’ and 1/74’. 


write eqn. (36) or (3c) for the diffusion thermoelectric power, but of 
course we must now write: 
dln7'(£) 
(“a ) 
as the last term in eqn. (3¢). 

The fact that 7’(H) varies now very rapidly with the electron energy (cf. 
eqn. (6)) might lead us to suppose that a very large thermoelectric power 
would result whenever inelastic scattering occurs, but if we examine eqn. 
(5b) more carefully we see that A, although highly energy-dependent, is 
symmetrical in y, and it appears (cf. de Vroomen and Potters 1961 ) that this 
is a rather general consequence when only inelastic scattering is involved. 


The contribution from this factor to the last term in eqn. (3c) when 
evaluated at H=¢ (y=0) will be zero, and therefore we cannot expect to 


4H2 
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encounter an anomalous thermoelectric power from inelastic scattering 
alone; this has also been pointed out by de Vroomen and Potters (loc. cit.) 
by considering the Peltier heat appropriate to the thermoelectric question. 

Recent work by Kasuya (1959), Bailyn (1961) and de Vroomenf (1959 ; 
see also de Vroomen et. al. 1960, de Vroomen and Potters 1961) on the 
anomalous transport properties of dilute alloys containing magnetic ions as 
solute has prompted us to consider a simple formal model where two groups 
of electrons contribute to electron transport. Let these two groups have 
distribution functions ¢ and %. We suppose that each group has a relaxa- 
tion time 7, and 7, independent of the other and for convenience let us 
assume these relaxation times are elastic, although this does not appear 
mandatory. We suppose in addition that the two groups are inelastically- 
coupled to one another; an electron transition from ¢ to % involves a gain 
of electron energy «, and a transition in the reverse direction a corresponding 
Joss in energy. Taking our cue from eqns. (4) and (5) we suggest as a 
crude approximation for the mutual inelastic scattering : 


Udy he — Px 2 
REE (7a) 
db, ob — dy . 


1/7,’ is similar to the first term on the right-hand side of eqn. (55), and 1/7,’ 
similar to the second term. Both 7;’ and 7,’ are highly energy-dependent 
relaxation times, but are asymmetrical about the Fermi energy (7 =0) 
(cf. the figure) ; however we note that 7,'( +1) =7,4'(—7) (ef. also de Vroomen 
and Potters, loc. cit.). Thus altogether we write: 


(St) ai Oe du — Po wx bx — Yt (Sa) 
dt coll. 77 T3. pes 
and 
dob, — bao te — be 
( dt ) coll. T3 T, ‘ °) 
Equations (8a and 6) can then be readily solved when an electric field 
(say H,) is applied to the common system and we find that: 


eH,7'. of 
= op = — = 0 
Se= Pu t+ te=ho m Ov, (9) 
where 
aale e ig ( 13 O(t 474 —TTs ) (10) 
2 T(TyTq +7973 +73 Ty ) tl iar 
where T=74+72, 
and 6= 71—T.. 


Since the overall system is characterized by a single relaxation time 7’ 
as given by eqn. (10) we can once more apply eqn. (3c) directly and expect 
to discover anomalously large thermoelectric powers when 7’ is both highly 
SG Se 
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energy-dependent in the Fermi tail and asymmetrical about n=O) We 
note that if 6=0, implying that the independent relaxation times of the 
two groups are equal, the inelastic coupling characterized by 7,/ and 7,’ will 
not contribute in any way to 7’ and no thermoelectric power anomaly will 
be expected. Also if the inelastic scattering is absent (75/, 74/—> 00) or 
overwhelmingly strong (73’, 74’>0) the model predicts again no highly 
energy-dependent term in 7’. Thus eqn. (10) shows directly within the 
limitations of our model that an anomalous thermoelectric power is only to 
be looked for in such a system if there is appreciable, but not overriding, 
inelastic coupling, and if the independent relaxation times of the two 
groups of electrons differ significantly from one another. 

Ifz,’ and 7,’ do contribute significantly then this will only be in the tem- 
perature region where inelastic collisions are frequent, i.e. when B~1, 
(when § <1, the collisions become essentially elastic and 73’ and 7,’ will cease 
to be highly energy-dependent, while when B>1 the inelastic collisions 
become very rare). But then, putting the last derivative in eqn. (3c) 
of order (but rather less than) 1/k7' (cf. eqn. (6)), we see immediately that 
a thermoelectric power somewhere approaching the order of magnitude 
(k/e) (about 86 wv/deg.) might arise. 

It is interesting to note that this is the same order of magnitude as one 
expects for the limit of ‘phonon-drag’ thermoelectric power on a simple 
argument, and it is not difficult to see that there is at least a formal analogy 
between the two processes; the present situation could perhaps be referred 
to as ‘mutual electron-drag’. 

The same conditions for the existence of an anomalous thermoelectric 
power (viz. presence of both elastic and inelastic collisions and unequal — 
elastic collision times for the two groups of carriers) have been reached by 
Kasuya (loc. cit.), Bailyn (loc. cit.) and de Vroomen (loc. cit.) for the 
system where we have conduction electrons being scattered in a metal by 
relatively dilute concentrations of magnetic impurity ions (transition 
elements). These ions are assumed to have their magnetic degeneracy 
split by some common ‘ pseudo-Weiss field’. In that case the two groups 
of electrons are those with plus and minus spins; the independent (and 
elastic) transitions are those in which an electron is elastically scattered 
by an ion without a ‘spin-flip’, while the mutual inelastic coupling arises. 
when an electron is scattered with a spin-flip, and correspondingly the- 
magnetic ion involved makes a transition from one energy level to another. 
We are certainly inclined to believe that the Kasuya—Bailyn—de Vroomen. 
model (which arises as a logical extension of an earlier model by Schmitt. 
(1956); cf. also Brailsford and Overhauser (1960)) is probably the right one 
for interpreting the highly anomalous thermoelectric powers observed. 
experimentally at very low temperatures (cf. e.g. MacDonald et al. 1958, 
1960, Gold et al. 1960). However, our elementary discussion of the 
problem leads us to ask also whether other situations could not arise in 
which two groups of conduction electrons, mutually coupled by some 
appropriate inelastic scattering (but not necessarily involving spin-flips), 
might lead to a similar result. More specifically we wonder whether the 
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parent transition metals themselves, involving s-band and d-band con- 
duction electrons coupled presumably by some inelastic transitions, might 
not give rise to a similar kind of situation. 

Note added in proof.—Kasuya (loc. cit.) has suggested that a ‘ spin- 
flipping ’ inelastic electron scattering mechanism within a single band 
may lead to an anomalously large thermoelectric power in pure transition 
metals. 

Finally we can see quite readily that this strong energy dependence 
of the overall relaxation time can also give anomalies in the electrical 
conductivity. Mott and Jones (loc. cit., p. 266) give for the actual 
conductivity, o,: 

* eye 2oal8) Pe Nee toe Ae Ts 
where o,,(£) is the conductivity we would expect if the Fermi ‘tail’ were 


quite sharp at energy ¢. 
On our present model this can be written 


8a ian Be, a) cast eee 
< 1+ ( ae | (12) 


Normally the second term in eqn. (12) is of order (k7'/¢)? and can therefore 
be neglected, but in the present situation it appears possible to get signifi- 
cant values of (1/7’)(0?7’(y)/0n?) for B~1. Furthermore this will give a 
temperature dependent term in o,, as the value of the derivative will 
depend on temperature through 8 (cf. the figure). 

Consequently eqn. (12) predicts directly on this model that a resistive 
anomaly should also be looked for arising directly from the highly energy- 
dependent relaxation time. 

Note added in proof.—However, more detailed estimates indicate that 
it is very difficult to account reasonably for the magnitude of the resistive 
anomalies typically observed. 
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Rotational Bands in the 1p-Shell 


By A. B. Curae 
Clarendon Laboratory, Oxford 


[Received March 9, 1961] 


ABSTRACT 


It is noted that certain states of lp-shell nuclei are strongly excited in 
the inelastic scattering of medium energy protons from these nuclei while other 
states are weakly excited. It is suggested that these observations can be 
understood if these nuclei are distorted, with rotational bands. The con- 
nection with other nuclear models of the 1p-shell is discussed. 


§ 1. INTRODUCTION 


In this note we propose that the collective model can be applied to the 
1p-shell in the same way as it has been applied to the 2s—1d shell (for a 
summary see Gove 1960). This suggestion is based on the observation 
that in the inelastic scattering of protons of 150-185 Mev certain states 
of lp-shell nuclei are excited strongly and others very weakly (Tyrén 
and Maris 1958, Clegg et al. 1961), noting that it has been shown by 
Clegg et al. that such inelastic scattering provides an almost quantitative 
means of picking out strong electric quadrupole radiative transitions to 
the ground state. Thus, in odd-mass nuclei in the 1p-shell, we identify 
low-lying states which are strongly excited by inelastic scattering as 
states of the same rotational band as the ground state while other, weakly 
excited, states are taken to be members of other rotational bands. In 
this way a good fit is obtained to the low-lying levels of these nuclei, the 
only intrinsic states being those expected from the filling of Nilsson’s 
energy levels in deformed nuclei (Nilsson 1955), taking prolate distortion 
in the first half of the 1p-shell and oblate distortion in the second half. 
Moments of inertia, intrinsic spacings and decoupling parameters involved 
in K =} bands are determined by empirical fitting; however, it is always. 
found that the decoupling parameter found empirically agrees very well - 
with the value calculated for the corresponding Nilsson wave-functions. 
In °Be and 1!B there are two bands, which one would expect to mix due 
to the rotation-particle coupling (Kerman 1956), but the apparent 
selection rule on inelastic scattering suggests that this mixing is not 
large. . 

The relationships between the model proposed here and other models, 
particularly the intermediate coupling scheme, which is so successful in. 
correlating so many properties of 1p-shell nuclei, are discussed in later 


sections. 
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§ 2. Exprcrep INTRINSIC STATES 


The Nilsson energy level scheme for the Ip-shell consists of three 
levels labelled level 2: pg/., K = 3; level 3: p3/2, K =4; level 4: pij,, K=3- 
The energies of each of these levels is a function of 7, the ratio of the 
spheroidal distortion to the spin-orbit coupling. 7 is positive for prolate 
distortion and negative for oblate distortion. 

Four particles, neutron or proton with Q= + K, can be placed on each 
level. Thus for mass 7, with 7 positive, three particles go in level 3, 
forming a low-lying K =} intrinsic state with a large gap before the 
next intrinsic state is reached. For mass 9 level 3 will be filled with 
four particles and the fifth will be on level 2 or 4; thus there will be two 
close intrinsic states with K=4, 3 and then a big gap before any more 
show up. For mass 13, with 7 negative, three holes are placed on level 4 
forming a low-lying K=4 intrinsic state with a large gap before the 
next intrinsic state. For mass 11 level 4 is filled with four holes and the 
fifth hole is on level 2 or 3 forming two low-lying intrinsic states of 
K =}, 3 with a large gap to the next intrinsic state. One notes that the 
intrinsic states proposed for 11B are formed in a different manner from 
those used in the fits made to the (4n+3) nuclei in the 2s—ld shell: ?*Na, 
"AJ, 31P (summarized by Gove 1960). Here we have chosen intrinsic 
states formed by a hole on otherwise filled levels, while in the fits to the 
level schemes of 28Na, 27Al, 34P a single particle has been promoted 
to unfilled levels. For example, if we had used this latter procedure in 
11B, we would have intrinsic states formed by single particles respectively 
on Nilsson levels 3 and 4, giving two K=4 bands, which do not fit the 
11B energy level scheme so well. One can argue that the intrinsic states 
- we have chosen would come lower in energy due to the pairing energy; 
on lifting a particle from a filled or partially-filled level to an empty level 
one is breaking a pair so that the corresponding state should come at a 
higher energy. That this is so for 27Al has been demonstrated by Brink 
and Kerman (1959). Further, experimental measurements of (p, 2p) and 
(p, pn) reactions on *4Mg, *8Si and °°S suggest that the low-lying intrinsic 
states in the (4n+3) nuclei are formed by holes in otherwise filled levels 
(Clegg and Foley, to be published). Thus we feel that the correct assign- 
ment of intrinsic states has been made in 'B, on other grounds than 
that they are the ones we need to fit the experimental energy spectrum. 


§ 3. Mass 11 


Tyrén and Maris (1958) see strong excitation of states of 1B at 4-7 + 0-2 
and. 6-6 + 0-2 Mev in the |B(p, p’) reaction. It has been shown that the 
lower of these two is the J =§ state at 4:46 Mev, and that the J =} state 
at 2-13 Mev is weakly excited in this reaction (Clegg et al. 1961). The 
J=}-— state at 5-03 mev is weakly excited (Clegg ef al. 1961) while the 
6-6 + 0-2 Mev state could well be the J = 3 state at 6-81 Mev so one suggests 
that the states excited are a K = rotational band based on the J= 3 
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ground state and the states not excited are a K =} band based on the 

=% state at 2-13 mev; these intrinsic states have the merit of being 
those predicted by Nilsson’s scheme. These two bands will mix due to 
the rotation—particle coupling (Kerman 1956), changing the level spacings. 
A fit has been made to the known level scheme and is shown in table 1 
where the energies quoted for the experimental levels are the means of 
the energies of corresponding levels in "B and 110. The decoupling 


Table 1 


Experimental Calculated 
State energy energy 
(Mev) 


0 
2-06 
4-35 
4-90 


3 
2 
Al 
2 
5. 
2 
3 
2 
ve 
2 
5 
2 
vi 
2 


parameter for the K=4 band is taken to be a= —0-8 which corresponds: 
to Nilsson’s level 3, 7=—6. Other parameters determined empirically 
are f?/2.% =0-85 Mev (we have throughout taken the same moment of 
inertia for both bands in the same nucleus, simply to reduce the number 
of parameters), 4,=0-60 mev (the matrix element of rotation—particle 
coupling) and the intrinsic spacing=1-71 Mev. In this fit the only 
requirement limiting the band-mixing is that it should be weak enough 
so that the individual states are still largely K = 3 or K=4 and then the 
simple features of the inelastic scattering are reproduced. An experiment 
is in progress to measure the ‘ forbidden’ transitions and so determine 
the band-mixing more accurately. A further contribution to the level 
energies is due to the rotation—vibration interaction which contributes a 
term to the energy of the form —B[J(J+1)+(—1)7t¥? b;;,.4(J + 3)]?. 
This term will obviously help to produce a better fit in the case of 1B; 
however, an actual numerical fit is postponed until the band-mixing is 
better determined. For a pure K = 3 band the ratio of inelastic scattering 
cross sections to the 6-8 and 4-46 mev states should be 0-56, which is in 
reasonable agreement with the experimental results. 


§ 4. Mass 9 


In the °Be(p, p’) reaction states strongly excited are at 2340-2 Mev 
and ~6mev. The first is presumably the J= % state at 2-43 mev and 
the other is most probably the state at 6-76 Mev. Again one supposes 
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that these statés are a rotational band based on the J = 3 ground state and 
it is immediately found that a good fit is obtained with h?/2.% = 0-524 Mev: 
see table 2. 

However, there are other well-established states in *Be which are 
presumably due to the expected K=} band. This band will mix with 
the K = 3 band, but there is not enone information about the spins and 
parities Sf the other states of *Be to justify a detailed fit. That these 
states in °Be form a rotational band has previously been suggested by 
Blair and Henley (1958) from the strong excitation of them in the inelastic 
scattering of 40 Mev alpha-particles. As in the 1B(p,p’) reaction the 
ratio of cross sections for exciting the two states in *Be(p, p’) should be 
0-56 for a pure K=3 band, a result in reasonable agreement with 
experiment. 


Table 2 


~ State 


S Model 
Experiment 


1b 
o> 00 


Energies 


oo 
i o> 
Ww bo 
enor) 


§ 5. Mass 7 


In the “Li(p ,p’) reaction the strongly excited state is at 4-4+0-2 Mev, 
which presumably corresponds to the state at 4-63 Mev, whose mirror 
state in 7Be has J=3—. It is easily found that the states of 7Li can be 
fitted with a K =4 band with a= — 1-3, calculated for level 3 with n= +8, 
and ?/2. =0-54 Mev. This gives an energy level scheme: 


State j=3- J=4— Je J=2 J=2 
Model 0 0-49 5:08 6-21 16-24 
Experiment 0 0-46 4-58 7-32 ~- 


The experimental energies quoted are the mean of those for 7Li and 
7Be. For such a K=} band the inelastic scattering cross sections will 
be proportional to: 


BEB —-L,) =[(3220 | 174)Qo + (32-31 | 1,3)Q, 7? 
-where Q), Q, depend on the intrinsic wave-functions. So we have 
I; B($—1,) 
[—(4)/?Qo + (35)7/7Q11? 
[(88)¥/2Q9 + (8) "20,7 
[(g5)/@o — (i) 7011? 


We see that if Qo, Q, are of comparable magnitude the transition to 
the J = § state at 4-63 mev in the 7Li(p, p’) reaction will be strong and the 
feenaition to the J= § state at 7-47 Mev will be weak, as is observed. 


boler pea 
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This predicts that the transition to the J=4 state at 0-48 mev will also 
be weak; measurements to check this prediction are in progress. 


§ 6. Masszs 8, 12, 13 
It has long been realized that *Be displays a rotational spectrum, 
though this has usually been described as an aspect of the alpha-particle 
model. 
Apparently h?/2.% =0-53 gives a good fit: 


State J=0+ J=2+ J=4+ 
Model 0 3°20 10-8 
Experiment 0 2-90 11-7 


Taking the J=2+ state of 12C at 4:43 Mev to be a rotational state we 
find h?/2. = 0-74 mev while similarly the J = 3— state of 13C at 3-68 Mev 
and of 15N at 3-51 Mev, corresponds to #?/2. = 0-96 ev, taking a=0-25 
calculated for level 4 and 7=—3. 


§ 7. CONCLUSIONS 
First we summarize the values of #?/2.% found: 


A 7 8 9 i 12 13 
h?/2F 0-54 0-53 0-52 0-85 0-74 0-96 


These agree very well among themselves, falling into two groups: the 
prolate nuclei in the first half of the shell and oblate nuclei in the second 
half of the shell. Also the values of the decoupling parameter that are 
required for an empirical fit agree very well with those calculated for 
Nilsson wave-functions for reasonable values of 7» and the intrinsic levels 
found to be needed agree with those expected from Nilsson’s energy levels 
so that one can imagine that this model has some physical significance. 
It presumably is a first approximation to the work of Kurath and Pi¢éman 
(1959) who show that wave-functions constructed by the generating 
procedure of Peierls and Yoccoz (1957) have a very considerable overlap 
with wave-functions from an intermediate coupling calculation which 
fits the energy level spectra of Ip-shell nuclei (Kurath 1956). It has been 
shown by Yoccoz (1957) that this generating procedure produces an 
energy level spectrum approximately the same as a rotational band if 
the nucleus is sufficiently distorted. It is perhaps therefore significant 
that the more distorted nuclei 7Li, *Be, °Be show closer approximations 
to rotational bands while in the case of the less distorted ‘1B there is 
need of the rotation—vibration interaction term to obtain a more exact 
fit to the energy level scheme. Thus it can be suggested that the work 
described here shows how medium energy inelastic scattering results 
support the picture of Kurath and Piéman and helps to identify the 
physical states corresponding to their wave-functions. Further, this 
work suggests that, in Be and ™B, there is very little mixing between 


the rotational bands. 
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§ 8. (p,2p) ReactTIONS IN THE 1p-SHELL 


This picture can apparently also explain the results of (p, 2p) reactions 
at 150-185 Mev in the lIp-shell (Tyrén ef al. 1958, Gooding and Pugh 
1960) where it is found that only a few states are excited strongly in the 
residual nuclei. These results are taken to imply some simple picture 
of nuclear structure, where the ground state of nucleus A has very few 
parent states in nucleus A-1, and have been interpreted by a simple 
jj-coupling picture. This interpretation would seem to conflict with 
the intermediate coupling picture which has been found necessary to 
explain most low energy data, for example, the energy level spectra of 
1p-shell nuclei (Kurath 1956), as the intermediate coupling wave-functions 
are apparently complicated so that the ground state of nucleus A would 
presumably have many parent states in nucleus A-1, contrary to the 
results of the (p,2p) reactions. However, the picture presented here 
would describe the ground state of nucleus A as an intrinsic state made 
up of Nilsson wave-functions. Hjecting a proton would leave only a 
few intrinsic states of nucleus A-1l, so that. this picture presumably 
shows the simplicity of the intermediate coupling wave-functions necessary 
to understand the (p, 2p) reactions. 


§ 9. COMPARISON WITH OTHER MODELS 


It has already been mentioned that the rotational band in Be has. 
been previously ascribed to the alpha-particle model, a model in which 
the *Be nucleus is described as made up of two alpha-particles in a 
dumb-bell arrangement. It has however been shown, by Perring and 
Skyrme (1956), that such a model is equivalent to the shell-model picture, 
if the nucleon wave-functions are fully antisymmetrized. Thus if our 
model is an aspect of the shell-model there should be no problem in 
reconciling it with the alpha-particle model. (One notes that the alpha- 
particle model for the §Be ground state is a prolate nucleus, and the 
alpha-particle model for the °C ground state—three alpha-particles at 
the corners of an equilateral triangle—is an oblate nucleus, in agreement 
with our model.) 

The alpha-particle model of SBe has been extended to 7Li and °Be, 
by adding respectively one hole and one particle to 8Be (Kunz 1960). 
In this way Kunz obtains distorted nuclei with rotational bands much 
_the same as those described in this paper. This picture is presumably 
the same as adding one hole or one particle to the shell-model wave- 
functions of *Be so we see that our model and that of Kunz are presumably 
just two equivalent ways of looking at the same wave-functions. 

One can see that the alpha-particle picture of 7Li is equivalent to the 
cluster-model picture due to Phillips and Tombrello (1960), where the 
“Li states are taken to be made up of an alpha-particle cluster and a. 
triton cluster, providing one adds to the cluster model the requirement. 
of full antisymmetrization of the nucleon wave-functions. So we can 
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take the evidence for the cluster-model as evidence for our model. A 
typical further example of agreement is the case where the 2-43 Mev 
state of *Be is shown to be made up of largely the 2+ state of *Be plus 
a neutron, in agreement with our assignment of both these states as the 
first excited states of rotational bands. 
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ABSTRACT 


From measurements of the radius of curvature of extended dislocation nodes 
observed by transmission electron microscopy, stacking-fault energies have 
been determined for a number of copper and silver base aluminium and zine 
solid solution alloys and also for a number of nickel-cobalt alloys. Rough 
estimates of about 40, 25 and 150 erg em~? respectively can be made of the 
stacking-fault energies of copper, silver and nickel. The accuracy of the 
measurement and the possibility of segregation of solute atoms to stacking 
faults are discussed, 


§ 1. INTRODUCTION 


In view of the importance of the stacking-fault energy of the face-centred 
cubic metals in determining the equilibrium ribbon width of extended 
dislocations (Heidenreich and Shockley 1948) and hence in theories of 
work hardening (Seeger 1956, Hirsch 1959, Basinski 1959) and its further 
potential importance in the theory of alloy phases because of its relation 
to the energy difference between the hexagonal close-packed and. face- 
centred cubic crystal structures, it is unfortunate that there are so few 
metals for which the value of the stacking-fault energy is accurately 
known. Two of the earliest estimates of stacking-fault energy remain 
the most reliable: that of » +40 erg cm~* for copper deduced by Fullman 
(1951) from his determination of the twin boundary energy as ~ 20 erg 
em-2 (on the assumption that the stacking-fault energy is twice the 
twin-boundary energy) and that of an internal energy U = — 20 erg cm~* 
(corresponding to a free energy y ¥ —7 erg cm~*) estimated for cobalt by 
Heidenreich and Shockley from the known heat of transformation in 
the f.c.c.>h.c.p. transition. The values for cobalt are taken to be 
negative since it is convenient to regard y as positive when the f.c.c. 
structure is stable. A value of yx 200 erg cm~? for aluminium was also 
estimated by Fullman. 
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Estimates of stacking-fault energy have also been made for a number 
of metals from observations of the mechanical properties interpreted in 
terms of a theory of work hardening (Thornton and Hirsch 1958, Seeger 
et al. 1959). Unfortunately, however, these results all depend on the 
validity of theories still held in question and the stacking-fault energy 
tends to appear simply as yet another adjustable parameter in the 
theory rather than as an independently known quantity. 

The effect of stacking faults on x-ray diffraction has been used as a 
means of determining stacking-fault densities but since the theory applies 
only to faults of large area and the effect of stacking-fault ribbons is 
unknown, the interpretation of the results in terms of stacking-fault 
energy (Christian and Spreadborough 1957) has only a qualitative 
significance. 


Fig. 1 


Diagram showing the arrangement of partial dislocations (for notation see 
Whelan 1959) to form a network of extended and contracted nodes. 
The shaded area is the stacking fault. 


The application of transmission electron microscopy to the study of 
dislocations in thin foils showed that extended dislocations and stacking 
faults can be observed in metals of low stacking-fault energy (Whelan 
et al. 1957). Many interactions between extended dislocations have been 
identified and in particular one which gives rise to the formation of 
hexagonal networks lying on the slip planes with alternate nodes in the 
network being of an extended or contracted type (see fig. lt) has been 


\ } Figures 5 and 6 are shown as plates. 
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used by Whelan (1959) to obtain an estimate of the stacking-fault energy 
of stainless steel. Whelan considered the partial dislocation at A to be 
in equilibrium under a force @b,2/2R per unit length, due to its line 
tension (the line tension is taken as Gb,?/2; b, is the Burgers vector of 
the partial dislocation and R its radius of curvature) and a force y per 


unit length due to the attraction of the stacking fault. From the simple 
formula 


y=Gb7/2R aa Te: a Paneer EL 
the stacking-fault energy of stainless steel was estimated to be ~13 


ergcm™*. The authors have determined stacking-fault energies in a 
number of alloys by this method. 


§ 2. ACCURACY OF STACKING-FAULT ENERGIES DEDUCED FROM 
MEASUREMENTS OF DisLocaTION NopEs 


Before describing the results obtained, some discussion of the limita- 
tions of eqn. (1) should be given. The line tension 7 of a dislocation 
depends on its character, i.e. on its edge and screw components and also 
on its configuration. It is only approximately described by the formula 
Gb,?/2 and a more accurate formula (see Friedel 1953) gives 

T=Gb? In (R/b)/4a0K 
for a dislocation of radius of curvature R, where K has a value of 1 for 
a pure screw dislocation and (1—v) for a pure edge dislocation (v being 
Poisson’s ratio). 6) has approximately the magnitude of the Burgers vector 
of the dislocation. Thus an improved version of eqn. (1) is 

y =Gb? In (R/b,)/40K BR. SR ee ee, 
The most serious error in eqn. (2) is probably that it neglects the stresses 
on the partial dislocation «B at A (fig. 1) due to the other partial disloca- 
tions «C and D«. At present there is no precise calculation of the effect 
of these stresses although computations are in progress (Bilby eé al. 
1960). It is possible however to make some estimates of the error 
involved in applying eqn. (2). The stress at A due to the two dislocations 
aC and De. will certainly be considerably less than that due to a long 
straight partial dislocation at a distance AB (fig. 1) from A. Since in 
practice AB x R/2 this indicates that the force of repulsion on «B due 
to «C and Doe is less than Gb,?/27R. The values of stacking-fault energy ~ 
derived from formula (2) might therefore have to be increased by a factor 
of up to 1-3. An alternative estimate of the error can be obtained by 
replacing the two dislocations aC and Dz by dislocation loops of radius R. 
The stress at A from such loops has been calculated using the formulae 
given by Kréner (1958) with some corrections pointed out by Kroupa 
(private communication) and indicates that the values of stacking-fault 
energy obtained from (2) should be multiplied by about 1-1. This second 
result probably underestimates the effect but is thought to be the more 
accurate of the two correction factors. 

The effect of external stresses on the node will also give rise to an 
error in eqn. (2) and can only be minimized by confining the measurements 
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to reasonably symmetrical nodes not too close to other dislocations. 
(Some asymmetry will usually be present because the nodes are almost 
always observed in projection.) 

Further errors in the determination of R could result from contrast 
effects since the image of a dislocation does not necessarily coincide with 
the dislocation (Hirsch et al. 1960). However, nodes are usually observed 
in fairly thick foils near the Bragg diffracting position for a low-order 
reflection and under these conditions the dislocation and its image are 
nearly coincident (Howie and Whelan, to be published) and the error in 
measuring R is not more than ~ 204, i.e. about 10% for the smallest 
nodes observed. 

§ 3. EXPERIMENTAL RESULTS 

Stacking-fault energies have been measured for several compositions 
of the primary solid solutions of copper and silver base alloys containing 
aluminium or zinc and also for the f.c.c. nickel-cobalt alloys near the 
transformation composition. The nickel—cobalt alloys were donated in 
sheet form by the Mond Nickel Co. Ltd. The other alloys which were 
made in the laboratory were homogenized and rolled into sheet 250 uw 
thick between several annealing treatments in an argon atmosphere. 
All the alloys were analysed by Johnson Matthey and Co. Ltd. after the 
final annealing treatment. The annealed sheet was then lightly deformed 
and electro-thinned for examination in the electron microscope. The 
electro-polishing solutions and conditions used in the final thinning process 


are shown in table 1. 
Table 1 


Electrolyte Voltage Temperature 


35°, HNO; in methyl! alcohol 9-15 —20°c 


6 g KCN in 100 cm? water 6-8 


70% ethyl alcohol 20-40 
20°, perchloric acid 
10% glycerol 


80% ethyl aleohol 
20°, perchloric acid 


Electron micrographs of several symmetrical isolated dislocation nodes 
were then obtained for each alloy using the Siemens Elmiskop I electron 
microscope operating at 100 kv. Examples are given in figs. 5 and 6. 
The projected radii of curvature R’ of the partial dislocations at each 
node were then measured and whenever possible corrected for inclination 
to give the true radius of curvature R from the formula 


R= R’ sec* 6 (cos? « + cos? 6 sin? «)3/2 
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where @ is the angle between the electron beam and the normal to the 
plane containing the node and « is the observed angle between the line 
of intersection of this plane with the foil surface and the normal to the 
partial dislocation at the point of measurement. The angles 6 and « 
were calculated from selected area diffraction measurements. For metals 
in which the nodes are small (R < 400 A) the dislocation arrangement is 
irregular and it is often impossible to be certain of the slip plane on which 
the nodes lie and hence to correct for the inclination. However, the 
correction factor is sometimes greater and sometimes less than unity 
and it does not seem probable that there is any large systematic error 
caused by ignoring it in these cases. 


Fig. 2 


Cu+Zn alloys 


: 


Cut Al alloys 


oO 
LOO HO 120 130 


Electron-atom ratio 


Stacking-fault energy in units of 10-3 Gb for Cu+Al and Cu + Zn alloys plotted . 
against electron-atom ratio. For pure Cu, 10-* Gb~ 10-5 erg em~*. 


The stacking-fault energy of each alloy has been calculated on the 
basis of at least ten different measurements. Since the values of the 
shear modulus G vary somewhat with alloy composition, the stacking 
fault energies are given in units of 10-* Gb where b is the Burgers vector 
of a whole dislocation. Using the values of shear modulus quoted by 
Seeger et al. (1959), the quantity 10~° Gb is found to be equal to 10-5, 
7-1 and 18-6 erg cem~ for pure copper, silver and nickel respectively. 
The results for the five alloys are given in tables 2-6 and are plotted 


A412 
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9.00 HO 120 30 


Electron-atom ratio 


Stacking-fault energy in units of 10-* Gb for Ag+ Al and Ag +Zn alloys plotted 
against electron-atom ratio. For pure Ag, 10-3 Gbx7-1 erg cm-?. 


Fig. 4 


SO 60 70 80 
Weight % Cobalt 


Stacking-fault energy in units of 10-*Gb for Ni+Co alloys plotted against 
weight % Co. For pure Ni, 10-? Gb= 18-6 erg em-?. 
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against electron-atom ratio (or composition in the case of the nickel— 
cobalt alloy) in figs. 2, 3 and 4. Since no preferred orientation of the 
nodes in the slip plane was observed, an average value of 5/6 for K 
(eqn. (2)) was employed throughout. 


Table 2. Cu+Zn alloys (maximum total impurity content 0-08°%) 


Wt. % Zn 


Electron-atom 1-10 1-14 1:20 1-30 
ratio 


103 y/Gb 1-68+ 0-15 1-21+0-10 0-86+ 0-09 0-67+0-05 | 0:53+0-06 


Table 3. Cu+Al alloys (maximum total impurity content 0-1%) 


Wt. %-Al 2°25 3°53 


Electron-atom 
ratio 1-11 1-16 1-20 1-31 1-34 


108 y/Gb 1-77+0-13 0-934 0-13 0-35+ 0-07 0-23+ 0-04 0-17+ 0:03 | 


Table 4. Ag+Zn alloys (maximum total impurity content 0-:08%) 


Wt. % Zn 1:94 3°30 5°00 15-00 20:00 
Electron-atom 

ratio 1-032 1-054 1-08 1:23 1-29 
103 y/Gb 1-63+ 0-12 1-11+0:19 0-914 0-07 0:53+ 0-07 0-36+ 0-06 | 


Table 5. Ag+ Al alloys (maximum total impurity content 0-08%,) 


Wt. % Al : 1-00 2. 5-00 


Electron-atom : 
ratio 1-078 ° 1-35 


103 y/G@b 55+ 0-5 1-1240-17 0-58+ 0-07 0-1440-07 


Table 6. Ni+Co alloys (maximum total impurity content 0-10%) 


Wt. % Co 54-50 59-40 64-29 69-48 


108 y/Gb 1-59 + 0-80 1:09+ 0-05 0-64+ 0-04 0-29+ 0-03 
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§ 4, Discussion 


4.1. Values of Stacking-fault Energy 


There is as yet no complete calculation of the stacking-fault energy 
of any metal or any theory of its dependence on composition for any 
alloy system. It is possible however, and convenient for the subsequent 
discussion, to list the various contributions to the stacking-fault energy. 


4.1.1. Energy of nearest-neighbour interactions at the fault 


As a result of the different arrangement of neighbours round the atoms 
on the fault compared with the other atoms in the crystal there will be 
a change in energy of the core electrons leading to a change in the total 
energy of the crystal. 


4.1.2. Energy due to scattering of conduction electrons 


The fact that stacking faults strongly reflect conduction electrons has 
recently been convincingly demonstrated (Cotterill, to be published) by 
measurements of the increase in resistivity due to the presence of tetra- 
hedra of stacking faults observed in quenched gold (Sileox and Hirsch 
1959). Attree and Plaskett (1956) have calculated the stacking-fault 
energy due to this scattering in a monovalent metal with almost-free 
conduction electrons and find that part of the energy comes from the 
electrons which are reflected by the stacking fault and part from electrons 
which are localized in surface states (Tamm waves) on either side of the 
fault. 

The results of figs. 2 and 3 show that the addition of zine or aluminium 
to either copper or silver lowers the stacking-fault energy and that for 
low solute-content alloys the stacking-fault energy depends mainly on the 
number of valence electrons in the alloy. This implies that it is the 
energy contribution due to the scattering of the conduction electrons 
which decreases with increasing number-of valence electrons. It is also 
striking that the slope of the curves for the copper and silver alloys should 
be similar in this region and this suggests, in agreement with recent 
experimental results for the pure metals (Shoenberg 1960), that the 
conduction electrons in these copper and silver alloys behave similarly. 

The similar slopes of the curves of figs. 2 and 3 in the region of low 
solute-content enable us to estimate with more confidence than would 
otherwise be the case that the stacking-fault energies of copper and silver 
are in the vicinity of 40 erg em and 25 erg cm respectively. This 
value for copper agrees well with Fullman’s estimate quoted earlier. 
From fig. 4 we estimate that the stacking-fault energy of nickel is about 
150ergcem™. Although these stacking-fault energies can only be 
regarded as rough estimates it does not seem possible to reconcile them 
with the figures of 169 erg em~*, 43 erg cm and 410 erg cm quoted 
for copper, silver and nickel respectively by Seeger et al. (1959). Even 
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when full allowance is made in our results for systematic errors of the 
type discussed in § 2 it does not seem probable that the figure for copper 
can be much higher than 60 erg cm-*. The curves of fig. 2 would have 
to rise very steeply indeed for low solute-content alloys. 

It does also seem probable that the method of estimating stacking-fault 
energies by observing whether loops or tetrahedra are formed on quench- 
ing is in error when based on simple energy considerations for defects of a 
given size. This method led to the high figure noted above for copper 
and also to an estimate of stacking-fault energy for Ni+60°% Co alloy 
of between 5 and 7 x 10-3 Gb (Mader and Simsch 1961). Both of these 
figures exceed our results by a factor of about 5 and it therefore appears 
that for some reason, probably connected with the rate of formation or 
shape of the dislocation loop, that the change to the tetrahedral configura- 
tion of stacking faults does not occur even though it is energetically 
favourable. 

In the case of the nickel—cobalt alloys where an f.c.c.—>h.c.p. transition 
occurs (at about 75°, Co at room temperature) the stacking-fault energy 
is directly related to the phase change and must pass through zero at 
this composition, becoming negative for alloys of higher cobalt content. 
The curve of fig. 4 would indeed appear to extrapolate to zero at this 
composition. No such simple considerations apply to the copper and 
silver alloys since beyond the solid solution limit (at an electron-atom 
ratio between 1-36 and 1-41 in all cases) complex phases are formed. It 
is possible of course that stacking faults act as nuclei for these new 
phases and that the stacking-fault energy plays some role in the transition 
but at present we have no evidence of this. In the case of all the copper 
and silver alloys, however, hexagonal phases are eventually formed at 
still higher solute contents so that the stacking-fault energy can still 
be related to the relative stability of the f.c.c. and h.c.p. structures. 


4.2. Segregation Effects 


The dependence of stacking-fault energy on composition implies a 
chemical interaction between the solute atoms and the stacking fault so 
that the question of segregation of solute atoms to the fault must -be 
considered. Segregation of this type was first proposed by Suzuki (1952) ~ 
to explain the initial yield points of various alloys and has been analysed 
further in this connection by Suzuki (1956) and Flinn (1958). It has also 
been suggested as a possible explanation of various other yield and 
strain-rate phenomena observed during the deformation of alloy crystals 
(see for example Bolling 1959). Such a hypothesis however implies that 
segregation can take place at temperatures well below those at which 
solute diffusion by the usual process of the creation and migration of 
vacancies is operative. Thus, although there is no doubt that segregation 
will tend to take place in principle, the interesting questions of how and 
at what temperature it occurs are still unanswered. 
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All the measurements of stacking-fault energy reported here refer to 
specimens deformed at room temperature with no subsequent heat 
treatment. The authors have carried out several experiments to investi- 
gate the effects of such treatment on the size of the nodes. Ageing at 
250°c for periods of 1 hour (after deformation at room temperature) 
produced no effect in any of the alloys studied (Cu+7% Al, Cu+4:5% Al, 
Cu+30% Zn, Ni+55% Co, Ni+60% Co). In the case of the copper- 
aluminium alloys this result is in strong disagreement with the recent 
suggestion by Cahn and Davies (1960) (on the basis of x-ray small angle 
scattering results) that segregation takes place in these alloys at a 
temperature of about 200°c. Ageing experiments with the same five 
alloys were also carried out for times of between a few minutes and 1 hour 
at temperatures of 450°c and only in the case of the Cu+ 7% Al alloy 
was any effect observed. With this alloy the nodes were observed to be 
smaller, indicating that the stacking-fault energy had risen during the 
ageing treatment from about 1-5erg cm to 5ergem~*. This change 
was observed in specimens furnace cooled from 450°c and in specimens 
of the alloy quenched from 450°c. A possible interpretation is that 
solute atoms (already segregated at the fault) are leaving the fault for 
entropy reasons. However, there is no explanation for the absence of 
the effect in the other alloys and the observation that the specimens, for 
which the node size changed, showed a tendency to form pits during 
electropolishing was disquieting. 

The apparent non-occurrence of segregation in the other copper alloys 
(Cu+4:5% Al, Cu+30% Zn) on ageing at 450°c where solute 
diffusion by creation and migration of vacancies is easy suggests that 
either segregation has already taken place at room temperature or else 
that the force of attraction between the solute atom and the fault is 
weak. Possible alternative processes of solute diffusion are: 


(a) Diffusion involving the migration of single vacancies created during 
the deformation. 

(6) Diffusion involving the migration of divacancies or other vacancy 
complexes created during the deformation. 


(c) Diffusion of interstitial solute atoms created during the deformation. 

(d) Pipe diffusion along dislocations. (It seems unlikely that this 
process would be an effective way of transferring solute atoms from 
the matrix to the vicinity of the stacking faults.) 


In copper where £,,, the activation energy for vacancy migration, is 
quoted as 0-7+0-lev (Kimura et al. 1959) segregation by process (a) 
could readily occur even at 0°c since a vacancy would make several 
jumps per second. It seems fairly improbable however that this process 
could be responsible for segregation at 0°c in nickel since even if HE. is 
as low as 1-0 ev (Sosin and Brinkman 1959) several hours would be 
needed for each jump and many of the specimens were examined in the 
microscope only 30min after deformation. It therefore appears that 
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process (a) can almost certainly be eliminated as a segregation mechanism 
in nickel since it would not be effective at the temperature of deformation 
and no segregation effects were observed even at higher temperatures. 
Segregation by process (a) may be ruled out in the case of copper and 
silver as well if we make the simplifying assumption that segregation 
occurs by the same mechanism in all the alloys. 

We may conclude therefore, that although the possibility of segregation 
of solute atoms to stacking faults at room temperature or below still 
exists, it is unlikely to occur by any process involving the migration of 
single vacancies. If segregation does not occur at these temperatures we 
can conclude that solute atoms are only weakly attracted to stacking 
faults since it does not occur at higher temperatures either. A definite 
conclusion about the possibility of room temperature segregation could 
perhaps be reached if stacking-fault energies were measured after 
deformation at low temperatures without ageing at room temperature. 


§ 5. CONCLUSIONS 


By measuring the radius of curvature of extended dislocation nodes 
formed during deformation the stacking-fault energies have been deter- 
mined for the solid solutions of aluminium and zine in copper and silver 
and also for nickel-cobalt alloys over a range of compositions. The 
results suggest that the stacking-fault energies of copper, silver and 
nickel are about 40, 25 and 150 erg cm~ respectively. The method is a 
powerful and direct way of measuring stacking-fault energies (provided 
they are not too high) and it should be possible to use it to obtain informa- 
tion about many other alloy systems. In practice, nodes with radii of 
curvature in the range 200 4—4000 A (or much more if the node lies in the 
plane of the foil) can be measured. 

It is not known at present whether segregation of solute atoms to the 
stacking faults takes place or not, but the experiments suggest that if it 
does not occur at room temperature, then it does not occur at all, even at 
temperatures where diffusion of solute atoms can take place. Further 
experiments (particularly at low temperatures) are needed to settle this 
question. 
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ABSTRACT 


Single crystals of aluminium were grown by the Czochralski method at 
various rates. The crystals were carrot-shaped, about 4 cm long, tapering 
from a diameter of about 3mm at the top to about 0-1 mm at the bottom. 
X-ray diffraction micrographs were taken of these specimens by means of 
Lang’s technique, which is capable of revealing the dislocation structure 
in the interior of a crystal. It was found that the density of dislocations 
present in an as-grown crystal increases with both increasing diameter and rate 
of growth. For sufficiently small diameters (below about 0-3mm) it was 
possible to grow crystals free of detectable dislocations. These experimental 
results are consistent with the idea that the majority of the initial dislocations 
are produced by the collapse of vacancy discs that are nucleated when a 
sufficient supersaturation of vacancies is reached in the crystal. Theoretical 
calculations of the supersaturation of vacancies that arise during crystal 
growth were made by considering the diffusion of these imperfections in a 
specimen cooled from the melting temperature. It is shown that the range 
of diameters below which the vacancy supersaturation is too small to 
nucleate a dislocation is of the same order of magnitude as the diameters of 
‘erystals found that contain no detectable dislocations. 


§ 1. INTRODUCTION 


Distocations do not exist in thermal equilibrium with respect to any 
‘conceivable sources and sinks in a crystal, particularly with respect to 
external surfaces. It is therefore necessary to think of their occurrence 
in terms of kinetic mechanisms that would operate during growth. Five 
reasonable mechanisms have been proposed (Elbaum 1959): 


(1) Propagation of dislocations into the growing crystal from the seed. 
crystal. 

(2) Shear strains arising from thermal strains or external mechanical 
constraints. 

(3) Shear strains resulting from abrupt change in composition accom- 
panied by changes in lattice parameter. 

(4) Dendritic growth. 

(5) The mechanism of collapsing vacancy discs. 


a is 
+ This work was sponsored by the Office of Naval Research under Contract 
No. Nonrr 562(27). 
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By a judicious choice of growth technique and the use of high purity 
material, the effects of some of these mechanisms may be minimized. 
A study will be described in this paper of the mechanism of collapsing 
vacancy discs, as it is thought to be fundamentally the most important 
in creating dislocations in metal crystals grown from the melt. 


§ 2. Cottapsine Vacancy Discs 


In a crystal at its melting point, a certain number or vacancies exist in 
equilibrium with external surfaces and any other sources or sinks (such as 
dislocations) that may be present. Let us consider a cylindrical aluminium 
crystal whose only source and sink for vacancies is the external surface ; 
no dislocations are assumed to be present. The concentration of vacancies 
in equilibrium with the surface, at the melting temperature, is given by 


Cy exp ("Gf RT 5) sue ieee ee 


where U; is the energy of formation of a vacancy (17500 cal/mole 
(Bradshaw and Pearson 1957)) and 7’) (= 933°xK) is the melting temperature. 
If the temperature is suddenly lowered to 7, <7), eqn. (1) states that a 
concentration of vacancies c, should now be in equilibrium with the surface. 
Since c, <c, it is apparent that at the instant the temperature is lowered, 
the concentration of vacancies in the crystal is greater than the new 
equilibrium value. The surface is now a sink for the excess vacancies 
which will begin migrating towards it. In a time ¢, vacancies will travel 
an approximate distance 

ema] (Dt) en We ee S| 


where D is the diffusion coefficient for vacancies. Thus, the concentration 
of vacancies in the crystal will begin to approach the equilibrium value. 
It is apparent that if the crystal is large, or if D is very small 
(D= Dy exp (— Um/RT) where Um is the activation energy for migration 
of vacancies) as will be the case if 7’, is rather low, then all of the excess. 
vacancies in the centre of the crystal will not be able to reach the surface 
in a reasonably short time. In this case there will exist a supersaturation 
of vacancies in the crystal, of magnitude c,/c,. If other vacancy sinks 
are present, of course (such as dislocations), they may be able to reduce 
the vacancy supersaturation by absorbing the excess vacancies. 

Frank (1950) was one of the first to describe how vacancies in super- 
saturation in a crystal could condense to form dislocation loops in crystals. 
Hitherto, the most direct support of the mechanism came from the 
experiments of Hirsch et al. (1958), who obtained electron micrographs. 
of polycrystalline aluminium foils 0-075 mm thick quenched from 600°c 
into iced brine and then thinned by electropolishing ; their pictures reveal 
many small dislocation loops. 

It is difficult to determine from thermodynamic considerations the 
vacancy supersaturation needed to nucleate a vacancy disc, for the 
uncertainty in the difference between the free energies of a vacancy disc 
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and a vacancy sphere is as great as or greater than the uncertainty in the 
value of the disc’s surface energy. Once a vacancy disc has formed and 
collapsed to form a dislocation ring (Schoeck and Tiller 1960 show that a 
vacancy disc of radius 5b or more may be considered to have collapsed), 
it is then possible to calculate the supersaturation necessary to keep the 
dislocation climbing. Since all dislocations formed by this mechanism 
must go through this early stage of growth, a lower limit can be obtained 
for the necessary vacancy supersaturation. The supersaturation that 
exists on cooling from a temperature 7’) to a lower temperature 7’, is given 


by 
ik aw Agtl 1 
Cy/¢; =exp = EUR -7) |. . . . . (3) 
: u 0 oil 


The assumption made here is that the equilibrium concentration of 
vacancies with respect to external surfaces, cy, that obtained at the melting 
point 7’, is frozen in; this estimate is, of course, an upper limit. Elbaum 
(1960) used the expression 


c Gb4 r 1 
De yf P36 SS at ge tae | 
"Gy, 4n(1—2vrkT, alg a ) * TEL 6b (4) 


for the supersaturation of vacancies required to nucleate a dislocation 
loop ofradiusr. Gis the shear modulus, vis Poisson’s ratio, kis Boltzmann’s 
constant, b is the Burgers vector and 7’, is the temperature at which a 
concentration of vacancies c, is in thermal equilibrium with the external 
surface; 7, (and thus c,) is a function of 7. Taking into account the 
temperature variation of the dislocation strain energy and expressing r in 
Burgers vectors, Elbaum (1960) finally derived an expression for the 
critical temperature 7’. at which the vacancy supersaturation is great 
enough to form a dislocation loop of radius r: 


3 
i AGP Ta 1 4 300(8 + 30) 

ope NL aE SSE ye 6 05 (6) 
eee re (B+ 3a) 
Ty Us: 
where A is a numerical factor that depends on the value chosen for 7, and 
(Gb?) Tr (which varies approximately as the strain energy per atomic plane 
of a dislocation) is taken at T7,=300°K. £ is the temperature coefficient 
of the shear modulus (which is considered to decrease with temperature) 

and « is the linear coefficient of thermal expansion. : 

Elbaum (1960) made estimates of the maximum size to which a loop can 
grow once it has formed for a given value of the ratio of the temperature 
drop from the melting point, 7)—7'c, to the melting point, Or 
several elements. The maximum size of a loop, rmax, calculated by 
the above method, is so large for the case of aluminium that the actual 
size of dislocation loops originating from collapsed discs will not be 
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determined by rmax- The size will probably depend on the presence of 
other vacancy sinks in the crystal, in particular, on the presence of other 
dislocations. On the other hand, in the case of silicon and germanium, 
the values of rmax are so small that dislocations produced by collapsing 
vacancy discs will not attain sizes visible in an optical microscope. This. 
suggests a reason why Dash (1958, 1959) was able to grow crystals. 
essentially free of detectable dislocations. 

In aluminium, the critical temperature 7’, for a dislocation loop of radius 
5b to grow is 778°K. Substituting this value of 7’c for 7; in eqn. (3) yields 
a minimum vacancy supersaturation in the crystal of 6-6. Allowing for 
errors in the calculations of the minimum size of a dislocation loop (since 
there are uncertainties in the value of the surface energy of a vacancy 
disc), one may say that the minimum supersaturation required to nucleate 
a dislocation by the process of vacancy disc collapse is somewhat less than 
10. The next step in showing that the collapsing vacancy disc mechanism 
is indeed possible is to calculate the vacancy supersaturation that exists. 
in a crystal as it is cooled from the melting point, and thus obtain a lower 
limit to the diameter of a crystal for which no such supersaturation could 
be built up. 

It is therefore desired to study the distribution of vacancies in a 
cylindrical aluminium crystal of radius A, whose temperature is lowered 
from the melting point 7’) to a lower temperature 7’ according to some 
function of time. 

The diffusion coefficient D is assumed to be a function of temperature 
only (and thus of time only) and not a function of concentration (this 
assumption is justified by the small overall concentration of vacancies). 
It is assumed that at any instant of time the temperature is the same 
throughout the crystal; that is, there are no temperature gradients. 
The initial concentration of vacancies in the crystal (when the temperature 
is 7'y) is that which corresponds to equilibrium with the external surface ; 
no other vacancy sources or sinks are assumed present. Thus, 
Cy =exp (— U;/RT 4), where c, is the equilibrium concentration of vacancies 
at 7), Us is the energy for formation of vacancies and R is the universal 
gas constant. As the temperature decreases with time in some fashion, 
the concentration of vacancies at the surface of the crystal is always assumed. 
to be equal to the equilibrium value at the particular instantaneous 
temperature prevailing, or, 


c(A, t)=c4(t) =exp (— U;/RT). ee ea ey 
The concentration of vacancies in the crystal as a function of distance 
from the origin and of time, for a given radius A for the above boundary 
conditions is the particular solution of the differential equation controlling 
the diffusion process. This equation, known as Fick’s second law, may be 
written in cylindrical coordinates : 


dc(r, t) e?c(r, t) in 1 dc(r, t) 


ee SD = 
at ag penser 


O<rsA - (7) 
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subject to the boundary conditions 


(8) 


The most likely functional dependence of temperature on time for crystal 
growth would be a relation of the form 


T =T, exp (— ft) a RO eaccr aey Ye)} 


where 7’, is the melting temperature, ¢ is the time and 8 is a coefficient that 
depends on the cooling rate. It was estimated that in 400 sec the 
temperature of the cylinder drops from 933°K (melting point of aluminium) 
to 588°K, a temperature at which the diffusion rate has dropped to a fairly 
low value; this yields a value of B=1-15 x 10-3 sec. 

An analytical approach to the solution was attempted, but with the 
temperature-time variation of eqn. (9), eqn. (7) could not be integrated 
in closed form. A linear temperature variation of the form 7’=7',—’t 
was tried; eqn. (7) could still not be solved unless it was assumed that 
U;:/Um=2, when it is, in reality 1-75 (Bradshaw and Pearson 1957). In 
view of all the assumptions made, it was decided that the analytical 
solution finally obtained was not a valid representation of the actual 
distribution of vacancies in a crystal cooled from the melting temperature. 
A finite difference solution was then evaluated in which the correct 
temperature-time relation was used with no further assumptions. These 
equations were solved using the [BM 7070 computer. The supersaturation 
of vacancies that exists in crystals of various diameters cooled from the 
melting temperature are plotted in figs. 1, 2 and 3+ as a function of distance 
from the centre of the crystal expressed in fractions of a radius for various 
values of D, (where the diffusion coefficient D=D,exp(—Um/RT) and 
Um is the activation energy for vacancy migration). D, usually lies some- 
where between 0-1 and 0-01; its exact value for aluminium is not known. 
If the minimum supersaturation needed to nucleate a dislocation loop is 
around 10, then one may define a critical crystal diameter (for the cooling 
rate assumed) such that the vacancy supersaturation exceeds 10 if 
d>dcrit where d is the diameter of the crystal. From figs. 1, 2 and 3 one 
may see that derit lies between 0-5 and 4 mm, depending on ike exact value 
of D,. If D, for aluminium is as low as 0-01, derit ¥ 0-5 mm. 

The ane considerations suggest that the mechanism of collapsing 
vacancy discs is possibly of prime importance in producing dislocations 
in aluminium. It is apparent that in a crystal, the density of dislocaitons 
originating from collapsing vacancy discs will be a function of the super- 
saturation of vacancies that exists in the crystal as it is cooled from the 
melting temperature. The vacancy supersaturation depends, as has been 
shown above, on the length of time during which vacancies can move 
over appreciable distances; for high growth rates involving high cooling 


+ Figures 8-12 are shown as plates. 
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SUPERSATURATION 


FRACTION OF RADIUS A 


Supersaturation of vacancies in a cylindrical crystal cooled from 933°K to 
588°K versus fraction of radius. (Radius of crystal in centimetres 
written beside each curve.) D)=1-0 cm? sec t. 


rates, the diffusion coefficient will fall rapidly with falling temperature to 
a low value, preventing further diffusion. A large crystal diameter also 
maintains a vacancy supersaturation by not permitting vacancies to 
escape from the interior to the surface of the crystal before the diffusion 
rate becomes negligible. The following sections describe experiments 
wherein the dislocation structure in aluminium crystals grown from the 
melt by the Czochralski method was studied using an x-ray technique 
that reveals dislocations in the interior of a crystal. 


§ 3. EXPERIMENTAL PROCEDURE 


The experiments consisted of two parts: growing the crystals, and 
examining their structure by means of x-ray diffraction pictures. 
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SUPERSATURATION 


FRACTION OF RADIUS A 


Supersaturation of vacancies in a cylindrical crystal cooled from 933°K to 
588°K versus fraction of radius. (Radius of crystal in centimetres 
written beside each curve.) D)=0-10 cm? sect. 


3.1. Growth of Crystals 


The choice of the particular crystal growing technique used in this 
research was dictated by certain requirements, the most important of 
which was the minimization of deformation due to constraints arising during 


P.M. 4K 
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Fig. 3 


200 


SUPERSATURATION 


FRACTION OF RADIUS A 


Supersaturation of vacancies in a cylindrical crystal cooled from 933°K to 
588°K versus fraction of radius. (Radius of crystal in centimetres 
written beside each curve.) D y=0-01 cm? sec. 


growth, and from handling. This eliminated any form of Bridgman’s 
technique which necessitates solidifying a crystal in a container, and is 
likely to introduce strains due to the constraint of the container walls. 
The Czochralski technique was chosen for several reasons: (1) The crystal 
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need never be touched. (2) Control of orientation is easy. (3) The diameter 
of the crystal can be easily varied. 

In order to study the variation of dislocation density with crystal 
diameter, tapered (carrot-shaped) crystals were grown; this shape could 
be obtained for a given pulling speed by correctly adjusting the melt 
temperature. The range of speeds used was from 0:4 mm/min to 5 mm/min; 
the growth process was performed in a helium atmosphere. The axis of 
the seed crystal was generally parallel to a (110) direction; this permitted 
diffraction from a {200}, {220} and two {111} planes. In some casesarandom 
orientation resulted when the oxide layer prevented contact between the 
molten part of the seed and the melt; in one case a (111) seed crystal was 
used. The purities of aluminium used were 99-999°% and zone-refined 
(12 to 13 passes) in order to study the effect of impurities on the 
dislocation density. In some cases, the specimens were electropolished 
in 95% methanol-5% perchloric acid cooled in a mixture of methanol-dry 
ice. 


3.2. X-ray Diffraction Technique 


An x-ray diffraction technique (Lang 1958, 1959) has been developed 
recently that enables one to see the dislocation structure in the interior 
of a crystal without damaging the specimen. In order to resolve individual 
dislocations by this technique the dislocation density must be of the order 
of 10°cm~ or less. In the present set of experiments, molybdenum 
Ka, radiation at a voltage of 40 kv and an intensity of 20 ma was used. 
Ilford L4 50, nuclear emulsion was employed to record the dislocation 
structure; with this film, the diffracted image could be magnified up to 
80 x without a loss of detail. The dislocations appear in the diffraction 
image as dark areas or lines on a lighter background. This effect arises 
from an enhancement of diffracted intensity from a small volume 
surrounding the dislocation. Dislocation densities were determined by 
counting the number of dislocations in a given region of a crystal and 
dividing this number by the area. The dislocation density was then 
plotted against diameter and rate of growth. 


§ 4. RESULTS 


Figures 4 and 5 show dislocation density plotted against rate of growth 
(with the thickness of the crystal written beside each point) for the two 
purities of aluminium used; both curves show that the dislocation density 
rises with increasing rate of growth. The actual values of the dislocation 
densities are not too meaningful; more significant is the rising trend of the 
curves. There is probably an error of perhaps a factor of two or three in 
the number of dislocations counted, due to clustering (where estimates 
had to be made of the number of discernible dislocations) and poor 
resolution. However, the relative error in counts from one area to another ~ 
of a given crystal or from one crystal to another is quite small, and was 

4K2 


1236 S. Howe and C. Elbaum on the 


Fig. 4 
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Plot of dislocation density versus rate of growth. Purity: 99-999% Al. 
(Diameter of crystal in millimetres written beside each point.) 


Fig. 5 
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Plot of dislocation density versus rate of growth. Purity: zone-refined Al. 
(Diameter of crystal in millimetres written beside each point.) 
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Fig. 6 


Dislocation Density 
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Plot of dislocation density versus crystal diameter. Purity: zone-refined Al. 
Rate of growth: 1:3 mm/min. 


Fig. 7 


‘2 re) 4 
Diameter mm. 


Plot of dislocation density versus crystal diameter. Purity: zone-refined Al. 
Rate of growth: 2-4 mm/min. 
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determined by having different individuals count the same area (this 
error is about 10% or less). 

Figures 6 and 7 show dislocation density versus thickness for two rates 
of growth of zone-refined aluminium crystals. In both cases the dislo- 
cation density increases with increasing crystal diameter in the range of 
0-3 to 0-4mm. For larger diameters, the dislocation density is too great 
to enable one to resolve and count individual dislocations. For smaller 
diameters, the concept of density loses its significance. That is, if a region 
of a crystal contains only several scattered dislocations, the density would 
depend on how one chose the area; the density might have any value from 
zero to a large number in adjacent regions. 

A study of many X-ray micrographs reveals the following pattern: in 
general, the thinnest region of a crystal contains very few dislocations in 
the form of isolated, dark, circular-shaped areas situated near the surface. 
With increasing diameter rather long dislocations appear in the centre. 
Away from the tip towards increasing diameter, their number increases ; 
the dislocations become shorter, spreading out from the centre until, 
in the thickest regions of the crystal, the dislocations appear as dots 
covering the whole area (see fig. 8). One thin crystal (axial orientation 
parallel to a (111) direction) of diameter 0-2 mm containing a few surface 
dislocations (fig. 9) was electropolished. After the treatment, two new 
X-ray micrographs were taken by diffraction from two {220} planes 60° 
apart. Near the tip region in both pictures one sees a dense network of 
dislocations. Towards the tip, the number decreases sharply to a few 
scattered dislocations and drops to zero several millimetres from the tip 
(fig. 10). It was not possible in this case to compare identical locations in 
the crystal in the micrographs taken before and after electropolishing, 
since the amount of material removed is not known. A ‘traverse’ 
micrograph was taken of this particular crystal (for a ‘traverse’, both 
the crystal and the film are translated back and forth across the x-ray 
beam; because the slit has a finite width, this technique will reveal 
individual dislocations only if their number is extremely low. Otherwise, 
overlapping of the dislocation images will result). The traverse shows 
essentially the same features as the section micrographs described above 
(which merely reveal dislocations lying in or intersecting the diffracting 
plane), but establishes most definitely that the tip region is free of 
dislocations detectable by means of this technique; the traverse would 
reveal dislocations situated throughout the volume of the crystal, provided 
the Burgers vector of the dislocation is not parallel to the diffracting plane. 

It is to be pointed out that this is the first recorded instance of a metal 
crystal grown from the melt having a large area free of detectable disloca- 
tions. Another crystal was grown of 99-999°% pure aluminium that 
revealed a tip region containing a few scattered dislocations lying on the 
surface (fig. 11). After electropolishing (fig. 12) it was possible to compare 
similar regions of both micrographs (arrows point to similar regions). 
One sees that some of the surface dislocations have disappeared, and that 
regions are present which contain no dislocations. 
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§ 5. Discuss1Ion AND CONCLUSIONS 


It has been shown above that a necessary condition for vacancies to 
aggregate into discs and collapse is that there be a sufficient super- 
saturation of vacancies in the growing crystal. A minimum value for the 
supersaturation was calculated to be about 10; from calculations of the 
vacancy supersaturations reached in an aluminium crystal cooled from the 
melting point under the conditions of crystal growth described above, a 
critical diameter of about 0-5 mm was deduced, below which no dislocations 
should be produced by the collapsing vacancy disc mechanism. This is 
of the same magnitude as diameters of crystals for which the density 
of detectable dislocations is rapidly decreasing with decreasing diameter 
(figs. 6 and 7). The fact that the observed dislocation densities (figs. 4 
and 5) increase with increasing rates of growth is consistent with the 
theory presented ; a high growth rate implies that vacancies present in the 
interior of a crystal would not have time to diffuse out to the surface before 
the diffusion coefficient drops to a low enough value for the rate of diffusion 
to become negligible. It must be pointed out that the results of the 
experiments support the validity of the mechanism, but do not deny 
that other ways of forming dislocations either may be operating simul- 
taneously or would operate under other growth conditions. We will now 
consider the operation of the other mechanisms mentioned in the 
introduction. 

(1) Propagation of dislocations into the growing crystal from the seed 
crystal: It is apparent from fig. 11 (and from other x-ray micrographs 
not included in this paper) that the dislocation density does not change 
suddenly with a change in crystal diameter. This is consistent with the 
idea of dislocations propagating directly or growing by climb from one 
region of a crystal to another. Ifa dislocation line is at an angle to the 
crystal axis, it will eventually grow out to the surface of the crystal; only 
dislocations parallel to the axis of the crystal will continue growing 
indefinitely (even these dislocations can in principle climb and ultimately 
reach the surface of the crystal). This explains why adjacent regions of 
a crystal of the same diameter may have varying densities of dislocations ; 
the last region to solidify will ‘inherit’ some of the dislocations from a 
previously solidified region. The propagation of dislocations from one 
region of a crystal to another is also an explanation of why the observed’ 
critical diameter (that diameter below which the vacancy supersaturation 
reached during growth is too small to nucleate dislocations by the 
mechanism of collapsing vacancy discs) is smaller than the theoretical 
value. Thus, a particular region of a crystal contains dislocations arising 
from two sources: those inherited from a previously grown part of the 
crystal and those created from collapsing vacancy discs. 

(2) Shear strains arising from thermal strains or external mechanical 
constraints: Thermal and mechanical strains were probably reduced to 
a minimum by the use of the Czochralski technique, but the many 
dislocations that are situated near the surface may have been produced 
by mechanical stresses present during growth. These stresses might 
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arise from drag against the liquid of the oxide layer attached to the 
rotating crystal during growth. A torsional stress, which reaches a 
maximum at the surface of the crystal is thus introduced. The fact that 
these dislocations disappear after electropolishing proves that they lie 
on the surface. 

(3) Shear strains resulting from abrupt changes in composition accom- 
panied by changes in lattice parameter: Qualitatively, regions involving 
large impurity content differences are less likely to occur in crystals grown 
from zone-refined (high purity) materials. In these experiments, no 
significant difference was observed between values of the dislocation 
density in 99-999% aluminium or zone-refined aluminium. It thus 
appears, that this mechanism was not playing an important role. 

(4) Dendritic growth: The crystals were grown by the Czochralski 
method, and did not solidify dendritically. 

In conclusion, the results of the experiments performed in the course 
of the research are completely consistent with the idea that the mechanism 
of collapsing vacancy discs is responsible for creating most of the 
dislocations found in a metal crystal grown from the melt. It was also 
found that the crystal diameter for which no dislocations were detected, 
was somewhat less than the theoretical critical value; this is interpreted 
to mean that for a given vacancy supersaturation, the density of dislo- 
cations present contains two contributions: those created from collapsing 
vacancy discs and those inherited from the previously solidified part of 
the crystal. 

Extensive regions of thin crystals were found to be free of dislocations 
when the diameter of the specimen was about 0-3 mm or less; this is the 
first recorded instance of the production of dislocation-free metal crystals 
by growth from the melt. 
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ABSTRACT 


Each magnetic grain in a rock contains ferromagnetic domains with a 
range of sizes, the small ones being magnetically harder and responding less 
readily to demagnetizing fields than the larger ones. By exposing a grain 
which has a thermoremanent magnetic moment to an alternating field, domains 
larger than a critical size are demagnetized, leaving the smaller domains un- 
affected. Domain size distributions are calculated for spherical and cubical 
grains and used to derive characteristic curves of the fraction of an initial 
thermoremanent magnetization which remains after demagnetization in a 
progressively increasing field. Coercive force, saturation remanence and 
anhysteretic magnetization are also calculated from the domain size 
distributions. Irregularities of grain shape are shown to lead to a pseudo- 
single domain behaviour which is responsible both for very hard components 
of thermoremanent moments and for the fact that small, randomly directed 
moments may be found to remain after rocks have been demagnetized in the 
highest fields. 


§ 1. INTRODUCTION 


THe technique of alternating field demagnetization of rocks is now an 
essential part of palaeomagnetic studies in a number of laboratories. Its 
use is based upon the observation that primary and secondary com- 
ponents in the natural remanence of rocks can often be distinguished. 
In most igneous rocks the secondary components are probably induced 
isothermally over long periods and are more readily destroyed than the 
original or primary thermoremanence (see for example Irving e¢f al. 
1961). In any physical description of the demagnetization process the 
important fact requiring explanation is that thermoremanent magnetiza- 
tion (TRM) is very resistant to demagnetization. 

This paper extends the multidomain theory of magnetic grains in 
igneous rocks (Stacey 1958, 1959) to give a semi-quantitative description 
of the physical process of demagnetization of TRM. It may be regarded 
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as a third-order calculation in the multidomain theory. To a useful 
approximation thermoremanent magnetization is a first-order effect, 
depending upon grain shape and the existence of a blocking temperature 
(Stacey 1958); grain size and domain structure have only secondary 
influences, provided the grains are true multidomains (as defined in 
§ 6). The variation of domain structure with grain size is a second-order 
calculation and coercive force is a property described by such a calcula- 
tion (Stacey 1959). To proceed further, specific consideration must be 
given to grain shape and domain structure, and in particular to the 
distribution of domain sizes within each grain. It is necessary for the 
explanation of alternating field demagnetization that the several domains 
in a single grain do not all have the same properties. 

The detailed magnetic behaviour of rocks in alternating fields is so 
sensitive to grain shape that there is no possibility of deriving theoretically 
a demagnetization curve which fits the measured curve for the TRM of 
any rock. Magnetic grains are irregular in shape, whereas only simple 
geometrical shapes can be considered theoretically, even approximately. 
The calculations therefore indicate the broad features of demagnetization 
and other alternating field effects in terms of idealized spherical or 
cubical grains. 

Irregularities in grain shape are particularly important when the 
grain size is small, i.e. less than about 1, but larger than 0-03 », which 
is roughly the upper limit for single domain grains. In this case magnetic 
properties may be substantially different from true multidomain behaviour 
and more like a superposition of single-domain properties upon the 
multidomain character. The essential feature is the appearance of 
permanent dipole moments which can be reversed but not destroyed by 
demagnetization, and the effect is therefore termed pseudo-single- 
domain behaviour (§ 6), which can occur in grains which are much too 
large to be true single domains. 

The inapplicability of the simple multidomain theory to grains smaller 
than about 1 was pointed out by Stacey (1958, §7(5)). However, 
Verhoogen (1959) has questioned the fundamental basis of the multi- 
domain theory on the grounds that some rocks are found to have single- 
domain-type properties; he suggested that the freedom of domain wall 
movements in a multidomain grain above its blocking temperature 
implies that energy barriers opposing domain wall motion must vanish 
and therefore that magnetocrystalline anisotropy and magnetostrictive 
energies also vanish, whereas it merely means that thermal activation 
provides sufficient energy to impel walls past the barriers (Stacey 1960 a). 
Verhoogen attributed the very hard components in the TRM of some 
rocks to the presence within normal, multidomain-sized grains of small, 
single-domain inclusions, which are chemically and physically indis- 
tinguishable parts of the grains except for the presence of very strong, 
local stress concentrations which cause them to act independently as 
regards their magnetic properties. Even if this is possible, the fact that 
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the single domains would be completely embedded in the magnetically 
softer matrix would neutralize their properties. The single domains 
would be magnetically screened by the matrix so that TRM would still 
be acquired by grains as a whole. Further, if alignment of the moments 
of the single domains in a grain were produced (the application of a strong 
field appears to be the only method) then the aligned moments would 
give rise to a self-demagnetizing field in the grain as a whole and so 
would tend to magnetize the matrix in the opposite direction. During 
demagnetization the domain structure of the matrix would adjust itself 
to produce zero internal field, i.e. zero total moment, without affecting 
the inclusions. Thus the remanence and magnetic hardness of the 
single-domain volumes would not be apparent. These difficulties are 
only partly avoided by supposing that the single domains represent a 
phase separation of finely divided but strongly magnetic material in a 
less magnetic matrix, as in dispersion-hardened alloys (Stoner and 
Wohlfarth 1948), and other difficulties are introduced. The hard com- 
ponents of TRM, which Verhoogen considered, are readily explained in 
terms of the theory of very small multidomain grains of irregular shape 
(§ 6). 

In this paper consideration is given only to the properties of strongly 
magnetic minerals, magnetite and the titanomagnetites, which normally 
occur in igneous rocks. Haematite, which is frequently a constituent of 
sedimentary rocks is quite different and cannot be treated in the same 
way. 


§ 2. Domain StRUCTURES OF SPHERICAL AND CUBICAL GRAINS 


The domain structure of a magnetic grain in its lowest energy state 
in zero external field can, in principle, be determined by minimizing the 
total magnetic energy. The more important energy contributions which 
must be balanced in such a calculation are magnetostatic (due to free 
poles), magnetocrystalline, magnetostrictive and domain wall energies. 
In practice the simultaneous consideration of all four energies is impos- 
sible, so that a probable model of the domain structure is first set up and 
can then be adjusted to minimize the energy contributions in pairs. 
Stacey (1959) showed that if the domain structures of magnetite grains 
are assumed to be controlled mainly by magnetostrictive and domain 
wall energies, then the domains take the form of rods of square cross _ 
sections. For irregular or spherical grains the magnetostatic and 
anisotropy energies must be appreciable but in view of the agreement 

between calculated and observed coercive forces (Stacey 1959, fig. 1) it 
is reasonable to assume that the rod model is nevertheless applicable to 
magnetite grains. 


2.1. Spherical Grains 


In spherical grains the rod model takes the form shown in fig. 1 (a), 
the spacing of domain walls being determined by balancing wall energy 
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with the energies of domain end-effects. A grain in the demagnetized 
state has no overall self-demagnetizing field, so that one part of the 
grain can have no influence on the domain structure of any other part 


Fig. 1 


yl 


(a) Rod model of the domain structure of a spherical grain of magnetite in the 
virgin state. 


(a) 


(6) Regions of net magnetization in a grain which has been magnetized to 
saturation in the upward direction and then placed in a downward 
directed field equal to the coercive force, producing zero total moment. 


more than one or two domain-widths away. Thus the energy balance 
of each part of the grain is locally adjusted, the widths of the domains 
being determined by their lengths, as in the original rod model calculation 
(Stacey 1959), the result of which may be used directly. The critical 
field, Her, required to move a domain wall past potential barriers whose 
heights are characteristic of the length of the wall, l, is given by 


Hep SC BRB yc 5 in ch aie a 


Compared with the central part of a grain the magnetic energy balance 
causes the domains of the outer parts to be narrower and magnetically 
harder, as a consequence of their shortness. Within each grain there is 
therefore a spectrum of critical fields controlled by (1). In a grain of 
diameter d and volume V, there is a volume v which contains domains 
shorter than /, where 

v = (77/6)1, 
so that 

vf Vs (t/a) 8 0 | Pe ee kee ee 


If the largest, central domains of length d have a critical field Her 0, 
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the spectrum of critical fields is obtained by substituting (1) in (2) from 
which 


V Her 9 
Equation (3) gives the volume fraction v/V of spherical grains for which 
the critical field is greater than Hey. For Helicrg= lh, vj V=1. 
Equation (3) is the starting point of the demagnetization calculations 
(§ 3). 

As a consequence of the domain structure indicated in fig. 1 (a), the 
central region of a spherical grain can be magnetized (and demagnetized) 
more easily than the outer parts. After magnetization to saturation and 
the application of a smaller field in the opposite direction, as in the 
measurement of coercive force, the grain therefore contains regions of 
net magnetization as shown in fig. 1 (b), in which the downward directed 
domains are larger than the upward directed domains in the central 
region and the upward directed domains are larger in the outer parts. 
However, this domain structure does not affect the variation of coercive 
force with grain size. 


H —4 
7-( = | Roe ee ee? os. (3) 


2.2. Cubical Grains 


The interest in grains of cubical shape is not that they represent any 
better approximation to the actual shapes of naturally occurring magnetic 
grains than do the spheres considered above but that they indicate very 
clearly the strong dependence of certain magnetic properties upon grain 
shape and even upon the orientation of the domains with respect to the 
shape. In a cube there are an infinite number of possible orientations 
the three principal cases being: 

(i) domains parallel to cube edges, 
(ii) domains parallel to face diagonal, 
(iii) domains parallel to body diagonal, 
Consideration of the geometry of each case gives the distribution of 
domain lengths, from which the spectrum of critical fields is obtained. 
These spectra are: : 
Case (i) all domains have the same critical field, 
@asev (iy iVee= egiltepjeo my as aoe ee (4) 
Case (iii) o/V= (Her/Her e243 ra 2(H er/Her Seen): . 5 Fy (5) 


§ 3. ALTERNATING FreELD DEMAGNETIZATION OF TRM IN SPHERICAL AND 
CuUBICAL GRAINS 


In this section it is assumed that each grain considered is demagnetized 
by the application of an alternating field in a direction parallel to its 
domains (which is the most effective direction). A rock specimen con- 
taining magnetic grains in which the crystallographic axes and domain 
directions are randomly oriented can be demagnetized by rotating it 
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about two mutually perpendicular axes simultaneously, at frequencies 
unrelated to one another and to the field (as in the experiments of Irving 
et al. 1961). In this way the grains are presented to the field in all 
possible orientations and the field therefore acts on each grain in the 
direction of its domains or very near to it. The calculations therefore 
refer to this method of demagnetization. Certain effects arising from 
the demagnetization of stationary specimens are considered in § 5. 


3.1. Spherical Grains 


When a spherical grain is given a small TRM the domains in the 
direction of net magnetization are all enlarged at the expense of their 
reversed neighbours, so that the net magnetization is uniformly distri- 
buted through the grain. Subjection to an increasing alternating field 


Partial demagnetization of a spherical grain. Heavy arrows in outer region 
indicate remaining TRM, light arrows indicate alternating magnetization 
induced by the applied field and show the form of its variation within 
the grain. The heavy arrows are not to scale as TRM is normally much 
smaller than the induced magnetization under these conditions. 


induces no irreversible domain wall movements until the effective field 
acting on the grain reaches Hery, when it becomes just sufficient to 
demagnetize the central region, i.e. to rearrange the domain walls, making 
the net magnetization of the region zero. With further increase in the 
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applied field to a peak value H the demagnetized region spreads outwards 
to a radius r, at which the critical field required to move the domain 
walls Hey’ is equal to the effective field acting at r (applied field, H, less 
self-demagnetizing fields due to polarization of the grain). 

A grain in such a partially demagnetized state is shown in fig. 2, in 
which the remaining TRM of the grain is carried by the volume outside 


radius 7, where 
Qp\ 2 y’\ 2/8 
rae) We 2g pee 5 
phe Obern Cuse. 


d being the grain diameter and V its total volume as before. Since the 
original TRM was uniformly distributed through the grain the fraction 
remaining, M/M,, is given by 


MM FOV, 5 le OR eC) 


The induced magnetization, indicated by the light arrows in fig. 2, is 
due both to reversible and to irreversible movements of domain walls. 
Thus the field He,’ acting at 7 is given by 


Te aReON GH = FP oe 5) ee a8) 


In eqn. (8) x is the low-field intrinsic susceptibility, i.e. that contribution 
to susceptibility which is due to small displacements of domain walls, 
which do not take them across potential barriers and do not produce 
remanence, and N is the demagnetizing factor appropriate to the reversibly 
induced moment, (yHer’) being the reversible moment. Although not 
strictly correct it is assumed for simplicity that the reversible moment 
is uniform over the whole grain, so that N=4{7 at all points. H), is 
that part of the self-demagnetizing field which is due to polarization of 
the central (demagnetized) region by irreversible movements of domain 
walls. Since, at any point inside radius r the critical field Hey is less than 
the effective field H,,’ acting at r, the irreversible polarization is con- 
trolled by (Her’—Her) at each point, and is superimposed upon the 
reversible polarization. This means simply that the irreversible polariza- 
tion has just such a magnitude that the self-demagnetizing field reduces 
the effective field at each point in the central region to the value of Hey 
at that point, this being the maximum field which the domain walls can 
withstand. Thus the field H,, acting at radius r is due to an alternating 
moment M’ in the central region, being given by ; 


y Her’ — Her 
Pass Vere CD) ere. ay ee oak (ee 
M'=[" S25 F2 (— dy (9) 
N’ being the demagnetizing factor of the central region, which is for 
this purpose adequately regarded as a uniformly magnetized prolate 


ellipsoid of length d and width 2r. Substituting for Hor in terms of v 
from (3) we obtain an expression for M’ which integrates directly to give 


4‘\ —1/4 4\ 3/4 
we “=| (5) +3(5) -4]. ee eet 0 
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Then the field H, due to M’ is given approximately by 


M’ 
= —__—__;; Pe ets Gree lg 
Hy ((d/2)? + 1?)3/2 ( ) 
Substituting for r from (6) and noting that V =$7n(d/2), we obtain 
fai, 47 M (12) 


Equation (8) may be rewritten 
1 Hy Hest 1 Hy 

LEN y | H@oee Hey 1+Ny Herg 
Now (Her’/Her 9) =(v'/V)-™4 from (3) and Hy is given by (12) where M’ 
has the form (10). Making these substitutions an expression is obtained 
for H/Hey » in terms of (v'/V), which is then expressed in terms of (M/J)) 
according to (7). Thus the final equation for the demagnetization of 
TRM in spherical grains is 


ee 
LEN a M, 3D'(1+ Nx) © (2—(M/M,)27]3?2 


(13) 
where, for convenience D'=N’'/47. Equation (13) gives the peak value 
HI of the alternating field required to reduce the TRM of a spherical 
grain to a fraction M/M, of its initial value, in terms of the constant 
Hey, and the intrinsic low field susceptibility x. It is convenient to 
consider the characteristic field for alternating field demagnetization to 
be 

Hy=Her » (1+ Ny) Sas nae Pa 
which is the field below which no demagnetization occurs. 

The two limiting forms of (13) are plotted in fig. 3, curves 2 and 3, 
The procedure for plotting curve 3 was to select a series of values of 
(d/2r) for which corresponding values of D’ have been tabulated by 
Stoner (1945). Values of (v’/V) and hence (17/M,) are given by (6) and 
used directly in (13) to find (H/H,). Physically the two limits repre- 
sented by these curves are that for very large y the self-demagnetization 
due to irreversibly induced moments is negligible compared with that 
due to the reversible moments and for very small y the reversible moments 
are negligible and the irreversible moments become important. The 
intrinsic susceptibility y, used in this section is the susceptibility in the 
direction of the domains, for which the symbol X, is used in §5. yx, is 
related in an inverse manner to Hey, and is hence not an independent 
constant, although no quantitative relationship between them is known. 


Both yx, and Her» are also related to the coercive force He, which is 
considered in § 4. 


3.2. Cubical Grains 
Demagnetization curves for cubical grains with y-> oo (negligible self- 
demagnetizing effects due to irreversible domain wall movements) follow 
directly from eqns. (4) and (5), if it is assumed that TRM is uniformly 


Magnetic Properties of Igneous Rocks in Alternating Fields 1249 


distributed through a cube, since in this case v/V=M/M,, as for 
spherical grains. This is the approximation used to plot curves 4 and 5 
of fig. 3. Curve 1 is trivial since it assumes simply that all the domains 
are alike. Equation (7) is not exactly valid for cubes because the factor 


Fig. 3 
Re) 
o-8 
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O-4 
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1 2 4 
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Alternating field demagnetization. Theoretical curves for (1) cube with domains 
parallel to cube edge; (2) sphere, assuming high intrinsic, low-field 
susceptibility, x; (3) sphere, assuming y+0; (4) cube with domains 
parallel to a face diagonal, y~°0; (5) cube with domains parallel to a 
body diagonal, y00. 


N which controls TRM (Stacey 1958) is not constant over the whole 
volume and TRM is therefore not uniformly distributed, but (7) is a_ 
useful approximation. Curves for cubes with y-+0, in which the self- 
demagnetizing fields lead to terms corresponding to the second one in 
(13) have not been plotted on fig. 3. The curves are readily obtained by 
approximating the central demagnetized regions to ellipsoids, as in the 
case of spherical grains, and lie above the curves for y— 00. However, 
the nature of the approximation is such that this result is only of qualita- 
tive interest. 

The variety of shapes of the curves of fig. 3 shows clearly that de- 
magnetization is very dependent upon grain shape. Demagnetization 
of the irregular grains occurring in rocks can therefore not be calculated 
quantitatively but the physical principles must be the same as for the 


P.M. 4L 
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simple shapes and the general form of curves of demagnetization of 
TRM must be similar. The more irregular a grain is, the less sharp 
becomes its demagnetization curve. Observed demagnetization curves 
are spread out by several other factors including grain size distributions 
which have the effect of superimposing curves with different characteristic 
fields H,. Many rocks contain intergrown and exsolved minerals which 
are even more complicated and can only be discussed qualitatively. 


§4. ComrcrveE Force, SaruRaATION REMANENCE AND ANHYSTERETIC 
MAGNETIZATION OF SPHERICAL GRAINS 


4.1. Coercive Force 


The measurement of the coercive force of a single grain, by the applica- 
tion of fields parallel to its domain directions, produces the state of 
magnetization indicated in fig. 1 (b), the total moment being zero, so that 


si lads i eee eee 


V 


where J is the magnetization at any point. Since the grain is subjected 
to a downwards directed field He, the effective field at any point is 


He=Ho—N f, . . . . . . . (16) 


N’ being the demagnetizing factor appropriate to the region in question. 
Equation (16) expresses the fact that at each point the domains are 
supporting a field equal to the critical field H.,y. Then 


T=(Hg- HaliN’ 0s ve eee 


It is difficult to specify an appropriate value for N’, but (15) can be 
integrated by assuming that N’ is effectively constant over the whole 
grain. Substituting for H,; from (17) and (3) and integrating (15) we | 


obtain 
- [a2 (7)-8 v\s/47 ae 
N' Her 9 V = V |_,- : 
He =H ¢ o. : . . . . . : (18) 


When an assembly of grains, uniform in size but with randomly 
oriented crystallographic axes is considered, an additional numerical 
factor is included. It is assumed that during the measurement of 
coercive force the domains of each grain are re-aligned to the easy direc- 
tion of magnetization which is nearest to the direction of the applied 
field. Consider a grain in which the domains make an angle « to the 


direction of the field H., equal to the coercive force of the whole assembly. 
Its magnetic moment will be 


M = B(H, cos «—$Her y), - <2 eee V0) 


where 6 depends upon grain volume, intrinsic susceptibility and the 


whence 
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demagnetizing factor, but is an arbitrary constant for the present purpose. 
The component of M in the field direction is 


M cos «=B(H, cos*«—4He, cosa), . . . . (20) 
but for the whole assembly 
M cos «=0, 
so that from (20), 
H.=$Her 9 ss (21) 
cos? a 


Taking values of cos « and cos? « from the Appendix, we obtain 
He=V-57 Hes . . . . . 
Slightly different factors would apply to other grain shapes. 


4.2. Saturation Remanence 


For a single grain, the saturation remanence per unit volume is 


1 f° 
I= 7 [ Tie de\e hn Se See (on 
where, at any point, 
hi — Hey|N. . . . . . . . . (24) 
Substituting for He; as before and integrating, we obtain 


iT, 
Ipg=4 ee 


For an assembly of grains the saturation remanence of an individual 
grain must be multiplied by the mean cosine of the angle between the 


individual moments and the total moment, cos «. Assuming these to 
be the minimum possible angles, since the grains have been exposed to a 


strong field to saturate them, cos « is given in the Appendix, so that 


Bee H 

Iaa=% 008 0 eo =0-72657. . . . (25) 
Equation (25) thus expresses the saturation remanence of an assembly of 
grains of uniform size in terms of the coercive force of the whole assembly. 


4.3. Anhysteretic Magnetization 


Anhysteretic magnetization is induced in a rock when it is exposed to 
a small, steady field h, superimposed upon an alternating field H which 
is gradually reduced to zero. The interesting case is saturation 
anhysteretic magnetization, induced when the alternating field is initially 
strong enough to saturate the grains, and is applied parallel toh. Figure 2 
may be used to represent a grain undergoing anhysteretic magnetization, 
when H has been reduced to the point at which irreversible domain wall 
movements no longer occur outside radius 7. The outer domains then 
carry an anhysteretically induced moment represented by the heavy 


4L2 
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arrows. The effective field acting at radius r oscillates between the 
values 

(1 +h)—Hp+ and (H —h)—Hp-_ 

1+ Nx 1+Nx : 
where H,, and Hy _ are the self-demagnetizing fields due to the moments 
induced irreversibly by external fields (H+h). The net effective steady 
field acting at r is thus 
1 

‘=4,—_— (2h-(Hy, -A se hie Sas eee 

=} ay (2h (Hp Hy) (26) 
The self-demagnetizing fields may be treated exactly as in §3, and 
examination of the numerical data shows that for spherical grains 


dH,/dH is approximately constant over the whole range of interest and 
equal to about 0-2, so that 


(Hp, —Hp_) = 9-2(2h), 


whence 
h’= 0-8h/(1+ Nyx). wend ee Rais Dome eae 
Thus the saturation anhysteretic magnetic moment per unit volume is 
Lip Sl [Ngo Ss te een 


where h’ is given by (27). It is therefore significantly smaller than the 
TRM induced in the same field h. 


§ 5. ANISOTROPY OF SUSCEPTIBILITY INDUCED BY ALTERNATING FIELDS 


Intrinsic magnetic anisotropy of rocks is due to alignment of grain 
elongations (shape anisotropy) and crystallographic axes (magneto- 
crystalline anisotropy); as a special case of shape anisotropy, fabric 
effects due to the stringing together of magnetic grains in lines or planes 
may also be considered. Stacey (1960 b) showed that such intrinsic 
anisotropies are conveniently measured by means of a torquemeter, with 
which it is possible to distinguish the types of anisotropy. Graham 
(1954) developed a sensitive apparatus for measuring anisotropy of 
magnetic susceptibility which is an alternative method of examining 
magnetic anisotropy and has been used by Balsley and Buddington 
(1960) and Fuller (1960). An important difference is that the suscepti- 
bility method is sensitive to induced anisotropy as well as to the intrinsic 
anisotropies. This is a disadvantage when slight intrinsic anisotropies 
are being measured; however it appears that useful information about 
the fundamental magnetic processes occurring in rocks can be obtained 
from measurements of induced anisotropy. 

Induced anisotropy results from the alignment of magnetic domains 
and is caused by the exposure of rock specimens to magnetic fields and 
possibly to several other effects including the application of stresses. It 
arises by virtue of the differences between the magnetic properties of 
grains parallel and perpendicular to their domain directions. A similar 


Magnetic Properties of Igneous Rocks in Alternating Fields 1253 


effect has been observed in bulk metals by Jackson et al. (1958), whose 
results give a useful qualitative agreement with the explanation here 
advanced more specifically for rocks. 

The intrinsic low-field susceptibility of the material of a grain measured 
parallel to its domain directions, y), is a measure of the ease with which 
the 180° domain walls may move. Such domain wall movements allow 
the growth of favourably oriented domains at the expense of their 
oppositely magnetized neighbours; in sufficiently low fields the movements 
are reversible and do not take the walls past potential barriers, but 
simply displace them from the centres of the potential wells between the 
barriers. x, is inversely related to the barrier heights and therefore to 
coercive force, but the relationship is not known quantitatively. The 
susceptibility measured perpendicular to the domain directions, x ,, is 
due to quite a different magnetization process. The application of a 
small field H perpendicular to the domains causes the magnetization of 
the domains to turn towards the field direction by an angle 6 while the 
positions of the walls remain stationary. This gives a magnetization 
I, sin @ perpendicular to the walls, [, being the saturation magnetization, 
and causes the appearance of magnetocrystalline energy by virtue of the 
deflection of the domains out of the easy directions. The total energy is 


E=|2K,+2K,|sin?9—HI,sin6, . . . . (29) 


K, and K, being the magnetocrystalline anisotropy constants. Since @ is 
determined by the condition dH/d@=0, this gives 


x I sin 6 poe thee 
UH RFK 
Normally y, is somewhat larger than y ,. 


In an assembly of grains with domain directions randomly oriented 
the effective intrinsic susceptibility is 


GaN aes Oy NO CE Ce ee LS) 
being the same in all directions. The effective susceptibility is changed 
in value and becomes anisotropic when a stationary rock specimen in 
which there is no alignment of grain elongations or crystallographic axes, 
is subjected to an alternating field which is gradually reduced to zero 
from a value high enough to saturate the specimen. After this treatment 
the domains of each grain lie in the nearest easy direction of magnetiza- 
tion to the direction in which the field was applied. The susceptibility 
in the direction of the alternating field becomes 


(30) 


X= Xj C08? | +X (1—cos* my), + + - + (32) 
and the susceptibility in any perpendicular direction is 


Xp = ${xj(1 — cos? aGy+x,(1+cos? mJ}, 2. 2: (33) 


where cos? «, is the mean square direction cosine of the directions of the 
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domains with respect to the direction in which the alternating field was 
applied. 

In magnetite the easy directions of magnetization are [111] crystal 
axes, and cos? «,=0-728 (see Appendix), so that (32) and (33) give the 
induced anisotropy of susceptibility y,/y, in terms of the fundamental 
ratio x,/x 1 > 

Xi _ 0°728 xy/x 1 +0-272 
Xp 0°136 y/x , +0862" 


Then using (31) and (34) absolute values of x, and x, may be calculated 
from the absolute intrinsic susceptibility y of the isotropic rock and 
the induced anisotropy ratio y,/y,. Thus the factor I,?/|$K, + ¢K 2| can be 
estimated in terms of susceptibility measurements only, or alternatively 
a known value of this factor may be used to estimate the absolute 
intrinsic susceptibility from the anisotropy ratio. 

In measuring the intrinsic susceptibility of the material comprising the 
magnetic grains of a rock, account must be taken of the demagnetizing 
factors of the grains. The susceptibilities considered above are all 
intrinsic to the magnetic material and are related to the observed suscepti- 
bilities by the usual expression 


(34) 


UX 

X04 Nyse 
where yx, is the observed susceptibility corresponding to an intrinsic 
susceptibility y, when the volume fraction of magnetic mineral present 
is v, which is assumed to be small so that grain interactions can be 
neglected. WN is the effective demagnetizing factor. For an assembly of 
spheres, N =$z but in rocks the grains are better represented by ellip- 
soids. Stacey et al. (1960) suggested that in igneous rocks a good 
approximation would be general ellipsoids with dimension ratios 
1-65 : 1-45: 1, although the middle figure is not very reliable. As the 
nearest simple case we may take an isotropic assembly of grains to be 
equivalent to an assembly of oblate ellipsoids with a dimension ratio 
b:a=1-5: 1 and randomly oriented axes. The demagnetizing factors 
are then N,=47 (0-466) and N,=47 (0-267) and as an approximation 
1/N=1/N,+2/N, so that an effective value of N is 3-91 compared: with 
4-19 for an assembly of spheres. It is doubtful whether this difference 
is significant in view of other uncertainties in determining intrinsic 
susceptibilities. 


§ 6. PSEUDO-SINGLE-DOMAIN BEHAVIOUR OF FINE GRAINS OF 
TRREGULAR SHAPES 


In the foregoing sections consideration is given only to spherical or 
other simple geometrical grains, which are much larger than the critical 
size for single-domain behaviour and have perfect multidomain properties. 
A perfect multidomain grain may be defined as a grain whose lowest 
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magnetic energy state in zero field is a state of zero total moment. Perfect 
multidomains may still have remanence but are not then in their lowest 
energy states. Naturally occurring magnetic grains are not regular in 
shape and frequently a single rock sample contains grains with a wide 
range of sizes. The irregularities of shape are particularly important 
in the smaller grains because they result in departures from perfect multi- 
domain character, i.e. the lowest magnetic energy states of the grains in 
zero external field are states of finite magnetic moments. Such grains 
are therefore considered to have pseudo-single-domain properties super- 
imposed upon their multidomain characters. 

A simple example is the wedge-shaped two-domain grain, shown in 
fig. 4, which is too large to be a single domain but only large enough to 


Fig. 4 


Magnetization in a wedge-shaped two-domain grain. 


have two domains with moments M, and M, as shown. The position of 
minimum energy of the domain wall is determined by the balance of wall 
energy and magnetostatic energy due to surfaces of polarity at the ends 
of the domains. Qualitative considerations are sufficient to show that 
the minimum energy occurs at a position such that M,>M,. The domain 
wall energy, being proportional to wall area, decreases as the wall moves 
to the left in the figure. The magnetostatic energy also has a minimum 
left of the position at which M,=,, because the shorter domain on the 
left experiences a stronger self-demagnetizing field and therefore has a 
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higher magnetostatic energy than the domain on the right; this energy 
is reduced by a limited leftward displacement of the wall. The effect is 
to give the grain a permanent dipole moment (J/,—M,) which is respons- 
ible for pseudo-single-domain properties. While this example is 
particularly simple, it is evident that permanent dipole moments must 
be normal in very small grains and may even be due to small irregular 
protuberances on larger grains; pseudo-single-domain properties, super- 
imposed upon multidomain properties, must be expected in rocks in 
which an appreciable fraction of the magnetic minerals occur as grains 
smaller than about 1 micron in diameter. The permanent dipole moments 
are necessarily very much smaller than the saturation moments of the 
grains, but can have an important effect on TRM. 

The fact that single-domain-type properties can be due to grains as 
large as one micron, while true single domains (in which the permanent 
moments are equal to the saturation moments) are smaller than about 
0-03 micron, is of considerable interest in rock magnetism. Verhoogen 
(1959) has pointed out that a compromise between the single-domain 
and multidomain theories appears to be necessary to explain the properties 
of some rocks; the data quoted in his paper is used here to support the 
pseudo-single-domain approach to this problem. For convenience in 
comparing it with the data, the multidomain theory of TRM (Stacey 
1958) may be rewritten as follows, using the original notation: 

TRM per unit volume of magnetic mineral is given by 


Tyee Nei ee ee 


The first term of (36) is equivalent to Verhoogen’s eqn. (5) except that 
in the latter the factor J,)/J, was accidentally inverted. The second 
term is the susceptibility correction (Stacey 1958, § 4), which allows for 
the fact that with magnetization J, a grain experiences a self- 
demagnetizing field (—NJ,) which produces a reversed moment. y is 
the intrinsic susceptibility, which is generally unknown, but 


I)(1+Nx) 


= > 
xy, 


(37) 
Q being the ratio of TRM acquired in a small field to the moment induced 
at room temperature in the same field. This is a quantity which is 


frequently measured in experimental work on rock magnetism. Combin- 
ing (36) and (37) we obtain 


faites AMEE SI 
N op Q+og/op 
The data in Verhoogen’s table 1 may be used as a test of the reliability 


of (38) and the relevant figures are given in the table below, together 
with the estimates of TRM from (38). The ratio (cp/o,) of spontaneous 


(38) 
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magnetization at room temperature to the value at the blocking tempera- 
ture is generally about 3, but for the material in column 4 of the table, 
the Curie point is so low that (cg/c,) could not be much greater than 2, 
as noted in the table. 


Properties of powdered magnetic minerals in controlled grain sizes. 
Figures quoted by Verhoogen (1959) from Roquet (1954) and Uyeda (1958), 
with calculated values of TRM by multidomain theory 


ee roost Nee reo. 


Material Fe,0, 0-89 Fe,0, 0-33 Fe,0, 
Grain size (u) 0-1 100 10 
Curie point (°c) 582 550 123 
TRM in 2 oersteds by experi- 

ment (e.m.u./em?) ol On77, 0-69 
@ (eqn. (24)) 17 2-2 4:8 
Assumed (cR/op) 3 3 2 
TRM in 2 oersteds from eqn. 

(25) (e.m.u./em) 1-2 0-61 0-70 


The TRM of the large grains listed in columns 3 and 4 of the table 
agrees excellently with (38); the uncertainty in estimating (c,/c,) is 
sufficient to explain the small differences. The very small grains of 
column 2 behave quite differently, as would be expected from their small 
size. The pseudo-single-domain nature of the TRM of the very fine 
material. is confirmed by Roquet’s (1954) observation that its TRM 
was not proportional to field but approached saturation in quite low 
fields. When an assembly of grains with permanent moments ». acquires 
TRM by alignment in a field H, the TRM is proportional to H for 
pH /kT,<1and saturates when pH/kT,> 1, k being Boltzmann’s constant 
and » is measured at the blocking temperature 7. The range of fields 
in which saturation becomes apparent depends upon the material but 
10 oersteds may be sufficient in many materials. The measurement of 
TRM therefore provides two tests for pseudo-single-domain effects, which 
must be invoked to explain any TRM which is much larger than indicated 
by (38) and is not proportional to the inducing field at low fields. 

The observations indicate that pseudo-single-domain effects are 
associated only with small grains. However, it is possible that permanent 
dipole moments occur on all grains more or less independently of their 
sizes, being relatively less important for the larger grains. 

During demagnetization of small grains the first process is the move- 
ment of domain walls to their equilibrium or minimum energy positions, 
that is the multidomain contribution to TRM is destroyed, leaving the 
permanent moments unaffected. The pseudo-single-domain TRM_ is 
only destroyed by randomizing the permanent moments. In the wedge 
shown in fig. 4 this would require the application of a field of strength 
sufficient to reverse both the domains of the grain, which is clearly greater 
than the field required to move the domain wall to the left or right. The 
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permanent moments of grains are therefore responsible for magnetically 
hard components of TRM which are of particular interest in palaeo- 
magnetism, although it must be expected that in most cases the multi- 
domain component is much larger. 

There is a further interesting consequence of the pseudo-single-domain 
behaviour which arises in demagnetization experiments. During com- 
plete demagnetization the permanent dipole moments are randomized, 
so that in a rock sample with n grains each of permanent moment p, 
one must expect to observe a randomly directed moment of average 
value about (n1/2) after demagnetization in the highest fields. Observa- 
tions of such random moments have evidently been made by Irving é¢ al. 
(1961). The effect is not consistent with perfect multidomain behaviour, 
because each perfect multidomain grain has zero moment in its 
demagnetized state. 

§ 7. Discussion 

There can be little doubt that magnetic grains in most igneous rocks 
are multidomains, but a quantitative description of their properties is 
difficult to give an account of the structural complexity of the grains. 
In this paper simple models are used to explain some of the important 
magnetic properties on the basis of conventional domain theory. Several 
predictions are made which will provide tests for the essential validity of 
the theory. An approximate magnitude is calculated for saturation 
anhysteretic magnetization; saturation remanence is shown to be deter- 
mined by coercive force (which is also related directly to the character- 
istic field of alternating field demagnetization). The induction of 
anisotropy of susceptibility by alternating fields is a consequence of 
multidomain structure and measurements of induced anisotropy of 
powdered magnetite, for which the factor (/,?/|$K, ++K,|) is known could 
provide an accurate test of the theory. Two further consequences of 
multidomain structure are considered qualitatively. 

Stacey (1960 b) and Stacey et al. (1960) found that the anisotropies of 
several types of igneous rocks which were examined were due to the 
alignment of elongations of magnetic grains and that the crystallographic 
axes of cubic magnetic minerals were randomly oriented and uncorrelated 
with the grain elongations. TRM in an anisotropic rock is stronger in 
the direction of the grain elongations than in the perpendicular directions 
because the demagnetizing factors of grains are smaller in the directions 
of their long axes. Also it follows from the considerations of domain 
structure in § 2 that grains in which the domains are parallel to the long 
axes are magnetically softer and are therefore more easily demagnetized 
than are grains of the same size and shape, but in which the domains are 
perpendicular to the long axes. Thus when a rock is cooled in a field at 
an angle @ to its axis of anisotropy it acquires TRM at an angle b <9, 
but by demagnetization ¢ is increased. A quantitative estimate of the 
change in ¢ is of considerable interest in palaeomagnetism, but it has 
not been found possible to make useful calculations of this effect. Firstly 
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it is not known whether the domain pattern in a cooling rock favours 
alignment with grain elongations; secondly calculations of changes in ¢ 
depend upon the differences in demagnetization curves parallel and 
perpendicular to the axis of anisotropy and these curves were shown in 
§ 3 to depend very strongly upon the precise grain shape. Only experi- 
ments will be able to show whether, during demagnetization, 4 reaches 
6 or overshoots it. However, if the multidomain theory is correct at all 
the change in direction of TRM occurring in an anisotropic rock during 
demagnetization must be such as to offset the deflection of TRM by the 
intrinsic anisotropy, at least qualitatively (Stacey 1960 b). 

Another interesting consequence of the multidomain structure of 
magnetic grains in rocks is the slight dependence of TRM upon cooling 
rate. Above the blocking temperature 7’, the domain walls in a grain 
are thermally excited past the potential barriers, which prevent move- 
ments below 7. Above 7, a grain thus has a magnetic moment in 
equilibrium with the external field (Stacey 1958). Quite apart from 
the fact that blocking temperature is actually a range of temperatures 
because the domain walls in a grain have different areas and activation 
energies, the value of 7’, is dependent upon the rate of cooling, because 
thermally activated domain wali motion is a rate process (Stacey 1960 a). 
The blocking temperatures 7, and 7’, corresponding to cooling rates 
T,'’ and T,' are approximately related by the equation 


1 { k (za 
a re ila os 7): ss eee oo) 
T 2 Ts E Ps 


where £ is potential barrier energy and k is Boltzmann’s constant. 
Slower cooling thus causes a lowering of blocking temperature, reduction 
of the ratio (o,/o,) and hence a smaller TRM. The effect may be appreci- 
able when the natural magnetic moments of very large igneous bodies 
are considered and brings an additional unknown into any attempt to 
guess the magnitude of the primeval earth’s field from the magnitudes of 
the natural moments of rocks. It may be noted that in single domain 
grains the effect is exactly opposite since lower blocking temperature 
results in higher moments » of the grains and therefore increases the 
factor (wH/kT,), producing stronger alignment. 


APPENDIX 

Mean and mean square direction cosines of the nearest easy directions 
of magnetization with respect to a fixed direction in an assembly of 
grains with randomly oriented crystallographic axes. 

Cubic grains with positive and negative magnetocrystalline anisotropy 
constants K, have easy directions of magnetization parallel to [100] 
and [111] crystallographic axes. In magnetite K, is negative and the 
easy direction is [111]. 
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The calculation of mean and mean square direction cosines of the 
angles between nearest easy directions of magnetization and a fixed 
direction in a random assembly of grains is equivalent to calculating the 
cosines for random directions with respect to the nearest easy direction 
in a fixed cube. We consider a cube of side 2d centred at the origin of a 
coordinate system, the sides being parallel to the x, y, z axes. A line 
emerging from the face (z=d) at (#, y, d) makes angles ¢ with the z axis 
([100] axis) and »% with the body diagonal (v=y=z) ([111] axis), 
where 


cos i ee 
in (d? + 02 4+ y?)1/2? 
cos f= TEE 


It is sufficient to integrate over one quarter-face since all other quarter- 
faces present simple repeats of the same geometry. The mean and 
mean square direction cosines are then: 


eee RS d 3 
cos? « a= Sf ne ab tes ==Q-701, 


For Kk, positive (+2? +y?) 
([100] easy) d 2 
cos “= a ib poe SOs een dy =0-818. 
(d? +a? + y?) 
2 
cos? a= vl ie ae ais pabtslls dx dy = 0-728, 
For K, negative (da? +2? +?) 
({111] easy) d 
cos a= — “{ [ s cons dx dy = 0-858;. 
(d24+ 22+ 42) y”) 
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ABSTRACT 


There is a close formal similarity between the equations of elasticity with 
anti-plane strain and the equations of electrodynamics in two dimensions, 
the screw dislocation being the counterpart of the line charge. This similarity 
has led to the suggestion that there is a specific ‘ Lorentzian ’ force on a dis- 
location which moves through a medium which is itself in motion. This is 
the analogue of the magnetic force on a moving charge. The equations of a 
linear continuum can be satisfied whether such a force exists or not. A non- 
linear model based on that of Peierls is constructed, and an approximate 
solution of a problem is presented which, if the approximations made are 
justified, demonstrates that this ‘ Lorentzian ’ force on a dislocation does not 
exist. 


§ 1. IntTRODUCTION 


FRANK (1949) first showed that there were strong formal analogies between 
the dynamics of a screw dislocation and the dynamics of a particle in the 
special theory of relativity, the velocity of light in the latter theory being 
replaced by the velocity of shear waves in the former. Eshelby (private 
communication 1950, 1953), by an extension of this formal analogy, 
suggested that a dislocation line which moved through a medium which was 
itself in motion might experience a force proportional to its velocity through 
the medium, analogous to the Lorentz force ona moving charge. A physical 
interpretation of this result was given by Nabarro (1951), and a simpler 
discussion. by Eshelby (1956). 

Interest in the nature of this force has been revived in discussions and 
correspondence with Gilman, Frank and Eshelby. Gilman has suggested 
that the existence of a Lorentzian force on a fast-moving dislocation may be 
a factor in explaining the observation of Johnston and Gilman (1960) that a 
dislocation loop may be expanded slowly in lithium fluoride without 
producing further dislocations, while an equal expansion in a shorter time 
leads to profuse multiplication of dislocations. On the other hand, Frank 
has given reasons for doubting the physical reality of the Lorentzian force, 
and has proposed a paradox which appears to show that the existence of the 
Lorentzian force on a dislocation leads to a conflict with the Principle of 
Special Relativity. The present approach is based on the idea that a 
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linear field theory cannot determine the equation of motion of its singu- 
larities. Accordingly, arguments based on a linear theory show that the 
motion induced by this force will be in accordance with the laws of 
mechanics, if it occurs, but leave open the possibility that it may not occur, 
or show that it is difficult to see what processes will induce this motion, but 
do not show that the existence of the motion is in conflict with the general 
laws of nature. The arguments in favour of the real existence of the force 
show that if the motion in the direction of the force occurs, external trac- 
tions have to be applied to maintain the prescribed velocities. When the 
motion occurs, these external agents do work, and we are accustomed to say 
that if a given motion of a dislocation allows external agents to do work, 
then a force acts on the dislocation in the direction of that motion. The 
arguments against the real existence of the force rest on the appearance 
of a term V, which is apparently the absolute velocity of the medium in the 
neighbourhood of the dislocation. Such a uniform velocity cannot affect 
the internal motions of a system. 

The nature of the force is such that if one part of a crystal is slipping 
over another by the spreading of loops of dislocations over a glide plane, 
and the crystal as a whole is moving in the same direction as one of the edge 
elements of each loop is moving, both screw elements of each loop tend to 
deviate towards the same side of the glide plane. Although the edge 
elements could leave the glide plane only by climb, a similar Lorentzian 
force acts on them, and, to simplify the sign conventions, we shall consider 
a two-dimensional case in which only edge dislocations are moving. 


§ 2. Toe Linear THEORY 


We consider the crystal of fig. 1. To avoid difficulties connected with 
the transport of matter in the climb of dislocations, we suppose the gliding 
of the upper part of the crystal over the lower with relative velocity v 
to be effected by the movement of equal numbers of positive and negative 
edge dislocations in opposite directions. We suppose that in a given frame 
of reference the upper and lower parts of the crystal have velocities in the 
direction of V+ 3v and V — }v respectively, while the glide plane is at a 
height ¢ above the middle of the crystal. The three quantities V, v and ¢ 
specify the motion of the crystal, but not its shape, which depends on the 
past values of v(t) and ¢(¢). The coordinates are thus rheonomic. 

We neglect any resistance to the motion of the dislocations. The initial 
rates of change of V, v and ¢ then depend on the constraints which are 
imposed on the motion. We shall consider constraints of three kinds: (i) 
free motion, (ii) linear motion with conservation of linear momentum and 
(iii) linear motion with v and V maintained constant. 

(i) We first consider the crystal to be moving freely in space. Then its 
kinetic energy, its linear momentum along Ox and its angular momentum 
around Oy must be conserved. These three conditions ensure that V,vand 
¢ remain constant, so that there is no Lorentzian force. | 
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When the crystal travels freely the conservation of angular momentum 
ensures that the relative velocity v continues to be directed along the x 
axis. If the conditions of constancy of kinetic energy and of linear 
momentum are replaced by other conditions, the condition of constancy of 
angular momentum can generally be satisfied only if the glide plane is 
allowed to rotate about the y axis. We therefore replace the condition of 
constancy of angular momentum by the condition of constant zero angular 
velocity of the glide plane. We do this by arranging that the upper and 
lower surfaces of the crystal are guided by frictionless planes at-z= + 4a. 
The forces exerted by the planes z= + 4a dono work and have no component 
in the Oz direction, so that kinetic energy and momentum along Ox are still 
conserved. They may introduce stresses into the crystal which exert 
forces on the dislocations, but these forces act in opposite directions on 
dislocations of opposite signs, whereas the Lorentzian force acts in the same 
direction on dislocations of opposite signs when they are moving in opposite 
directions. 


Fig. 1 


Vt lev - 


of 


V-‘av —- & 


Glide of a crystal. The mean velocity of the crystal is V, the relative velocity 
of the upper and lower parts v. The thickness of the crystal is a, and 
the glide plane lies a distance ¢ above the middle of the crystal. 


(ii) Having replaced the ‘natural’ condition of constant angular 
momentum by the condition of constant angular velocity which was 
implicit in the statement of the problem, we proceed to replace the natural ; 
condition of constant kinetic energy by the condition of constant relative 
velocity v which is stated in the problem. We may do this by means of @ 
mechanical device which is attached to the two parts of the crystal, but is 
not anchored to a fixed frame of reference. This device can alter the- 
kinetic energy of the crystal, but cannot alter its linear momentum. The 
linear momentum is proportional to ($a—¢)(V+4v)+(ga+ ¢)(V — 30). 
Thus, if the initial values of V, v and ¢ are Vo, v) and Co, the conditions of the 
experiment are now given by 

C=0, 
(Ja 0)(V + 30) + (a+ L(V — $e) = (4a Lo)( Vo + 300) + (b+ Ns 
(1 


VOCE Wal daseog Sue ba IA (2) 


These lead to 
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and a kinetic energy 7’ given by 
21 |pS = (ha —0)(V +40)2+ ($a-+.0(V — 40)? 
r = 

: 2aT oS =a? V 2 —2aloV oto + Corte + (407 — 07) ee) 
where p is the density of the crystal and S the area of the glide plane. 

In accordance with the usual definition of the force on an elastic 
singularity (such as a group of dislocations), we say that if the movement 
of the glide plane in a given direction allows the applied forces to do work, 
then a force acts on the glide plane in this direction. Inthe present problem 
the work done by the applied forces appears as kinetic energy. The only 
term in (3) which depends on ¢ is the last. This does not contain V5, so 
that the motion of the glide plane cannot be used to determine the absolute 
velocity V,. In fact the glide plane is in stable equilibrium in the middle of 
the crystal, and, if displaced, tends to move back to the middle, a satis- 
factorily ‘relativistic’ result. 

(iii) Let us finally replace the ‘natural’ condition of constancy of linear 
momentum by the condition V =constant which is stated in the problem. 
This condition cannot be expressed in terms of the internal state of the 


Fig. 2 


(a) (b) 


The upper half of a crystal glides with velocity U to the right with respect to 
the lower half, the motion being caused by a procession of positive edge 
dislocations moving from left to right. In (a) the bottom of the crystal 
is fixed and top moves to the right, while in (b) the top of the crystal is 
fixed while the bottom moves to the left. 


crystal, but only with respect to an external reference system. It is now 
possible to vary ¢ if linear momentum is no longer conserved. This 
requires a force in the Ox direction to act on the crystal as a whole, and this 
force must react on the fixed reference system. Once we recognize that 
the imposed conditions of the experiment require external forces to act, we 
see that the situations (a) and (b) in fig. 2 are not identical, any more than 
the conditions are identical if we require a rifle not to recoil when a bullet 
is fired at 7 x 104cm sec“ or if we require the rifle to recoil at 7 x 10¢em sect 
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while the bullet is held fixed. In fact under the imposed conditions V = V 
v =v, the kinetic energy is given by 


27 [pS = (4a—C)(Vot hup)®+ (a+ £)(Vo— 40»)? 

er neta aC eae une ean eae ee eee ge (4) 
If V, and v have the same sign, 7’ is increased by decreasing ¢, and con- 
versely, so that the system gains kinetic energy from its constraining 
mechanisms when the slip plane moves to increase the volume of material 
which is moving fastest with respect to the system of reference in which V 
is prescribed. By expressing vy in terms of the number and velocity of the 
dislocations moving in the glide plane it is easily shown that 07'/d¢ repre- 
sents a force on each dislocation equal to beV qv), where b is the Burgers 
vector. . 

As was to be expected from the fact that only linear mechanics is used, 
this calculation does not lead to a unique conclusion. We have shown that 
if the glide plane does change in the way predicted by the force which has 
just been calculated, the equations of mechanics will be satisfied. If the 
glide plane does not change, the equations of mechanics are still satisfied. 
Similarly we may suppose that the free motion of fig. 2 (a) is perturbed by a 
force applied to the upper block in that direction of its motion U. Then 
the equations of mechanics may be satisfied either by the glide plane moving 
downwards to increase the mass of the moving block, or by an acceleration 
of the block rigidly over the existing glide plane. Conversely it is strange 
that the internal behaviour of the crystal should depend on its mean velocity 
with respect to the reference system against which the applied force 
reacts, but the crystal is not an isolated system to which a principle of 
relativity can be applied. If the crystal and the massive reference system 
are treated as a single closed system, V is an internal coordinate, and neither 
the presence nor the absence of a force depending on it violates the ideas 
of relativity. 


0° 


§ 3. THE Non-Linear THEORY 


We attempt to resolve the question by constructing a non-linear model 
in which two screw dislocations of opposite sign lie in the plane y= vt, and 
the Peierls condition is applied between planes y=vt+4b. The disloca- 
tions are held in equilibrium against their natural attraction by an applied 
shear stress p,,. There is an analogous electromagnetic problem in which 
two line charges are held in equilibrium by an electric field #,. In both 
problems we may calculate the external field which is necessary to maintain 
equilibrium when the singularities are at rest, and then enquire how this 
stress is altered when the singularities move in the y direction with speed ». 
In both cases the Lorentz contraction increases the stress or electric field 
with which each singularity attracts the other by a factor (1 —v?/c?)~"?, and, 
if there were no other effects, the external stress would have to be increased 
by this amount in order to maintain equilibrium. In the electromagnetic 
case the moving charges each move in a magnetic field produced by the 
- motion of the other. This results in a repulsion between the charges which, 
to lowest order in v, is — v2/c? times their attraction when at rest. To this 
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order the total attraction between them is multiplied by a factor 1 —v?/2c?, 
and the external field HZ, required to maintain equilibrium is therefore 
reduced by this factor. The velocity-dependent terms in the interaction 
reverse the effect of the Lorentz contraction of the static field. 

The function w which exactly satisfies the Peierls condition between the 
planes y =vt + 4) and exactly satisfies the elastic equation of motion outside 
these planes is given by 


Qaw _ ees 4b(sec 6 — 1) sgn (y —vt)—vt 


b p(w — 1) 
ie y + 4b(sec 6 — 1) sgn (y —2t) —vt 
p(x +1) 
+ (y/b) sind +4(5—sind)sgn(y—vt), . . . . . (5) 
where 
p=(l—v/c?)¥2 Pecan em ad (i) 
and 
5=4sin— (b/pl). 2 oS ae cs 


The distance between the dislocations, 2/, is independent of v. The 
external stress at large distances which maintains equilibrium is 


Pyz= (u/27) sin 6. 

If the dislocations are far apart so that their cores do not overlap, 6 is small, 
and sin 6=6/2pl. Thus the external stress is proportional to p, which 
corresponds to the effect of the Lorentz contraction alone and the absence 
of velocity-dependent forces which would lead to a total equilibrating 
stress proportional top. The calculation is not wholly convincing, because 
it fails to consider the mechanisms by which the dislocations can move in 
the y direction and, more seriously, by which the state of being non- 
Hookean between the planes y=vi + $b can propagate. 


§ 4. CoNCLUSION 

The existence or non-existence of a Lorentzian force on a dislocation 
which moves in a medium which is itself in motion can be proved only by 
constructing anon-linear model of a dislocation. We have obtained an 
approximate solution for such a model, and, in this solution, the Lorentzian 
force is absent. It is possible, but hardly likely, that it would appear in an 
exact solution. 
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ABSTRACT 


An elastic calculation shows that there is an interaction between two 
jogs or kinks on an infinite straight dislocation line. For jogs sufficiently 
far apart, the interaction potential varies as the reciprocal of the distance 
between them, and as the product of their lengths. Formulae are also 
obtained for the self energies of jogs. 


RECENT papers on plasticity have emphasized the role of jogs en dis- 
location lines in limiting the velocity a dislocation can achieve under a 
given applied stress (Mott 1960, Hirsch 1960, Gilman and Johnston 
1960). Such theories are based on calculations of the rate at which 
jogs are created by forest cutting, and the rate at which they are 
annihilated by glide or climb. Two important quantities are required 
for these calculations: the self-energy of a jog, Hg, and the interaction 
energy between two jogs, H;. We shall calculate these quantities using 
the model of a dislocation familiar in the classical isotropic theory of 
elasticity. We shall assume that both jogs are perpendicular to an 
infinite straight dislocation line. 

The method is illustrated in fig. 1. When the loop 2 is added to a 
straight dislocation 1 with the same Burgers vector b, the resulting 
configuration is that of a dislocation line containing two ‘ antiparallel ’ 
jogs. The total elastic energy H of the dislocation with two such jogs. 
is given by 

H=#H,+H,+£,, Ne ent aaa ae Ly 
where H, and #, are the self-energies of the dislocation line and loop— 
respectively, and H,, is the interaction energy between them. 

The ‘ double jog energy ’, Hjj;, can be defined as 

Byala Hot las. oe (2) 

It can be divided into two parts. That part which does not depend on 

the distance r between jogs, is twice the self-energy of one jog Hy. The 

part dependent on r is precisely the interaction energy between the two 
jogs, H;. Thus 

Hjj=2Hg + Hy(r). Se nis rt Be CO) 


+ On leave from the Institute of Physics, Prague. 
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The Burgers vector b can be resolved into three components (fig. 1): 
b?, perpendicular to the plane of the jogs; 6® in the direction of the 
straight dislocation (such that it is screw); and b" in the plane of the 
jogs but perpendicular to the straight dislocation (such that it is edge). 
We shall treat these cases separately. H, and H, for a general Burgers 
vector b are given as sums 


Ey=EP + HE + HS; Hy= HP + H+ Hs 
for 
b=ib8+jpe+kb®, . |. 4 ad) 


Fig. 1 


The self-energy of a rectangular loop can be calculated by the method 
of Yoffe (1960). For a prismatic loop, with its Burgers vector b? normal 
to the loop plane, we have: 


(OF 2 2a ee (a? +r?) +a 
Gi galas ole ih loo bank et Sa ENS Ae eae egies 
: {4 2e € melee € 31 \/(a?+r?)—a 
r, a/(a?+r7)+r 
— Fog Sarpy Teveterty— mata t 


where yp is the shear modulus, v is Poisson’s ratio, and « is the core cut-off 
radius, «xb. The length of the loop is r, and its width is a. 
For a glide loop, with its Burgers vector parallel to the side a, we find 


p_ Bb")? 2a_ 1, V@ir)ta  Y(a?+r?) a+r 
E, arr a} log log ee 
e 


2 4/ (a? +r?)—a@ a a 
rete ot — Slo sete V(a+r?) at 
Ged Rewrite secrete || 


The ee E,§ of a similar loop, with its Burgers vector b8 parallel 
to r, can be found by interchanging r and a in the above formula. 

The interaction energies are easily found by computing the work 
required to make the loop in the stress-field of the infinite dislocation. 
In the prismatic case we find for instance: 

p(OP)* 


Bh iP te eS : 
ae Selice <td ee 
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The interaction £,," is given by eqn. (7) when the superscript P is 
replaced by E. H,,§ is given by (7), replacing P by S and omitting the 
factor (1—v)-1. The double jog energy, Hjj®, is given by eqn. (2), 
using eqns. (5) and (7), and similarly for #j;® and Hj;%. 

Dividing Hj; up as in eqn. (3), one finds for the self-energies 


P)2 9 
e* = AS | tog = 2 | . 


W4r(l=y € 
bE)2 2 2— 
ji ee log = — if : Mt ats) 
47 € l-—yp 


s_ HOS) P, 2a, 
ONS ean ae 

Note that these elastic energies all become zero for a jog length ax 4e. 
Thus, for small jogs, or jogs of from one to four or five Burgers vectors 
long, the elastic approximation breaks down, and most of the jog energy 
must be considered to come from atomic arrangements in the core of the 
dislocation. Our formulae are therefore useful only for long jogs. 

When a/r<1 so that powers of a/r greater than one can be neglected, 
we find for the interaction energies: 


) 
Ti 


87 l—yp 

pea — HON ee, Bee sy 
87 l-—-yvy | 7 

psa — Wey? eae | 

87 Ll—vjr° | 


The total force on a single jog can be found from differentiating the 
above expressions with respect to 7. In all cases, the jogs attract one 
another with a force which falls off with distance as r~?._ The force also 
increases as the square of the jog length. Two antiparallel jogs behave 
thus as equal and opposite electrical charges. 

The method used above cannot be applied directly to calculate the 
interaction energy between antiparallel jogs of different length, or between 
parallel jogs. For the general case, when both jogs are perpendicular 
to the dislocation line, but have different lengths and lie on different 
planes making the angle « with each other (fig. 2), the dislocation with 
two jogs can be composed from a straight dislocation 0, and two semi- 
infinite loops 1 and 2. The total energy E of the jogged dislocation is. 
the sum of the self- and interaction energies 


H=H,+H,+#,+ By + Lo.+ Fi. 


Only the interaction energy H,, between two semi-infinite loops 1 
and 2 depends on the distance 7 between the two jogs. Thus 


Ey= (7). 


To calculate this interaction energy, we must know the stress field of 
the loop 1. This can readily be found using the expressions for stresses 
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of an infinitesimal loop (Kroupa 1962) at the point (0, 0, z) and integrat- 
ing over the area of the loop. Since we are only interested in the case 
where r>4a,, r>d@ , we can treat the widths of the loops as infinitesimal, 
effectively replacing all double integrals by single integrals. In this 
way, the interaction energy is found to be, for general b 


=B,.= — — — J [(bP)2 + (68)2(1 + v) + (0F)2(1 — 2) ] cos « 
By=B yy — f(b) + OS) +0) + OI 
— 200° sin a bate gt ode We eeeann 10) 
Fig. 2 


For «=0 and a,=a,, eqn. (10) reduces to the special case of anti- 
parallel jogs of length a (eqns. (4) and (9))._ When the jogs are parallel, 
%=7, coSa=—I1, sin x«=0; the interaction energy is positive, and the 
jogs repel one another. Thus the analogy with electric charges is very 
close in that unlike (antiparallel) jogs attract, and like (parallel) jogs 
repel, provided that the jogs are in the same plane, and far apart. 
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ABSTRACT 


The cosmic ray intensity was measured in a mine using a combined 
scintillation counter and Geiger counter telescope. It was shown that this 
apparatus could discriminate effectively against local y-radiation. The 
intensities at 1660 m.w.e. and 3280 m.w.e. were 3x 10-5 and 10-6 of that 
at the surface ; no cosmic ray events were observed at 5050 m.w.e. The 
results suggest a steepening of the energy spectrum in the region of 10! ev, 
The most probable energy loss of y-mesons which had a mean energy of about 
3 x 104+ ev was found to be not significantly higher than that of minimum 
ionizing particles. 


§ 1. INTRODUCTION 


PROGRESS in the measurement of the intensity of cosmic rays deep 
underground has been remarkably slow, although it has always been 
recognized that such information would be of considerable interest. 
Since the review article of Barrett et al. (1952), which showed that data 
on the intensity were then established down to a depth close to 2000 m.w.e., 
no further satisfactory results have been published. There appear to be 
two major difficulties which have to be surmounted. Firstly, there are 
only a small number of mines which are of sufficient depth and are sunk 
in rock of low radioactive content. (The alternative technique of working 
in the deep oceans presents considerable technical problems which have 
not yet been solved in an economically feasible manner.) Secondly, it is 
necessary to measure a very feeble flux of mesons in the presence of 
y-radiation from the surrounding rock. Even in a mine of relatively low 
radioactive content the ratio of the flux of mesons to that of y-rays is less 
than 10~® at a depth equivalent to a few thousand metres of water; such 
aratio would be regarded as troublesome in an accelerator or pile experiment 
and is even more so when working away from a laboratory. An attempt 
made using a Geiger counter telescope in 1956 and reported recently 
(Barton and Michaelis 1961) failed for this reason. A new apparatus was 
therefore designed which would make the identification of cosmic ray 
events more definite, and readings were obtained with this apparatus in 
the summer of 1960. So far, only a small number of readings have been 
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obtained but these are being reported now as it may be some time before 
a more satisfactory set of data is available. 


§ 2. APPARATUS 


This was a modification of that used earlier (Barton 1956) with the 
middle tray of Geiger counters replaced by a large scintillation counter, 
fig.1. The block of plastic phosphor, of size 53-0 cm x 15-0 cm x 3-0 cm, 
was viewed through a short light guide (at one end) by two type 51 AVP 
photomultipliers. The optical arrangement was dictated by the need to 


Fig. 1 


PHOSPHOR 
IScm. 
Section through apparatus. 


fit the: scintillation counter into the same space as the Geiger counters 
which had been replaced. As a result the uniformity of response was 
slightly poorer than that obtained with the rather similar counter for which 
full information was published recently (Barnaby and Barton 1960). 
A block diagram of the circuits used, which were much more complicated 
than those used before, is shown in fig. 2. The fast comcidence between 
the two photomultipliers was included as an additional precaution against 
spurious events. 

The magnitude of the pulses from the last dynodes of the two photo- 
multipliers, connected in parallel, was measured by a ‘ringing coil’ type 
pulse height analyser. The damped oscillation of a tuned circuit, excited 
by the incoming pulse, was amplified and then passed to a simple integral 
discriminator which triggered once each cycle until the oscillation had 
decayed below a certain fixed level. In this manner a chain of pulses was 
produced whose number was proportional to the logarithm of the height 
of the original pulse. This number was converted into binary form and, 
together with the hodoscope data from the eight Geiger counters, was 
recorded on fourteen tracks on magnetic tape moving at 2-5 em hour. 
With this type of analyser the width of each channel is proportional to 
the pulse height accepted by it. In this particular analyser the incremental 
channel width, which could be controlled by the resistance across the tuned 
circuit, was set at 10% and was found experimentally to be uniform from 
the second to at least the fortieth channel. 
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Block diagram of circuits. 


The e.h.t. supply for the counters and photomultipliers was provided 
by a voltage converter circuit derived from that given by Hooker and 
Light (1956). A corona stabilizer was included in a feed-back loop to 
control the oscillator and resulted in a stability of 0-1% for the changes 
in battery voltage and ambient temperature which were encountered 
under operating conditions. 

The apparatus contained a total of 80 transistors which consumed in 
all rather less than 100mw. All the circuits were operated from 
‘Mallory’ type dry cells which provided a satisfactorily uniform voltage. 
The most critical stage, the discriminator in the pulse height analyser, 
was further stabilized by a Zener diode. Both batteries and tape recorder 
could operate for up to five weeks at a time without attention. Laboratory 
experiments showed that the counting rate of the apparatus was indepen- 
dent of temperature within the range 10°c to 40°c. There was a slight 
dependence of pulse height on temperature, and the results have been 
corrected for this effect whenever the temperature underground was 
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different from that at the surface: the correction necessary was a 
difference of only one or two channels. 


§ 3. OBSERVATIONS 


Most of the data was obtained during a visit to the Lake Shore Gold 
Mine, Northern Ontario, in the summer of 1960. This mine is situated 
at latitude 48° 9’ N, longitude 80°3’ W; the surface of the mine is 320 m 
above sea level and does not vary by more than 30 m for several miles 
in any direction. 

The rock of the Canadian Shield in the region of this mine, except for 
the thin ore-bearing fissures, is almost entirely syenite-porphyry of 
density 2°73, mean Z=11-7 and mean Z/A=0-50 (Thomson 1948); 
one foot of rock is thus ‘equivalent’ to 0-83 m of water. 

In the first set of experiments the apparatus was operated at the surface 
and at depths of 2000 ft and 6075 ft. Histograms of pulse height data 
from all the events recorded are presented in fig. 3, where the data have 


Fig. 3 
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Pulse height histograms for first series of experiments. 


been grouped into pairs of channels. The data obtained at sea level and 
at 2000 ft show the characteristic ‘Landau’ shape, though the asymmetry 
is less striking than if the pulse heights had been plotted instead of their 
logarithms. On the other hand the 15 events registered at 6075 ft appear 
to be of an entirely different character since the mean pulse height is only 
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of the order of one-third of that recorded elsewhere. Further information 
can be obtained by examining the Geiger counter hodoscope records. 
If only one counter has been discharged in each tray the inclination to 
the vertical of the particle trajectory can be determined. Those events 
in which more than one counter is discharged in either or both trays are 
referred to as multiple events, although in a proportion of cases a single 
particle may be responsible. The resulting analysis is given in table 1. 

Allowing for the appropriate statistical and geometric weighting factors 
it is found that the results at the surface are consistent with the expected 
approximately cos? @ distribution. 

The results for 2000 ft indicate a more collimated distribution and a 
higher proportion of multiples. Both these effects are to be expected and 
should be even more noticeable at 6075 ft. Instead the events at the 
latter depth appear to be more widely spread in angle and include no 
example of a multiple event. This strengthens the suggestion made above 
that the events here are of a different nature. 


Table 1. Hodoscope results for first series of experiments 


Mean projected inclination to vertical | Multiple Total 

0° 23° 41° 53” events : 
Surface 88 94 28 3 17 230 
2000 ft 33 30 + 1 15 83 
6075 ft 2 ad 5 1 0 15 


Mean projected inclination Maleole 
to vertical Pp 
OF 237 41> 253" 


Total 
events 


Surface 76: 92 
3950 ft logh <13 93 131 
3950 ft log h>=18 1 3 


In order to investigate the nature of these other events more closely 
the apparatus was adjusted so that all the scintillation counter pulses 
were larger but the minimum signal required to trigger the master 
coincidence circuit was smaller. Histograms then obtained at the 
surface and at a depth of 3950 ft are shown in fig. 4. 

The distribution for the events recorded at the surface is not significantly 
different from that obtained before, but is shifted to a higher mean value. 
The two events with a very small pulse height are probably due to particles 
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passing close to the end of the apparatus as the Geiger counters were 
several centimetres longer than the block of scintillation phosphor. The 
events at 3950 ft clearly divide into two groups of which the lower one 
can be compared to that registered at 6075 ft but, owing to the different 
discriminating requirements, is very much more striking. The hodoscope 
data are presented separately for each group in table 2. 


Fig. 4 
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Pulse height histograms for second series of experiments. 


It can be shown that the smaller pulse height events at 3950 ft are 
compatible with an isotropic distribution. These events are clearly 
similar to those recorded in the earlier experiment (Barton and Michaelis 
1961) using a three-fold Geiger counter telescope. The evidence presented 
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there leads to the conclusion that such events were due to y-radiation 
from the surrounding rock. Only the other nine events are to be inter- 
preted as due to cosmic rays, for which the high proportion of multiple 
events is not unexpected. By making a similar analysis of the results at 
6075 ft it follows that all the events at that depth can be attributed to the 
effects of y-rays, and that there is no evidence for any cosmic rays there. 
This means that it is only possible to deduce an upper limit for the cosmic 
ray intensity. 

There is one other point of interest about table 2. Among the y-ray 
events recorded at 3950 ft there were two events in which two counters 
were discharged in one tray. This means that, including the scintillation 
counter, a total of four counters was discharged. In the earlier investigation 
of Barton and Michaelis (1961) no evidence was found for four-fold coinci- 
dences due to y-rays but it would not have been possible to exclude a very 
low rate of such events. 

The results for the cosmic ray events in the two sets of experiments are 
combined in table 3. The mean channel number (i.e. the mean of the 
logarithm of the pulse height) has been adjusted for the second experiment 
so that it can be compared directly with the first series of experiments. 
Also the number of multiple events has been reduced to take account of 
the proportion which can be attributed to a single particle passing through 
two counters in a tray. 


Table 3. Summarized results 


Corrected | Proportion 


: Mean number of 
Time |Counts| Normalized rate} channel of multiple 


aes multiples| events 


60 sec 448 1 ; 2 22 5% 
4 days 83 | (8-2+0-4) x 10->| 16- : 16° 

15 days 9 (O23) sc gk: 45% 

21 days 0 107 — 


§ 4. ANALYSIS 


The mean channel number, used in table 3, is convenient as an easily 
observed quantity but is not of direct physical significance. A more 
detailed analysis of the histograms of pulse heights, in the manner 
suggested by Barnaby (1961), does give rather similar results. The 
increase in the most probable pulse height between the surface and 2000 ft 
was found to be 54+ 3% and between the surface and 4000 ft 1145%. 

There is a possibility of a small additional systematic error so that any 
rise in the energy loss might be between 0 and 10% at 2000 ft and 
between 0 and 20% at 4000 ft. The minimum energies of the mesons 
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which penetrate to these depths are at least 5 x 104 ev and 1-4 x 10” ev 
respectively. The mean energy at the underground location of those 
mesons which succeed in penetrating can only be estimated very crudely 
since the values depend on the shape of the energy spectrum; estimates 
of 3 x 10" ev and 6 x 104 ev are in accord with the spectrum deduced here. 

Recently Barnaby (1961) has shown that between 10° ev and 101° ev 
there can be no increase in the most probable energy loss greater than 
1%, in agreement with the predictions of Sternheimer’s (1956, 1959) 
theory. The present data suggest that any increase at higher energies 
must be very gradual. 

The result is not at variance with the increase in the proportion of 
multiple events underground. These would be expected to contribute 
mostly to the high energy tail of the Landau distribution and have relatively 
little effect on the most probable energy loss. There are not sufficient 
of these events to enable any detailed deductions to be made but the 


Fig. 5 
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increased number can readily be explained by the increasing importance 
of electromagnetic interactions by ,.-mesons in this energy range. (The 
contribution from multiple penetrating particles, either surviving from 
air showers or else produced locally, is extremely small.) Figure 5, which 
includes additional results obtained at shallower depths, shows that the 
cross section for these interactions is consistent with a dependence on the 
logarithm of the meson energy, as predicted theoretically (see, e.g., 
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Fowler and Wolfendale 1958). Not all the soft component at sea level 
is in equilibrium with the ,.-meson beam so that the results obtained at 
the surface are not directly comparable and have been omitted from fig. 5. 

The total intensity data of table 3 need two corrections in order to 
construct a range spectrum of u-mesons. Firstly, mesons which miss the 


Fig. 6 
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apparatus may nevertheless be recorded through secondary particles 
produced in the rock close to the apparatus. The importance of this effect 
will be proportional to the number of multiple events but it is hard to 
calculate the necessary numerical factor; it must be between 0 and 1 
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so that to adopt a value of } will be reasonable and sufficient for the 
immediate purpose. Secondly, the vertical intensity is not quite pro- 
portional to that measured with a wide angle telescope because of the 
increasing collimation of the »-meson beam with depth underground. 
The angular distribution results of Bollinger (1951) are sufficient to enable 
an adequate estimate of this factor to be made. 

The corrected results are plotted as a function of “water equivalent’ 
depth in fig. 6 together with other available data for the intensity at 
great depths. The upper limit for the intensity at 5050 m.w.e. has been 
chosen so that the probability of the true rate being greater is 16% 
(Regener 1951). The agreement is seen to be good except that the 
measurement of Miyazaki (1949) at 3000 m.w.e. is not consistent with the 
new results. There are several reasons for not being surprised at this 
discrepancy. Miyazaki’s results were obtained with a simple telescope of 
Geiger counters which, as shown by Miesowicz and Massalski (1950) and 
by Barton and Michaelis (1961), could be sensitive to local radioactivity. 
The reported 75% decrease in intensity on inserting 100 cm of lead in 
the telescope suggests that there was some such effect. Also the data 
are only reported very briefly and there appears to be a discrepancy 
between the results for two separate experiments. 

The intensities measured at various inclinations to the vertical at one 
depth by Bollinger (1951) can be plotted on an intensity—depth curve, as 
is done in his thesis. However, estimates made in this way cannot be 
compared easily with results obtained directly by measurements at different 
depths because the effects of meson interaction and decay are different 
in the two cases. Allowing for this factor, and the uncertainties in it, 
there seems to be adequate agreement between the two sets of results. 
Neither set has sufficient accuracy for it to be possible to make any useful 
deductions about the problem of the proportion of K to 7-mesons among 
the parents of the .-mesons. 

To convert the intensity—depth curve to an integral energy spectrum 
it is necessary to assume the form of the range—-energy relationship for 
high energy ,.-mesons, which in turn is obtained from their rate of energy 
loss. The results for depths less than 1000m.w.e. have been analysed 
recently by Ashton (1961) so will not be considered here. 

It has been customary (Barrett et al. 1952, Ashton 1961) to assume 
that the rate of energy loss is something like 


dH 
=- =2:05+0-175 logy, 


= 708 + 3-0 x 10-§ H Mev g-! cm?, 


103 + 11-3 x 

Unfortunately there is considerable uncertainty about the correct value 
for the numerical coefficient in the third term. Although the total cross 
sections for the various processes (bremsstrahlung, pair production and 
nuclear) are fairly well established experimentally (Avan 1956, 1958), 
neither theory nor experiment provides reliable estimates of the partial 
cross sections when the energy transfer is large. These rarer high energy 
interactions do not contribute appreciably to the total cross section 
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but may make a large difference to the mean rate of energy loss. This 
applies especially to the theoretical estimate by Kessler and Kessler (1957) 
of the nuclear interaction cross section. It is therefore possible that the 
third term may need to be anything up to twice as large as given above. 

A further uncertainty arises because the energy loss can no longer be 
considered as a continuous process once the additional energy loss terms, 
particularly the bremsstrahlung, become important. The problem is then 
analogous to that of B-ray absorption and it is necessary to consider the 
effects of straggling. This has been treated theoretically for .-mesons by 
Mando and Sona (1953) and by Rozenthal and Strel’tsov (1958); their 
assumptions and numerical results are rather different but they agree in 
predicting that the intensity at a given depth will be rather greater than 
that calculated without allowing for straggling. 

If the energy loss relationship is given correctly by the expression 
adopted by Barrett et al., any high energy transfers are sufficiently rare 
that it should be reasonable to ignore the correction for straggling. If, 
however, the second term has been underestimated then the number of 
high energy transfers is correspondingly larger and the effect of straggling 
may no longer be negligible. Assuming that at least the sign of the effect 
is given correctly by the theories, then it will tend to mitigate the larger 
rate of energy loss. 

The data have been plotted as an energy spectrum in fig. 7 for both 
low and high values of the energy loss expression. The true spectrum 
should lie between the two curves and, as explained above, is unlikely to 
be very close to that drawn for the higher value. It is therefore probable 
that the exponent of the energy spectrum is increasing at energies above 
5x10“ev. This result is not in disagreement with the sea-level measure- 
ments of the momentum spectrum (Ashton et al. 1960), where the exponent 
of the spectrum was found to be approximately constant, since the 
statistical accuracy of their results is not sufficient to rule out an appreciable 
change at energies greater than 5 x 10“ ev. An indication of the steepening 
_ of the meson energy spectrum at very high energies has also been reported 
by Duthie et al. (1961) from the analysis of the spectrum of 7° mesons 
which produce electron—photen cascades in nuclear plates exposed near 
the top of the atmosphere. The results in the paper suggest an absolute 
magnitude of the spectrum exponent which is larger than that which would - 
be deduced from fig. 7, but subsequent work by the same group (Perkins, 
private communication) leads to a meson spectrum which is in. agreement 
with that deduced from the underground observations. 

The lack of any cosmic ray events during the 21 days of observation 
at the greatest depth is of interest in setting an upper limit to the cross 
section for the interaction of neutrinos. As far as neutrinos emitted by the . 
sun are concerned the limit is much higher than the recently observed 
value, but for neutrinos produced by cosmic ray decay processes there is 
the possibility that, due to a resonance in the cross section at high energy, 
the cross section may be very much higher. An estimate of the possible 
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rate which would be expected on account of anti-neutrinos has been 
given by Glashow (1960). It has already been pointed out (Barton 1960) 
that the present observations put the upper limit within the region predicted 


Pah e: 
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Integral spectra of «-mesons. b is the coefficient of the third term in energy 
loss expression. 


by Glashow. Unfortunately neither the statistical weight of the 
observations nor the uncertainties in the theoretical treatment allow any 
definite conclusion to be drawn. 
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§ 5. CONCLUSIONS 


1. The intensity of the cosmic radiation at a depth of 3280 m.w.e. is 
about one-millionth of that at the surface and at 5050 m.w.e. it is 
unlikely to be more than one ten-millionth. 


2. The most probable rate of energy loss of .-mesons is not more than 


10% greater at 3 x 10!4ev and 20% at 6 x 10! ev than at 10%ev, and may 
be the same. 


3. The proportion of secondary particles increases slowly with the 
energy of the mesons. 


4. The energy spectrum of ,.-mesons appears to become steeper a 
energies greater than 5 x 10" ev. 


} 
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ABSTRACT 


By employing a simple model of point charges in a uniform electron gas, 
it is shown how interaction energies between imperfections (for example, 
impurity—vacaney or vacaney—vacancy interactions) may be evaluated 
within the framework of low-order perturbation theory. Close contact may 
be established throughout with the earlier investigation of Alfred and March 
(1957), m which, however, a restrictive assumption that the perturbing 
potential was slowly varying had to be made. Our work now removes this 
assumption completely and provides a firm basis for the so-called ‘ electro- 
static model’ of the interaction between imperfections. In view of the 
long-range oscillations in the electrostatic potential around a defect, the 
interaction energy also has an oscillatory behaviour as a function of distance. 

The relevance of our work to the Lazarus theory of impurity diffusion in 
metals is discussed. In this case, the numerical results are in fair agreement 
with those of the earlier treatment, and lend further support to the Lazarus 
theory. Also the investigation of Seeger and Bross (1956) on the di-vacancy 
may be critically examined by comparison with the present findings. 
It then appears that the terms in the total interaction energy which these 
earlier workers have omitted may be large. Indeed, in our calculations they 
are sufficient to change the sign of the interaction. It is finally pointed out 
that the present theory may have relevance in the interpretation of the 
results of Biondi and Rayne (1959, Rayne 1961) on optical absorption in 
Cu-Zn and Cu-Ge alloys. 


§ 1. IyTRODUCTION 


We develop here a theory of the interaction between defects in metals, 
based on careful study of an admittedly over-simplified model. Thus we 
adopt the uniform gas approximation for the matrix metal, and assume 
that the point defects may be simulated by embedding appropriate point 
charges in the gas. Departures from such a model, particularly the use 
of Bloch waves rather than plane waves, and hence a non-uniform periodic 
density, for the conduction electrons in the unperturbed metal certainly 
require full investigation, but at present appear to lead to major 
complications. ite 

We shall refer briefly here to a few of the numerous applications to 
which our interaction energy calculations are highly relevant. Lazarus 
(1954) developed a theory of impurity diffusion in metals based on a 


1286 G. K. Corless and N. H. March on 


calculation of the interaction energy between an impurity and a vacancy 
(for a review of diffusion in metals, see Lazarus 1960). While Lazarus’s 
earlier treatment was open to certain objections, particularly with regard 
to the neglect of correlation effects in diffusion (see, for example, Bardeen 
and Herring 1952, Manning 1958) and to the calculation of the energy of 
movement, the recent reformulation of the theory by Le Claire (1961) 
appears to go some way towards removing these unsatisfactory features. 
The other application of our work which we consider in some detail is 
to discuss the stability of the di-vacancy and to make an analysis of the 
earlier calculation of Seeger and Bross (1956). 

The present work has been made possible by recent progress in the 
general theory of defects in metals (March and Murray 1960, 1961). While, 
in this earlier work, explicit results were obtained only for a. single 
imperfection in an infinite metal, a full perturbation development of the 
Dirac density matrix was presented, and this is applicable to a general 
perturbing potential. We show here how considerable progress can be 
made in the two-centre caset (for example, impurity—vacancy or 
vacancy—vacancy complexes), using a first-order approximation. Our 
results reduce to the first-order treatment of Alfred and March (1957) 
when the assumption of a slowly varying potential is made. While these 
workers obtained the same final result as Lazarus for the interaction energy, 
it was clear that his original argument was valid only because of a major 
cancellation of the kinetic energy terms with part of the potential energy. 
We shall show by a quite general first-order argument that this cancellation 
occurs, independently of any assumption concerning the spatial variation 
of the potential. 


§ 2. DENSITY Matrices anp Kinetic ENercy DENSITY 
Following March and Murray (1961), we write for the Dirac density 
matrix p(r,ro,k), where k?=2H, with H denoting the Fermi energy: 

p(t, Fo, &) = py + px + pp + 0(V®) 2 te meee ed 
where V is the perturbing potential energy. We then have the equations 
defining po, p, and p, as 
k2j,(k[r— rol) 
meal 

oa fan bern 
—Fo||r—ry| 
p= Ty ei] dr, dr, V(ry)V( he klry—rol+lro—ry|+k|r— rel) 
us [ry —rollro—ry||r— rol 
We note first that the kinetic energy 7’ is given by 


ype Het | dtlV 2p] eee: 5 te ee tan) 


} The generalization to M centres is readily achieved, as shown in Appendix 1. 


= (2.2) 
(2.3) 


(2.4) 
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Using (2.2)-(2.4) we then find by straightforward calculation : 
KA jy" (kr — rol) 
7 me rol 
2h? V(r)j4(k|r — rol) 
7 ene rol 

ibe V(r) jy" (Ale — ry [+ &[r tony 
eee a wJ1 1 1— Fol) 

ra | 2 [r—ry [ri — rol ae 


team ay | | desde, Meo edi ele elt eens Bead 


lr, —Frollro— r,||r—ro| 


Mas fer V (ry) V(r)47j,(k[ry—ro|+h|r—r,]) . 
ek LS oe ) ed EN 


Ant lr; —ry|[r—r,| 


Evaluating (2.6), (2.7) and (2.8) at ry>=r, and inserting the results in (2.5) 


we find 
arise == | Jas sa | anv (r;) 


(2k|r.—ry|) , 4):(2k|r.—r,]) 
+ 5y | | dridraV(r) V(r Jr (2k\re—ril) , ii(2k\re—m) ] og 
=a) Gra (ri) (eee “2lrs—r] * 2kry—F,)2 on 


At this stage we make use of the relation that to 0(V?) the displaced 


charge p— py is given by 
ae oa | anv (r,) Jt ae 
r,| 


V.Po= ; pee eo (2.6) 


(2.8) 


es Ifenan? ieyV (ro)ji(k|ry — r|+4|r.—ry|+ Ar — rol) (2.10) 
lr, —r||ro—ry||r—ro| 


and integrating (2.10) with respect to r we find 


{ dr(p—p.)=- | dr,V(r,) 


+ =i | [ arvana (rs) V(re) ut Vrs) V(ra)jol2elre—"al) 4 o¢yay, (2.11) 


Thus, fe eer we may rewrite (2. in the form 
V (ry) V (re) ja V (ry) V(2))1(2k] 2 — Fs) ra) 2.1 
roa) ng 3 dre pont =l{ ded ~) [e@kr=r[P 7 Gag 


§ 3. InTERACTION ENERGIES BETWEEN DEFECTS 


We turn now to deal with the energies of interaction between two 
defects, carrying charges Z, and Z,. We obtain the interaction energy by 
considering the changes in total energy as we pring the imperfections from 
infinity to the near-neighbour distance a. 
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The basic equation of the perturbation theory of March and Murray 
(1960) for the perturbing ere V, to first-order in V, is 


V2V(r =f arr ee ana eS) 


where k is the magnitude of the wave vector at the Fermi surface. This 
equation reduces to the Mott form 
: VV HRV, g=4h/a es 2ew: 
when V is slowly varying in space. It will be seen that the superposition 
property of (3.2), exploited by Alfred and March (1957) in their two-centre 
work, is also valid for (3.1), and thus, within the first-order approximation, 
no new calculation of the density and potential for the defect complex 
under consideration is required. Thus, if charges Z, and Z, are situated 
at a, and a, respectively, and if V,, and V,, represent the one-centre 
potentials as given by March and Murray (1961), then the two-centre 
potential may be written as 
V(r) =V 7,(r—a,) + Vzo(r — a9). PO Sete em ONE 

We now proceed to obtain the desired kinetic energy change from (2.12), 
and in doing so we note that if we measure changes from the unperturbed 
gas state and use charge normalization requirements, only the third term 
on the right-hand side of (2.12) needs to be considered. We have then, for 
the complex Z,Z, at separation a=|a,—a,|, the contribution from this 
an as 


S| | aradeatValer— ay) + Vealey— a)I(V alte as) + Vato ae)] 
x jx(2k|ry—ry|)?. (3.4) 
The one-centre terms are now seen to correspond to the kinetic energy 


changes at infinite separation of Z, and Z, and thus the kinetic energy 
change AT’ is given by 


AT = ks {| dr, dr, Vzalty — 81) V zo(tg— a5)j,(2k|r.— r,|) 
7 2G Tie ae 


(2k[r,—r,|)? 
ke Vaal — a2) Vza(he— a) jr(2k|re—ry|) | 
Uae aeud Z2\0y 2)V zi\lo wJ1 ie 
+S ff i ad em erseey) Agra jue erties NE 


We turn next to the evaluation of the potential energy change AU in 
this approximation. As we are working in a Hartree framework, we have 
that the potential energy of the gas with charges Z, and Z, at a, and a, 
is 


| rote rate +4 [ar dr, ae (rs )e(rs) + Meo (3.6) 
Tro—ri| |a,—a,| ; 
where 
Z Z 
V Ss 1 oS es a 
o(F) Reece = V,(r—a,) + V,(r—a,) 


and V (ry) =Vo(ry) + ae en 
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Clearly (3.6) may now be expressed as 
2 [ arpcryr( r)+4 [ aror, BER. (3.7) 
* fay —a,| 
Using (2. 1) for p(r), adopting the superposition potential energy (3.3), 


and omitting unperturbed and one-centre terms we immediately find to 
second order in V, 


AU=}3 4 drpyz1(r—a,) V go(h — ag) +4 | drpyz.(¥ — ay) V7,(r —a,) 


S far FAA) sre 2 2 [ ar Przo(t — a.) ; (3 8) 

a “Sy io; lr—a,| ’ pat 

where p,,, and p,z, are related to V,, and V,, by the Poisson equation 
V?V 7i(r — aj) = —4rpyz(r—a)), i=1, 2. (3.9) 


Using (3.9) is it easily shown that (3.8) reduces to the form 


ZL, Zs 
AU=— = 5 Vaal |a,—a,|)-— = 5 Vaal |a,—a,|) 


228s | drpyzi(" — a,)Vzo(r — ay) +3 i drpyzo(t—a,)Vz,(r—a,). (3.10) 


But from (2.3) we may write 


Viz (2k|r, — 
Pigi(l — a; ~ 5a [an aed (3.11) 


and it may then be shown tect that (3.5) has the form 


AT =~ drpya(r—a)V alta 


—t [ a puzal— ay) Veale —ay) Be EOS) 

Thus, adding (3.10) and (3.12) we are left with the final simple result 
Z Z 
AB=AT +AU =— = V 72(|a1— el) — 5 V p1(|a2— a))). (3.13) 


We now demonstrate directly the connection with the earlier treatment 
given by Alfred and March (1957). As pointed out in §3, the basic 
eqn. (3.1) of the first-order treatment reduces to the Mott eqn. (3.2) in 
the slowly varying potential limit, and the appropriate | one-centre 
solution is 


Z; 
Viilr)= — — exp (—9"). so ee Le ESS 
Thus, from (3.13) we then find 
rae ere mB STADE 2 oS 8 3.16) 
a 


which is just Lazarus’s formula. A numerical comparison of the predictions 
of (3.13) and (3.15) will be made later. 
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§4. Impurtry Dirrusion IN METALS 


In the Lazarus theory, as modified by Le Claire, the energy of movement 
is calculated using a model for the saddle-point configuration which was 
originated by Huntington (1942). The model consists of a linear three- 
centre configuration of point charges, —$, Z, —4, the half-vacancies being 
separated from the impurity by distances (11/16)a in the Huntington 
approximation (see Le Claire 1961 for detailed reasons for this choice). 
The M-centre calculation of Appendix 1 immediately gives the desired 
result, which, in Le Claire’s notation, may be written quantitatively as 


AH, AB = —4V-(11a/16) ¥Z Vogl lojra) 


Substituting V(r) as given by the slowly varying potential approximation 
in (4.1), we find 


AH,+AH= -— exp (— 1lqa/16). 7 (Sima Se) 


Z 
(11/16)a 
This differs from Le Claire’s result only by a factor a of the order of unity, 
a being Z dependent (see Alfred and March 1956). The introduction of a 
may be thought of as an attempt to introduce higher order perturbation 
corrections (see also the remarks of §6) and is necessary if curvature in 
the activation energy versus Z plot is to be explained. In our approxima- 
tion, « is simply unity. 

We record in the table results for divalent impurities in copper, as 
calculated from (3.13) and (3.15) for AH, and from (4.1) and (4.2) for 
AH,+AE. It will be seen that the results of our treatment provide a 
fair numerical justification of the formulae used previously. In view of 
this, and the crude approximation involved in the choice of the saddle- 
point model, it does not seem necessary at this stage to rework the 
predictions of the theory of impurity diffusion as presented by Le Claire 
(1961), although eventually numerical revision may become desirable. 


Results for various energy terms for divalent impurities in copper 

Impurity-vacancy interaction energy AE (a=4-8 atomic units) 
(1) From eqn. (3.13): AH = —0-0028 (—0-076 ev ) 

(2) From Lazarus formula (3.15): AH = —0-0021 (—0-057 ev) 

Hinergy AH,+AE, for saddle-point configuration described by Huntington’s 
model (see § 4) 

(1) From eqn. (4.1): AH,+ AH = —0-0095 (—0-26 ev). 
(2) From eqn. (4.2): AH,+ AH = —0-0128 (—0-35 ey). 


§5. Srapmiry or Di-vacaNcy : COMPARISON WITH THE CALCULATION 
OF SEEGER AND Bross 


The existence of the di-vacancy is important in a number of fields (see, 
for example, Lomer 1959). As Lomer points out, the stability of the 
di-vacancy does not seem established beyond doubt, although recent 
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work (Weizer and Girifalco 1960) suggests on both theoretical and 
experimental grounds that it should be stable. 

The only calculation which considers the electronic distribution round 
the di-vacancy, prior to our work, is that of Seeger and Bross (1956). 
However, their work on the di-vacancy is incomplete, in the same sense 
that the work of Fumi (1955) on the formation energy of a vacancy was 
incomplete. While some justification can be given for the partial 
cancellation of the terms in the energy omitted by Fumi (Lidiard and 
Tharmalingham, private communication) for the single-centre problem, no 
consideration has as yet been given to this point for the di-vaeancy. 

In order to throw light on this point, we shall now attempt to cast our 
treatment into a form which parallels the calculation of Seeger and Bross. 
We note first that these workers, following Fumi, focus attention only 
on changes in eigenvalue sums for model potentials chosen to satisfy the 
Friedel sum rule. It is clear then that their calculations are not complete, 
in that Coulomb corrections are not properly accounted for. In our work, 
‘on the other hand, we have used the full energy expressions within a Hartree 
framework. 

We consider then, at this stage, the change in eigenvalue sum using 
-our model of the di-vacancy. The potential energy term in the eigenvalue 
sum for a complex Z,Z, is 


| erence) Kode, Stee Te ate tae seme CORD 


and we now invoke the superposition theorem to obtain 


far [pr gi(t — ay) + pize(r — Ay) — po] [Vzi(r — a,) + Vzo(r — ag). (5.2) 


Again, separating out unperturbed, one-centre, and two-centre terms, 
the potential energy contribution to the change in eigenvalue sum is 
clearly 


AU eigenvalue = | dr|pyz,(" — a,) V z2(r — ag) + Przo(¥ — ag) Vaa(r — a4) ]. (5.3) 


From the Mott approximate solution (3.14) we obtain the result, for 
—qexp (—qa). io tel ee hee 4) 

(5.3) has been obtained numerically, as it seemed of interest to know how 

the energy change was built up from potential and kinetic contributions. 

We obtain 

AU eigenvalue = — 0:0074 atomic units = — 0-20 ev, 
whereas (5.4) yields — 0-26 ev. 
Also, from (3.12) it follows that 


AT = —4$AU eigenvalue = 9° 10 ev, 
-and hence 
AE eigenvalue = — 9° l0ev, 
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whereas using the approximate solution (3.14) the result is —0-13 ev. 
Since AZ eigenvalue is just the quantity which Seeger and Bross take as 
a measure of the binding energy, these calculations would then indicate 
stability of the di-vacancy, and our result of 0-1 ev is to be compared with 
their value of 0-3ev. Unfortunately, our estimate of the total energy 
change AZ is given by (3.13) and is AZ =0-076ev. This indicates repulsion 
and we see that in the present framework even the sign of the interaction 
is changed when we use the eigenvalue sum as a measure of the binding 
energy. Our work therefore must cast serious doubt on the validity of the 
Seeger—-Bross energy; as even if we adopt our Coulomb correction of 
0-18 ev and their eigenvalue sum, their binding energy is reduced to 0-12 ev. 


§6. HiagHER APPROXIMATIONS 


There is little doubt that, even within the framework of the present 
model, higher approximations should eventually be considered. Thus, 
the one-centre calculations of March and Murray (1960, 1961) show that 
the electron density given by the first-order eqn. (3.1), while containing 
all the essential features of the problem, is not capable of complete 
quantitative accuracy. The framework developed by these workers is 
general enough however so that the second-order corrections can be 
estimated in principle, and this is under discussion for the one-centre case 
(March and Murray, to be published). The calculations for the present 
problem are, unfortunately, very formidable indeed, the main obstacle 
being the loss of the superposition property when terms of 0(V?) are 
included in (3.1). Thus, entirely new calculations of the two-centre 
potential and density would be required. However, such a calculation 
was in fact carried through completely by Alfred and March (1957) for 
a divalent impurity and a vacancy in Cu, in the simpler case when the 
potential is assumed slowly varying, and this enables us to draw some 
tentative conclusions. 

Thus, the full computations of Alfred and March referred to above 
showed that while a result was obtained in quite reasonable numerical 
agreement with the first-order interaction energy, the inclusion of exchange 
was important. What emerged was that while the kinetic energy was 
given reasonably satisfactorily by the first-order treatment, the summation 
of the slowly varying potential series to infinite order in V led to a rather 
different potential energy, but the inclusion of exchange brought the result 
back to approximately the first-order value. Also, as pointed out by 
Alfred and March (1957), inclusion of exchange in first order completely 
changes the results, and appears to lead to unreasonable conclusions}. 
We infer then that, at least for the impurity-vacancy problem, the full 
calculation of Alfred and March is most probably much better than might 


ri Nevertheless, for completeness it seemed worthwhile to include in Appendix 
2 a brief derivation of the expression for the exchange energy using the density 
matrix formulation employed here. 
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have been supposed from the apparently severe approximation of slowly 
varying potentials. 


§ 7. CoNCLUSION 


Within a Hartree approximation and assuming the validity of low-order 
perturbation theory, we have shown how a rather general treatment of 
interaction energies between point charges in an electron gas may be 
developed. In particular, the so-called ‘electrostatic’ model of inter- 
actions between defects in metals follows as a direct consequence of our 
perturbation approach. However, in view of the long-range oscillations 
in the electrostatic potential around a defect, it follows that the 
interaction energy also oscillates with distance. 

It turns out that in the theory of impurity diffusion, proposed by 
Lazarus and recently modified by Le Claire, these oscillatory effects 
should not be very important, and in particular the major term, the energy 
of movement, should not be greatly affected. Thus, even though detailed 
numerical revision will eventually be called for, our work seems to 
substantiate further the current theory. On the other hand, we feel that 
stability of the di-vacancy cannot be explained in terms of the point 
charge model employed here. The work of Seeger and Bross (1956) has 
been shown to be open to criticism on the grounds that changes in sums of 
eigenvalues are not a proper measure of the interaction energy. Indeed, 
in our calculations, even the sign is changed by the terms Seeger and Bross 
omit. 

Undoubtedly, in the applications referred to above, reservations must 
be made concerning the limitations of the model. No doubt relaxation 
effects around vacancies and impurities will eventually have to be included. 
Within the uniform gas approximation, a crude account could perhaps 
be taken of nearest-neighbour displacements using a model like that of 
Lee and March (1960). 

Finally we might briefly refer to some recent measurements of Rayne 
(1961) on optical absorption in Cu-Zn and Cu-Ge alloys. These 
experiments may be thought of as measuring the change in the Fermi 
level with impurity concentration, and are in serious conflict with Friedel’s 
earlier theory (Friedel 1954). However, since the change in the Fermi 
level is dependent on the interaction between the impurities, that is 
on the overlap of the screening distributions, it seems very likely that the 
present treatment, with long-range oscillatory terms, will lead to results 
markedly different from those based on essentially exponential screening. 
It is hoped to report on a quantitative treatment of this problem at a 
later stage. 

ACKNOWLEDGMENTS 


It is a pleasure to thank Dr. A. M. Murray for valuable discussions and 
material help with some aspects of the work, and Dr. W. H. Young for 
reading a preliminary manuscript and, by means of numerous suggestions, 
considerably influencing the final presentation of the theory. 


1294 G. K. Corless and N. H. March on 


The work reported here was performed under contract with the United 
Kingdom Atomic Energy Authority, A.E.R.E., Harwell. One of us 
(N. H. M.) would particularly like to thank Mr. A. D. Le Claire and the 
other members of his group for the stimulus and encouragement which 
he derived from a brief visit as a Vacation Consultant in the summer of 
1959, during which period the need for the present investigation became 
clear. 


APPENDIX 1 


INTERACTION ENERGY IN MULTI-CENTRE PROBLEM 


We outline here the generalization of the result of §3 for two centres to. 
M-centres. We now write 


V(r) = — — - —— —...=Vo71+ Vogat---- (A 1.1) 
Obvious generalization of the method of §3 then gives for the potential 
energy 
PE = 1 | delpo+ pr. ae a ey ra ee 
ZG: 
+ Vogilr — a;)+ Voze(r—az)+...]+ > —=4 
iS5 |aj—a,| 


and again omitting one-centre and unperturbed terms we have for the 
potential energy contribution to the interaction energy 


AU=-1>4Z, dr C1zilt — 34) ry 44 = 
Peau ar Meena eruiae yy eee? 
Syd Blaise dr piz;V z;(r —a)). (tay AZ) 
VA) tAJ 


Similarly, we have for the kinetic energy contribution to the interaction 
energy 


=| | dr, dr, Pods ad¥ era apt Hlrs= ta) 
me (2k|r, —r])? 
and using the result that 


~ ae | an See ro|) 


Pizj— 
=r 
r,—r,| 


we find 
ar=—} y | ar sesPalr a. ee ee) 
CF) : 
Thus from (A 1.2) and (A.1.3) we have 
ed soak Aaa een Erbe 4) 
i4j 
which is the desired generalization of (3.13). 
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APPENDIX 2 
EXCHANGE ENERGY IN PERTURBATION APPROACH 
BASED ON THE Density Matrix 


In terms of the Dirac density matrix p(r,r j,k) we may write for the 


exchange energy 
; TER | | dear, FAGuipiales pee Ao. 1) 
[r—ro| 


We now insert the perturbation expansion (2.1) for p in (A 2.1) and we 


‘find 
Dah : 
= £| [ dear, [po” + 2popi + px? + 2pops | (A 2.2) 


[r—ro| 


where we retain only second-order terms in V. We may then show that 


a 1/3 
- tf Jarare es [a i(=] | dr[p(r,r)]*®, (A238) 
0 


which is the usual result for free electrons. For the next term in (A 2.2) 
we may write almost immediately 


— 3] | drdr, Oo = dr[p(r,r) —po(r, r)]+0(V2). (A2.4) 
[r—ro| rr 
There are three contributions to the term of 0(V?), as listed below : 
i ic V(r,)V(rg) . 
0) a4 | Jan: drs Te pdel?H rs rl), 


fi) al te Far, Meira 


167° [r—ro| [ry —rol|r—rs| 


(iii) = | | dear te a | i Gree Vicnie® 


i (k|ry—ro|+ cee r,|+4]r— ral) 
[ry —Frollre—ril|r — ral 
Tn the case when V is constant, it is of interest to note that the sum of the 
above second-order terms reduces to | 


yp 
- 55] a 


or in terms of the density difference 


1 3 1/3 
ped uals Amp 
;(=5) | (p—po) 


a result which is easily shown to follow from the Dirac expression (A 2.3) 
for the exchange energy. This, in fact, is also true for the contribution 
(A 2.4). 
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Quench Softening of Cu Whiskers 


By F. W. Scuaprnk and H. B. M. Wo.rtsrs 


Koninklijke/Shell-Laboratorium, Amsterdam, 
(Shell Internationale Research Maatschappij N.V.) 


[Received August 9, 1961] 


It has been shown recently by electron transmission microscopy that 
prismatic dislocation loops may be formed in f.c.c. metals by vacancy 
condensation after quenching (Hirsch and Silcox 1958). On subsequent 
ageing of the specimens the loops anneal out with an activation energy 
which for Al is approximately equal to the energy of self-diffusion (Sileox 
and Whelan 1960). 

We have examined the influence of quenching upon the elastic strength 
of Cu whiskers, as it may be expected that prismatic dislocation loops so 
produced act as dislocation sources resulting in a lowering of the high 
elastic strength of near dislocation-free Cu whiskers. 

Whiskers were prepared by hydrogen reduction of cuprous bromide, 
a method described by Brenner (1956). The main impurity in the 
specimens was silver in a concentration of about 30p.p.m. A special 
search was made for the presence of residual bromine in the lattice by 
means of radioactivation analysis. It was established that the bromine 
concentration of the whiskers was certainly less than 0-5 p.p.m. Crystals 
were selected for perfection under the microscope. Only those with 
apparently perfect side faces and a diameter not exceeding 8 » were 
used in the mechanical experiments. In general whiskers of larger 
diameter do not possess a high elastic strength (Brenner 1958). 

The specimens were annealed in a miniature furnace in high vacuum 
(10-§ mm of mercury) at 950°c and subsequently rapidly cooled by turning 
off the current. In this way a cooling rate of about 10° °c/s was achieved. 
In later experiments the quenching rate was further increased by causing 
the crystals to jump out of the furnace by means of electrostatic charging. 
This resulted in a quenching rate estimated to be 104 °c/s. Afterwards 
the specimens were mounted in a micro-tensile tester developed for 
straining whiskers at a rate of 1 mm/h (Wolters and Schapink 1961). 
About 30 whiskers of different orientations were tested in this way. 
It was found that the yield strength invariably dropped to 10-20% of 
its original value, the lowest values corresponding to the highest cooling 
rates. The decrease in strength is not due to a stress concentration at 


P.M. te 
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surface steps produced by evaporation during vacuum annealing, as 
slowly-cooled specimens did not show the effect. 

The diameter of the crystals ranged from 4 to 8p. As it was very 
difficult to handle still thinner whiskers of, say, less than 1 x, a possible 
dependence of this effect on diameter has not been checked. Such a 
dependence might be expected as during quenching vacancies may travel 
an appreciable distance and escape through the surface. For example, 
according to Kimura et al. (1959), the number of jumps of a vacancy 
in Au during quenching from 900°c is about 10° for a cooling rate of 
2-5x 104 °c/s. If the same is true for Cu, the vacancies move over a 
few microns and consequently the effect should decrease with decreasing 
whisker diameter. 

If the condensation of vacancies during quenching causes the observed 
decrease in strength, the crystals should recover on ageing. This was 
verified on specimens which were cut from the same whisker, each of 
which received a different ageing treatment in order to keep the inevitable 
scatter at a minimum. It was established in this way that recovery does 
occur, and that it is almost complete after 30 min annealing at 150°c. 
Preliminary experiments to determine the activation energy for the 
process gave a value of about 0-8ev. Although this is not an accurate 
determination, it is beyond doubt that it is much lower than the 
activation energy for self-diffusion in Cu, reported to be 2-05 ev (Kuper 
et al. 1954). It may be significant that the observed value is close to the 
energy of vacancy migration in Cu, which is about 0-7 ev (Kimura 
et al. 1959). 

It should be noted that the effect of quenching on the yield strength 
of whiskers is just the opposite of that observed in bulk crystals, where 
quench-hardening is known to occur (Cottrell 1957) and which is 
attributed to the blocking of existing dislocations by either single 
vacancies or vacancy clusters. 

A possible explanation of these results is furnished by the formation 
of Frank sessile dislocation rings originating from the collapse of vacancy 
discs. It was suggested by Read and Kuhlmann-Wilsdorf that such 
rings of Frank partial dislocation may be transformed into rings of total 
dislocation (Read 1953, Kuhlmann-Wilsdorf 1958). Such dislocation 
rings were in fact observed in quenched foils of Cu (Hirsch and Sileox 
1958), although their origin was somewhat uncertain. They may act as 
dislocation sources under influence of a shear stress which is of the order 
7~pb/r where b is the Burgers vector of the dislocation, the shear 
modulus and r the diameter of the loop. Inserting the values for Cu : 
w=4x 10" dyne/em?, b= 2-54 and r=250 A, we obtain 7~40 kg/mm?. 
This crude estimate is to be compared with the experimental values of 
6-12 kg/mm? for quenched specimens. 

The recovery of elastic strength should be governed by the rate of 
shrinking of the prismatic loops during ageing. The rate of shrinking 
of loops in Al has been observed to be determined by the activation energy 
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for self-diffusion (Silcox and Whelan 1960). Our value for the temperature 
dependence of the recovery process of yield strength in Cu whiskers is, 
however, much smaller than the activation energy for self-diffusion, which 
constitutes an obstacle to application of the ‘ vacancy-—cluster collapse ’ 
mechanism to the present results. 

It should be noted that an alternative explanation is possible, because 
it is not known for certain whether the excess vacancies condense on 
single close-packed planes. If, for instance, condensation should occur 
on several adjacent {111} planes, stable voids might be formed. The 
stress-concentration factor near the edge of such a void is given by : 
1+ 24/(t/s), where s is the semi-axis in the direction of the applied stress 
and ¢ the semi-major-axis of the hole (Neuber 1946). The creation of 
such voids would also explain the observed decrease in strength of the 
whiskers, if ¢/s reaches a value of about 10. The significance of such a 
model is that it makes the low activation energy for recovery plausible, 
because on annealing the voids probably become spherical, thus reducing 
their surface energy. The change in shape causes a decrease in stress- 
concentration factor and consequently hardening is observed on ageing. 
This process requires an energy not much larger than the energy for 
vacancy migration, in accordance with the ageing experiments. 

However, the final step in the recovery will perhaps be the disappearance 
of the spherical voids which is controlled by self-diffusion. The present 
experiments are not sufficiently accurate to determine whether a larger 
activation energy is associated with the later stages of recovery. 
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On the Creep of Metals at very Low Temperatures 


By P. FELTHAM 
Department of Metallurgy, The University, Leeds 


[Received June 28, 1961] 


Ix 1956 Glen described a rapidly decaying type of creep occurring in 
cadmium crystals at 1-5°K. His analysis of the process, based on Mott’s 
(1953) theory of low-temperature logarithmic creep, indicates that the 
energy barrier to the flow is of the order of 0-1 ev. Mott (1956) advanced 
the hypothesis that the barrier may be surmountable at such low tempera- 
tures because of the quantum mechanical tunnel effect ; Taubert (1959) 
believes however that the temperatures are still too high for quantum 
effects to assume significance in creep. The object of the present note 
is to show that a plausible explanation of the observations is possible also 
on the basis of classical dislocation theory. 

The figure shows part of a ‘ jogged ’ dislocation loop, with jogs indicated 
only in the predominantly ‘ screw ’ segment. Under a small shear stress 
the loop would expand locally, as indicated by shading ; a component 
of the line tension of the dislocation would then force the jogs aside. 
Their movement, being conservative, will be comparatively easy until 


M 


they reach the part of the dislocation which is of the mixed edge/screw 
type, at points such as M. Here they are stopped, for further movement 
could not take place without the generation of point defects, for which the 
effective stress is assumed too small. The displacements of the jogs are 
accompanied by slip arising from the increase of the area of the dislocation 
loop. Creep would thus occur at a rate determined by the migration 
velocity of the jogs ; the extent of the creep would be controlled by the 
mean free path of the jogs between barriers (e.g. at M), and would not be 
expected to be significantly temperature dependent, in agreement with 
experiment. 

An analogous creep mechanism, involving the conservative movement 
of jogs in edge dislocations, and leading to logarithmic creep, has been 
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described by Feltham (1961 a,b). The energy barrier in the latter 
process appears to be about 2:5k7'm, where T’'m is the melting point 
of the metal. In cadmium this would yield approximately 0-13 ev. 

The corresponding energy barrier for the conservative migration of 
jogs in screw dislocations would be expected to be rather less, as the 
local lattice distortion is smaller than near jogs in edge dislocations 
(Maddin and Cottrell 1955); but as the jogs are inflexible they would 
at least have to overcome the Peierls force. In the absence of other 
estimates we shall take the energy barrier in the case of screw dislocations 
equal to H=kTm, which for cadmium yields a value close to the latent 
heat of melting, i.e. approximately 0-065 ev. 

Now, if we write for the mean velocity of conservatively moving jogs 


v=vbexp {(—HA/kT)[1 —(c/o9)]}. 


where v is an atomic frequency equal to about 101! per second, b the 
Burgers vector and o, the stress at which the jogs would move without 
the aid of thermal activation, we obtain for the stress 


o= 0,1 —(kT/H) In (vb/2)]. 


Assuming that if creep is readily observable the jogs will migrate with a 
velocity of about 0-1 ./sec, then with b=3 x 10-§em and H=kTy, = 0-06 ev, 
the term in square brackets is zero at a temperature of 35°K. Hence, 
irrespective of the magnitude of o9, this process could no longer be rate 
controlling at higher temperatures ; above 35°K jogs in edge dislocations 
might control the creep rate. 

The results of Glen (1956) show that the energy barrier to transient 
creep in cadmium is in fact considerably higher at 77 and 90°K than at 
1—5°K, indicating a transition to a new energy barrier somewhere in the 
range 5-77°K, in agreement with the above inference. The preceding 
reasoning could, of course, be applied to creep at higher temperatures 
as well, and therefore leads to the conclusion that similar increases in 
the energy barriers should occur on raising the temperature further. 
Thus, apart from creep at very low temperatures, one might expect to 
observe energy barriers associated with 


(a) the conservative movement of jogs in edge dislocations, 

(6) vacancy formation by jogs in rapidly gliding dislocations, as 
discussed by Seeger (1955), 

(c) vacancy migration, giving rise to climb by jog displacement in 
dislocations of the edge or mixed types carrying a high concentra- 
tion of vacancies absorbed by elastic interaction (Friedel 1956), 

(qd), vacancy formation followed by migration, in a cooperative 
glide/climb process of mixed dislocations, in which the vacancies 
generated in the non-conservative displacement of the jog subse- 
quently lead to its migration resulting in climb, and 

(e) vacancy formation followed by vacancy migration, in the 
movement of jogged screw dislocations (Mott 1956). 
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Expressed as fractions of the activation energy of self-diffusion the 
barriers in close-packed metals would then be expected to be equal to 
about 0-13, 0-4, 0-6, 1-0 and 1-0 respectively. In high purity aluminium 
Sherby et al. (1957) established the existence of a series of activation 
energies in creep in the range 77—880°K ; again, in terms of the activation 
energy of self-diffusion these appeared to be equal to about 0-14, 0-5, 
0-75 and 1:0. The intermediate values may need correcting, as was 
pointed out by Basinski (1957) ; nevertheless agreement with the assumed 
values of the corresponding fractions may be regarded as being satis- 
factory. The strong elastic interaction between point defects and 
dislocations, particularly vacancies, may also give rise to creep barriers. 
Friedel (1956) estimates the interaction energy between vacancies and 
edge dislocations in copper to be equal to 0-35 ev ; the corresponding 
energy for interstitial atoms would be expected to be lower by an order 
of magnitude. Jogs, as well as interstitials formed in the early stages of 
creep, may thus control the creep rate at very low temperatures. 

Finally, it is of interest to note that (a) the to-and-fro movement of 
jogs on screw dislocations and (b) the pinning and un-pinning of edge 
dislocations at interstitial sites, may be responsible for the two ‘ Bordoni ’ 
peaks observed for example in the damping capacity versus temperature 
curves of copper subjected to oscillating stresses at very low temperatures. 
Both effects could account qualitatively for the observation that in slightly 
strained metals, in which jogs and interstitials could be assumed to be 
abundant, the damping is much more intense than in well-annealed ones 
(Niblett and Wilks 1960). 
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The Stacking-fault Energy of Graphite 


By C. Baker, Y. T. Cuou and A. Ketiy 
Department of Metallurgy, University of Cambridge 


[Received August 21, 1961] 


STACKING-FAULT energies can be obtained from electron transmission 
micrographs by two methods. The first, due to Whelan (1958), involves 
measurement of the radius of curvature of a partial dislocation at an 
extended node. The interaction of the various partial dislocations at the 
measuring point is neglected and the value to be taken for the line tension 
of the partial dislocation is uncertain, in cases where the radius of curvature 
of the partial dislocations is comparable with the thickness of the foil. 
The second method is to measure the equilibrium separation of the partial 
dislocations in along straight extended dislocation. This has the advantage 
that the solution of the elasticity equations is exact provided the thickness 
of the foil is much larger than the equilibrium separation of the partial 
dislocations. We have applied the second method to graphite using 
anisotropic elasticity. 

In the hexagonal lattice the elastic constants are invariant for a rotation 
about the hexagonal axis, so that for a few special problems the use of 
anisotropic elasticity becomes much simpler than for a cubic crystal. 
Using anisotropic elasticity and the Fourier synthesis method (Chou and 
Eshelby 1961) the width of a long straight extended dislocation with its 
line and Burgers vector both lying in the basal plane is given by: 


2 
W= FAC + 2 cos 2a) + K,(1—2 cos 2«)}. 
8iry 


Here bis the Burgers vector of the partial dislocations, wis the angle between 
the dislocation line and the total Burgers vector and y is the stacking-fault 


energy ; : 
Ks= {3Cu(Cu 7A C1,)}12, 


a C (Gc ie C3) as 

Ke=(Cy4+C {ee oe ; 
aii tanitiss) C3 (C3 + Cig + 2C ga) 
where Cie = (Ci, Gagy": 
Thus for a total dislocation in the edge orientation we have 
b2 
We = (3K,—Ks) 
Say 
and for a screw 


b2 
Ws= —— (3Ks—Ke). 
Sary 
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The best values of the elastic constants for graphite are those derived by 
Bowman and Krumhansl (1958) from spectra, compressibility and specific 
heat measurements. The values obtained are: 


C,,= 1113 10 dynes cm ; Ci, =0-28 % 10 dynesiem a2. 


C2 1°38 x10" dynesem -2; C= 2310 dynescma. 


The accuracy claimed for these figures is 20-50%. A value of Cj; is not 
obtained. In determining the stacking-fault energy of graphite we have, 
therefore, two unknowns: the stacking-fault energy and one of the elastic 
constants, i.e. Ke is unknown. Figure 1} shows a loop of extended dis- 
location lying in the basal plane of a graphite crystal. The partial dis- 
locations are clearly resolved and their Burgers vectors lie along ¢1010) 
directions: These loops of extended dislocation are found in material 
irradiated in BEPO to a dose of ~5x10!8n.v.t. The radiation induced 
defects have pinned dislocation sources operating in the crystal during 
cleavage. Figure 2 shows graphically how the separation of the partial 
dislocations varies around the loop. The maximum separation of 15004 
occurs when the total dislocation is in edge orientation and the minimum 
of 800 4 in the screw orientation. We have then, 


From this, and K,=3-12 x 104 dynescm~?, we find K,=4-24 x 10" dynes 
cm~*, whence C,,=7-7 x 10‘dynes em. 

Since K, and Ks are now known y is obtained from either W, or Ws, 
giving in both cases y=0-5lergscm-?. In the above calculations 
C33 = 1-8 x 10. dynes cm~? has been used. The value of y is subject to the 
error in the elastic constants and there may also be a systematic error in 
the measurements of W, and Ws due to the fact that the image of a dis- 
location does not coincide with the position of the dislocation. There may 
also be an error due to the fact that in graphite the thickness of the foil is 
not many times greater than the separation of the partial dislocations. 

Preliminary estimates of the stacking-fault energy of graphite have been 
given, Williamson (1960), Amelinckx and Delavignette (1960) using either 
an assumed value of the shear modulus or treating C,, as the relevant shear 
modulus and assuming isotropic elasticity. Siems etal. (1961) haverecently 
submitted for publication a comprehensive treatment of stacking-fault 
energies in layer structures also using anisotropic elasticity. Their treat- 
ment for graphite is slightly different from that used here in that Oj, is 
not obtained explicitly, an effective value of Poisson’s ratio being derived 
instead. The value obtained for y is 0-58 erg cm~* in excellent agreement 
with that obtained here. 

The values obtained using anisotropic elasticity differ by more than an 
order of magnitude from those obtained using isotropic elasticity. For 


ne SS eee 
{ Figure | is shown as a plate. 
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instance the radius of curvature of partial dislocations at a node in graphite 
is 0-95. Using Whelan’s formula 


y= Gb?/2R 


with G = 2-3 x 10!°dynes cm~ gives y=0-025ergsem-2. Ifthe equilibrium 
separation of partial dislocations in a long straight extended dislocation is 
used to find y employing the formula 


_ Gb? (2-7 1 2v : 
eee Sa Ny [ = 5 00820 |, 


12) 60 120 180 240 300 360 


—> Angle® 


Variation of the equilibrium partial dislocation separation with the angle to the 
total Burgers vector (angle not equal to «) for a dislocation loop. 
Maximum separation 1500 4 in edge orientation and minimum 800 A in 
screw orientation. 


where v= Poisson’s ratio, then one finds from the value of We and W; given 
above that v=0-265. If G is assumed to be 2-3 x 10!°dynescm~, then 
y=0-038ergsem~*. It is quite clear that in graphite and possibly in many 
other materials which are elastically anisotropic that the anisotropy must 
be considered in deriving values of the stacking-fault energy. All the 
stacking-fault energies reported here have been determined on irradiated 
crystals. Irradiation to these doses produces no measurable change in 
the radius of curvature observed at an extended node and hence no 
measurable change in stacking-fault energy. 
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REVIEWS OF BOOKS 


Progress in Low Temperature Physics. Vol. III. Edited by C. J. Gorrer. 
(Amsterdam: North-Holland Publishing Co., 1961.) [Pp. xiii+495.] 
Price £4 16s. Od. 

In his preface Professor Gorter tells of the doubts that have assailed him, 
whether there was a continuing need for reviews of the progress in low tem- 
perature physics. I have shared and given expression to these doubts, but am 
for the moment silenced by his latest production. If he can continue to give 
us such magnificent value for our money, may he long continue to edit: this 
series. What is characteristic of the present volume is that the topics dis- 
cussed, though having their origin in the realm of low temperatures, involve 
problems that are the very stuff of physics. Vinen on vortex lines in Helium IT, 
Grilly and Hammel on ?He, Bardeen and Schrieffer on superconductivity, all 
these illuminate the central theme of the quantum mechanical many-body 
problem; Azbel’ and Lifshitz, on electron resonances in metals, are concerned 
with metals and not just with metals at low temperatures; and one could 
extend the list to cover almost every one of the eleven articles. Here indeed 
we see low temperature physics in its proper perspective as part of a much 
more important whole, and Professor Gorter is to be warmly thanked for 
revealing the truth so clearly, while providing a series of reviews that will have 
an enduring worth. AS BoE: 


Mechanical Radiation. By Rosert Bruce Lrypsay. (McGraw-Hill Book 
Co., 1960.) International Series in Pure and Applied Physics. [Pp.ix+415.] 
Price £3 17s. 6d. 

THIS is a substantial textbook on the theory of waves and radiation, of a 

mechanical nature, in solids, fluids and gases, designed for advanced under- 

graduate work or first-year post-graduate work at American universities. It 
is well written and, although comprehensive in scope, deals with a great variety 
of topics in detail. There are copious references, many of them to recent work, 
both theoretical and experimental. The treatment is theoretical, but the point 
of view is a thoroughly physical one. There are three introductory chapters 
on the general theory of waves and vibrations. The general theory is then 
applied to a variety of subjects, in chapters dealing with waves in strings, 
membranes, elastic media, rods and plates, incompressible fluids, sound waves 
in fluids, elastic waves in the atmosphere, the sea, and the earth, and acoustic 
radiation in closed spaces. As may perhaps be expected in view of the author’s 
own work, the chapter on sound waves in fluids is the most substantial of these. 

The book concludes with a useful account of the relation of mechanical radiation 

and the properties of matter. EO GE. 


Proceedings of the Fourth Conference on Carbon held at the U niversity of Buffalo, 
New York. (Pergamon Press, 1960.) [Pp. 778.] Price £7 10s. Od. 
Tux Biennial Conferences on Carbon have since their inception in 1953 provided 
a meeting place for industrial and academic scientists interested in carbon in all 
its aspects. The Proceedings of the Fourth Conference, held in 1959, contain 
over 70 papers from laboratories in many countries. These cover a wide 
variety of topics ranging from theoretical papers on electronic structure, energy 
levels and Fermi surfaces to the technology of carbon brushes and fillers. The 
book probably provides the most comprehensive available survey of current 
researches on carbon and as such is invaluable to anyone concerned with the 
detailed properties of carbon. aN ‘ 
The book is superbly produced by Pergamon Press, and is unusual in including 
colour reproductions of photomicrographs taken in polarized light. 
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Crystal Structure Analysis. By Martin J. Burercer. (New York and London: 
John Wiley & Sons, Inc., 1960.) [Pp. xvii+668.] Price £7 8s. Od. 


Tue author, who is Director of the School for Advanced Study at M.I.T. and 
Professor of Mineralogy and Crystallography there, is well known for his earlier 
textbooks on structural crystallography and for the fact that he is co-editor of 
the Zeitschrift fiir Kristallographie. He has always had the direction and 
supervision of the work of many research students and others engaged in 
learning the methods of crystal-structure analysis and the present book has 
grown out of the notes that he prepared to help them. 

As would be expected, the style is excellent. One can only regret that the 
price may prevent purchase of this book by the students who would most 
benefit from having a personal copy. This is all the more reason why it should 
be available in all scientific and technical libraries. Some of the material given 
here is to be found equally well treated in other textbooks but attention may 
be drawn especially to the discussion of twinning in relation to the collection 
of data for-erystal structure analysis, and to the detailed description of intensity 
measurement. Throughout the book there is a precise and detailed account 
of recommended techniques, together with many useful tables and informative 
diagrams. 

Subjects not treated in detail are the Patterson synthesis and vector space 
(which are the subject of a companion volume published in 1959), reflection 
statistics, optical analogue methods of analysis and the determination and 
study of anisotropic thermal motion, bonding anisotropy and disorder. K. L. 


X-Ray Microscopy. By V. E. Cosstert and W. C. Nrxon. (Cambridge: 
University Press, 1960.) Cambridge Monographs on Physics. [Pp. 406.] 
Price £4 Os. Od. 


For years Cosslett, with his collaborators, has excited considerable interest in 
the development of x-ray microscopes and we welcome indeed this first com- 
prehensive account by Cosslett and Nixon of the use of x-rays for microscopical 
investigations with both low and high magnification. Both have contributed 
much in this field. This is a technique, which 20 years ago appeared to offer 
very real possibilities, but it has taken many years to overcome the formidable 
difficulties encountered in connection with focusing x-rays. Resolutions 
obtained are perhaps only twice that of the optical microscope, but of course 
with x-rays internal x-ray photographs are obtained with all the fascinating 
additional information these offer, especially for opaque materials. Then 
there are X-ray emission analysis, micro-diffraction and stereo-microscopy 
and quantitative analysis as additional attractive features. The book is very 
comprehensive indeed. First, contact microradiography is examined in detail 
and its ultimate resolution discussed carefully. Then x-ray microscopy by 
point projection (to which the authors have contributed much), is examined. 
This is followed by a chapter on reflexion methods and mirror systems, a section 
which is very well illustrated. Reflexion from curved crystals follows. 
Chapters follow on the biologically important methods of x-ray absorption and 
emission micro-analysis, and here we have chapters packed with practical 
‘know-how’ of immense value to the experimenter. Chapters 8, 9, 10, 11 
are all also mines of know-how. Eight gives an account of methods of exciting 
X-rays and X-ray sources and their characteristics, 9 describes specimen pre- 
paration techniques, 10 is devoted to techniques of contact microradiography 
and 11 to techniques of projection microscopy. Applications to biology and 
to metallurgy receive their full treatments and a final chapter on very recent 
techniques is included. Some valuable tables on absorption and emission 
data are added in appendices. A close packed bibliography fills no less than 
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23 pages. The index is rather brief. Although the book is naturally expensive 
it is generally admirably produced and certainly not over-priced. There is one 
weakness in production, for although there are 32 pages of delightful plates, 
these are all fastened together in one group in the centre of the book. How 
much better it would have been to have distributed these plates to be near the 
pages of their corresponding text. 

We are convinced this will remain an authoritative text on the subject for 
quite a long time. Sarl 


Thermodynamics. By P. T. LanpsperG. (Interscience Publishers, 1961.) 

[Pp. viii+499.] Price £5 9s. Od. 

THE author of this textbook has for some time made notable contributions to 
the more abstract aspects of thermodynamics; as befits a professor of applied 
mathematics he is here mainly concerned with these aspects and with their 
mathematical formulation. To quote from the author’s preface: ‘‘ The abstract 
nature of the subject is always allowed to become fully evident in this book ”’, 
“ General principles rather than detailed applications or physical mechanisms 
are emphasized’, and “ Questions of historical detail and of experimental 
results are omitted ”’. 

Although the author claims that “ this text should be suitable for anyone who 
is not deterred by abstract thought ”’, this reviewer feels that only the exception- 
ally brave physicist will venture to use this volume. If he does, he will be 
rewarded as the reasons for the rather general fear of physicists for thermo- 
dynamics will become clearer to him, and his understanding of this extremely 
difficult subject should greatly improve. 

Apart from the basic laws of thermodynamics this text also gives a few 
basic applications and a thorough discussion of the relation between thermo- 
dynamics and statistical mechanics. The text is generously interlaced with 
problems of a great variety. D. TER Haar. 


Close Binary Systems. By ZDENEK Kopau. (London: Chapman & Hall, 
1959.) [Pp. xiv+558.] Price £5 5s. Od. 


In this comprehensive book on a very fascinating astrophysical domain the 
author provides a well-written, informative and in most respects exhaustive 
treatment. There could hardly have been anyone better qualified to write it 
than Professor Kopal, who could draw on his wide personal and pioneering 
experience in a field to which he has contributed since the 1930’s. Six main 
sections illuminate the subject from many angles: we start with a development 
of a celestial mechanics of two closely associated stars, and follow the effects of 
various dynamical properties of bodies distorted by axial rotation and mutual 
tidal action on the observable orbit about their common centre of gravity. 
The evaluation of the potential and kinetic energy terms leading to the 
Lagrangian equations of motion, internal structure, precession and nutation of 
fluid components, etc., are some of the items dealt with. Observational 
features are introduced in the subsequent two sections on the theory of the 
light-curves and of the radial-velocity changes. The treatment is aimed to be 
as complete as possible to first-order terms, but we find frequent suggestions 
about the fruitful possibilities of an extension to higher orders. More than 
170 pages are devoted to the determination of the geometrical orbits, derived _ 
from the preceding analysis. Here, the progress made within the last two 
decades is impressive: the previous semi-empirical approach is now replaced by 
the introduction of exact-mathematical functions. However, in spite of the 
complexity of this new apparatus, problems like those of extremely close dis- 
torted systems are still beyond our full grasp. The author asks the question 
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‘“ where do we go from here? ” and indicates that the future might perhaps see 
the replacement of the present brightness-time analysis of the light-changes by 
that of their Fourier transforms, i.e. by the study of their frequency-spectrum. 
In conclusion, we should mention that the importance of these stars for our 
general ideas on stellar evolution is emphasized by numerous interesting 
interpretations. A. B: 


Techniques of High Energy Physics. By D. M. Ritson. (Interscience Pub- 
lishers, 1961.) [Pp. 540.] Price £6 6s. Od. 


Here’s how to do it: distillate of practical physicist, a careful comprehensive, 
and for the most part accurate survey of the present techniques for detecting 
and measuring elementary particles. Bubble chambers, cloud chambers, 
emulsions, ionization counters, scintillation and Cerenkov counters are all 
exhaustively discussed. If you want to design a beam, write a computer 
programme, or use a transistor, you can dip into this collection and find the 
snags, the-tricks, and the design factors all carefully laid out. Once or twice 
the pale ghost of a particle appears. 

This collection of essays is clearly intended as a compact reference work for 
research workers, and as such may well leave the ordinary reader rather breath- 
less. However, worthy of special mention are a delightfully written chapter on 
Scintillation and Cerenkov counters by Ritson and Weinstein, and a thoughtful 
and. critical little monograph on Emulsions by Widgoff. The first chapter is a 
valuable summary of formulae concerning the general properties and kine- 
matics of particles and y-rays, although a number of printing errors and occa- 
sional textual errors are confusing. Reading about autocode in the chapter on 
digital computers only strengthens my belief that it is easier to programme in the 
basic machine code. However, it is a pleasure to read that the design of mag- 
nets, like that of other experimental apparatus, leaves room for personal 
preference. Da Vas: 


BOOK NOTICES 


Noise in Electronic Devices. Papers based on material presented at a conference 
held by the Electronics Group of The Institute of Physics at the Services 
Electronics Research Laboratory, Baldock, Hertfordshire on October 2-3, 
1959. (London: Chapman & Hall, 1961.) [Pp. 100.] Price £1 15s. Od. 


Classical Mechanics. Second edition. By H. C. Corpen and Paine STEHLE. 
ee York and London: John Wiley & Sons, 1961.) [Pp. 389.] Price 
8. Od. 


Physical Chemistry. Second revised edition. By E. A. Momtwyn-Hucuss. 
(Oxford: Pergamon Press, 1961.) [Pp. 1333.] Price £4 4s. Od. 


Fast Reactors. By R. G. Patmer and A. Piatt. Nuclear Engineering Mono- 
graphs No. 12. (London: Temple Press, 1961.) [Pp. 93.] Price 12s. 6d. | 


Turbulence: Classic Papers on Statistical Theory. Edited by S. K. Frrep- 
LANDER and L. Topper. (New York and London: Interscience Publishers, 
1961.) [Pp.187.] Price £2 5s. 0d. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. ] 


A. HOWIE and P. R. SWANN Phil. Mag. Ser. 8, Vol. 6, Pl. 147. 


Extended and contracted nodes in a Cu+8%, Al alloy (Rx3000 4). One of 
the nodes lies on a different plane from the others. 


Fig. 6 


Nodes in Ag +2°%, Zn alloy (R300 A). 


S. HOWE and C. ELBAUM Phil. Mag. Ser. 8, Vol. 6, Pl. 148. 
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X-ray micrograph of zone-refined Al crystal. Axial orientation parallel to 
<110> direction. Growth rate: 1:3 mm/min. Reflecting plate: (111). 
Magnification: 24x. In the thin region (right of photograph), note 
arrows pointing to scattered, rather long dislocations in interior of 
crystal. The thick region (left of photograph) contains many short 
dislocations. Note great number of dislocations lying near surface of 


crystal. 
] ( 
Fig. 9 


micrograph of zone-refined Al crystal. Axial orientation parallel to 
<111> direction. Growth rate: 2mm/min. Reflecting plane: (220). 


Magnification: 47x. 


X-ray 
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X-ray micrograph of crystal shown in fig. 9 after electropolishing. Arrow 1 
indicates isolated dislocations. Arrow 2 indicates scratch in emulsion. 
Note that region to right of arrow 1 is free of visible dislocations. 


Magnification: 28 x. 
Fig. 11 


X-ray micrograph of 99-999% pure Al crystal. Axial orientation parallel to 
<110> direction. Growth rate: 4 mm/min. Reflecting plane: (111). 


Magnification: 22x. 
Fig. 12 
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X-ray micrograph of crystal in fig. 11 after electropolishing. Arrow 1 indicates 
the same dislocation-free area in both pictures. Arrow 2 in fig. 11 
indicates dislocation near surface that disappeared after electropolishing. 
Magnification; 22x , 


GHBAKER, vat. CHOU and A. KELLY: Phil. Mag. Ser. 8, Vol. 6, Pl. 149. 


Transmission micrograph of graphite, irradiated to a dose of ~5 x 1018 n.v.t., 
showing a dislocation loop with its Burgers vector in the basal plane. 
The partial dislocation separation in edge orientation of the total dis- 
location is nearly double that in screw orientation. 


